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The book is a monograph written as a result of research by the 
author. The diffraction of plane electromagnetic waves by ideally 
conducting bodies, the surface of which have discontinuities, 1s in- 
vestizated in the book. The linear dimensions of the bodies are 
assumed to be larrve in comparison with the wavelength. The method 
developed in the book takes into account the perturbation of the field 
in the vicinity of the surface discontinuity and allows one to sub- 
stantially refine the approximations of geometric and physical optics. 
Expressions are found for the fringing field in the distant zone. 

A numerical calculation is performed of the scattering characteristics, 
and a comparison is made with the results of rigorous theory and with 


experiments. 


The book is intended fer physicists and radio engineers who are 
interested in diffraction phenomena, and also for students of advanced 
courses and aspirants who are specializing in antennas and the 
propagation of radio waves. 
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FOREWORD 


First of all, one should explain the term “physical theory cf 
aiffraction". In order to do this, let us discuss brlefly the histo- 
rical development of diffraction theory. 


If one investigates, for example, the incidence of a plane elec- 
tromagnetic wave on a body which conducts well, all the dimensions of 
which are large in comparison with the wavelength, then the simplest 
Solution of this problem may be obtained by means of geometric optics. 
It is known that in a number of cases one must add to geometric optics 
the laws of physical optics which are connected with the names of 
Huygens, Fresnel, Kirchhoff and Cotler. Physical optics uses, together 
with the field equations, She assumption that in the vicinity of a 
reflecting body geometric optics is valid. 


At the start of the Twentieth Century, a new division of mathe- 
matical physics appeared — the mathematical theory of diffraction. 
Using it, rigorous solutions to the problem of diffraction by a wedge, 
sphere, and infinite cylinder were obtained. Subsequently, other 
rigorous solutions were added; however, the total number of solutions 
was relatively small. For sufficiently short waves (in comparison 
with the dimensions of the body or other characteristic distances) 
these solutions, as a rule, are ineffective. Here the direct 
numerical methods also are unsuitable. 


Hence, an interest arose in approximation (asymptotic) methods 
which would allow one to investigate the diffraction of sufficiently 
short waves by various bodies, and would lead to more precise and 
reliable quantitative results than does geometric or physical optics. 
Obviously, these methods must in some way be considered the most 
important results extracted from the mathematical theory of diffraction. 
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In the "geometric theory of diffraction" proposed by Keller, the 
results obtained In the mathematical theory of diffraction of short 
waves were exactly the ones which were used and gencralized. Here, 
the concept of diffraction rays advanced to the forefront. This 
concept was expressed rather as a physical hypothesis and was nct 
suiltatle for representing the fleld in all of space: it was not 
usable where the formation of the diffraction field takes place (at 
the caustic, at the boundary of light and shadow, etc.). Here it is 
impocsible te talk about rays, and one must use a wave Interpretation. 


What hac been said above makes it clear why a large number of 
works appeared in which the diffraction of short waves was investi- 
rated by other methods. Among those applied to reflecting bodies with 
abrupt surface discontinuities or with sharp edges (strip, disk, 
finite cylinder or cone, etc.) one should first of all mention the 
works of P. Ya. Ufimtsev. These works began to appear in print in 
1957, and it is on the basis cf them that this book was written. 

\ 

P. Ya. Ufimtsev studied the scattering characteristics by such 
bodies by taking into account, besides the currents being excited on 
the surface of the body according to the laws of geometric optics 
(the "uniform part of the current" according to his terminology), the 
additional currents arising in the vicinity of the edges or borders 
which have the character of edge waves and rapidly attenuate with in- 
creasing distance from the edge or border (the “nonuniform part of the 
current"), Cne may Zind the vadiation fleld created by the additional 
currents by comparing the edge or border with the edge of an infinite 
wedge or the border of a half-plane. In certain cases, one is 
obliged to consider the diffraction interaction of the various edges 
— that is, the tact that the wave created by one edge and propagated 
past another edre is diffracted by it (secondary diffraction). 


Such an approach to the diffraction of short waves has rreat 
physical visuallizability and allows one to obtain rather simple 
approximation expressions for the field scattered by various metal 


bodies. This approach may be called the physical theory of diffraction. 
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This name is applied to many works on the diffraction of short waves 
{n which the mathematical difficulties are bypassed by means of physi- 
cal considerations. 


It 1s clear that the physical theory of diffraction 1s a step 
forward in comparison with physical optics, which in general neglects 
the additional (edge) currents. The results obtained in this book 
show that with a given wavelength the physical theory of diffraction 
gives a better precision than physical optics, and with a given pre- 
cision the physical theory of diffraction allows one to advance into 
the lonrer wave region and, in particular, to obtain a number of 
results which are of interest for radar where the ratios of the dimen- 
sions of the bodies to the wavelength do not reach such large 
values as in optics. 


In addition, the physical theory of diffraction encompasses a 
number of interesting phenomena which are entirely foreign to physical 
optics. Thus, in a number of cases the additional currents give, not 
a small correction to the radiation field, but the main contribution 
to this field (see especially Chapters IV and V). If a plane wave 
is diffracted by a thin straight wire (a passive vibrator), then the 
additional current falls off very slowly as one goes further from the 
end of the wire. Therefore, the solution is obtained by summinz the 
entire array of diffraction waves (secondary, tertiary, etc.) which 
successively arise as a consequence of the reflection of the currents 
from the ends of the wires. It has a resonance character. Thus, the 
problem of the scattering of the plane wave by a finite length wire 
which is a diffraction problem of a slightly unusual type is solved 
in Chapter VII. The resulting solution is applicable under the condi- 
tion that the diameter of the wire is small in comparison with the 
wavelencth and length of the wire, and the ratio of the length of the 
wire to the wavelength is arbitrary. 


The final equations which are derived in this book and are used 
for calculations are not asymptotic in the strict sense of the word. 
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“herpofore, Lt is matural to pose the question: in what way will the 
subsequent asymptotle equations differ from them when at last one 
ehtains them in v.e mathematical theory of diffraction? One can say 
beforehand tha: the main term of the asymptotic expansicn will not, 

in tne general case, agree with the solution obtained or the basic of 
physical considerations: other (as a rule more complicated) slowly 
varying functions which determine the decay of the fields and currents 
as one coed’ further from the edzes and-borders, and also the diffrac- 
tion interaction cf the edges and the shadowing of tne edge waves will 
fisure in the main term. However, the refinement of the slowly vary- 
‘ine functions in the expression for the diffraction field is not able 
to seriously influence the quantitative relationships. This is seen 
from a comparison of the results obtained in this book with calcula- 
tions based on rigorous theory and other approximation equations, and 
also with the results of measurements. 


The rélationships obtained in this book also should help the 
development of asymptotic methods in the mathematical theory of dif- 


raction, since they suggest the character of the approximaticns and 
the structure of the desired solution. 


L. A. Vaynshteyn 


SLE Ged Sayles t ae 


INTRODUCTION 


In recent years, there has been a neticeable increase of interest: 
in the diffraction of electromasnetic waves by metal bodies cf compiex 
Shape. Such diffraction problems with a rigorous methematical formu- 
lation reduce to an interpretation of the wave equation or Maxwell 
equations with consideration of the boundary conditions on the body's 
surface. However, one cannot succeed itn finding solutions in the 
case of actual btedles of a complicated configuration. This may be 
done only for bedies of the simplest reometric shape — such as an 
infinitely long cylinder, a sphere, a disk, etc. It turns out that 
the resulting solutions permit one to effectively calculate the dif- 
fraction field only under the condition that the wavelenrtn 1s larrer 
than, or comparable to, the finite dimensions of the tbody. In the 

"quasi-optical"case, when the wavelenzth is a great deal less than the 
dimensions of the body, the rigorous soiutions usually lose their 
practical value, and it is necessary to add to them laboricus and 
complicated asymptotic studies. Here, the numerical methods for the 
soiution of boundary value problems also become ineffective. There- 
fore, in the theory of diffraction the approximation methods which 
allow one to study the diffraction of sufficiently short waves by 
various bodles acquire sreat importance. 


The field scattered by a given body may be calculated apprexi- 
mately by means of geometric optics laws (the reflection equations, 
see,for example [1-3]), from the principles of Huygens-Fresnel an} 
from the equation. of Kirchhoff and Cotler [3-6]. 


The most common method of calculation in t’.e quasi-optie case 
is the principle of Huygens-Fresnel in the formulation of Kirchhoff 
and Cotler — the so-called physical optics approach. The essence 
of this method may be summarized as follows. 
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Let a plane electromagnetic wave fall on some ideally conducting, 
body which Is found in free space. In the physical optics approach, 
the surface current density which ts induced by this wave on the 
{rradiated part of the body's surface is (in the absolute system of 
units) equal to 


j=, £- [nl]. (A) 


where c is the speed of light in a vacuum, n is the external normal 
to the body's surface, Hy is the mapnetic field of the incident wave. 
On the darkened side of the body the surface current is assumed to be 
equal to zero (3° = 0). Equation (A) means that on each element. of 
the body's irradiated surface the same current is excited as on an 
Ideally conducting surface of infinite dimensions tanrent to this 
element. The scattered field created by the current (A) is then 
found by means of .Maxwell's equations. 


It is cbvious that in reality the current induced on the body's 
surface will differ (as a consequence of the curve of the surface) 
from the current 3°, The precise expression for the surface current 
density has the form 


sa HS = 
where jt is the surface density of the additional current which 
results from the curve of the surface. By the curve of the surface, 
we mean any of its deviations from an infinite plane (a smooth curve, 
a sharp bend, a bulge, a hole, etc.). If the be *s convex and 
smooth and its dimensions and radii of curvature are lnarre in compari- 
son with the wav>-length, then the additional current is concentrated 
mainly in the vicinity of the boundary between the 1lluminated and 
shadowed parts of the bocy's surface. But if the body has an edge, 
bend, or point, then the additional current also arises near them. 
The additional current density is comparable to the density a as a 
rule, only at distances of the order of a wavelength from the corre- 
opouding edee, bend, or point. Thus, if the body's dimensicns 
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slenificantly exceed the wavelength, the additional currents occupy 
a comparatively small part of its surface. 


Since the current excited by the plane wave on an ideally cun- 
ducting, surface lo distributed uniformly over it (the absolute magni- 
tude of its surface density is constant) then the vector 3° may be 
called the "uniform" part of the surface current. The additional 
current jt which is caused by the curve of the body's surface we will 
henceforth call the "nonuniform" part of the current. In the physical 
optics approach, only the uniform part of the current is considered. 
Therefore, it is no wonder that in a number of cases it gives unsatis- 
factory results. For a more precise calculation, it is necessary to 
also take into account the nonuniform part of the current. 


In this book, the results of the author relating to the approxi- 
mation solution of diffraction problems are discussed and systematized. 
Essentially, these results were briefly discussed in a number of 
papers [7-14]. Roughly at the same time, the works of other authors 
devoted to similar problems appeared. We will discuss them in more 
detail (in §25) after the reader becomes accustomed to the concepts 
being used in diffraction problems of this type. For the present, let 
us only note that in these works, as a rule, other methods are used. 


In the book, problems of the diffraction of plane electromarnetic 
waves by complex metal bodies, the surfaces of which have discontinuti- 
ties (edges), are investigated. The dimensions cf the bodies are 
assumed to be large in comparison with the wavelength, and their 
surface is assumed to be ideally conducting. 


Obviously, if the edges are sufficlently far from one another, 
then the current flowing on a small element of the body's surfacc in 
the vicinity of its discontinuity may be approximately considered to 
be the same as on a corresponding infinite dihedral angle (a wedre). 
In fact, in Chapter I it is shown (see also [5] §20) that the nonuni- 
form part of the current on a wedge has the character of an cdre 
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wave which rapidly decreases with the distance from the edge. 
Therefore, one may consider that the nonuniform part of the current 

is concentrated mainly in the vicinity of the discontinuity. By means 
of this physically obvious assumption, the field scattered by a strip 
(Chapter I), by a disk (Chapter II), by a finite length cylinder 
(Chapter III) and by certain other budies of rotation (Chapter IV) 

is calculated. 


For a more precise calculation, however, it is necessary to keep 
in mind that the actual current distribution in the vicinity of the 
body's edges differs from the current distribution near the edge of 
the wedge. Actually, the edge wave corresponding to the nonuniform 
part of the current, propagated along the body’s surface, reaches the 
adjacent edge and undergoes diffraction by it, exciting secondary edge 
waves. The latter in turn produce new edge waves, etc. If all the 
dimensions of the body are large in comparison with the wavelength, 
then as a rule it is sufficient to consider only the secondary dif- 
fraction. This phenomenon is siudied in Chapter V using the example 
of a strip and disk. 


In the case of a narrow cylindrical conductor of finite length, 
the edge waves of the current decrease very slowly with the distance 
from each end. Therefore, here it is impossible to limit oneself to 
a consideration only of secondary diffraction, and it is necessary to 
investigate the multiple diffraction of edge waves. Chapter VII is 
devoted to this problem. 


The uniform and nonuniform parts of the current are more than 
auxiliary concepts which are useful in solving dciffraction problems. 
In Chapter VI it is shown that one is able experimentally to separate 
from the total fringing field that part of it which is created by 
the nonuniform part of the current. There, it is also shown that the 
depolarization phenomenon of the reflected signal ifs caused only by 
the nonuniform part of the current. 


PTD=11C-2 3+259-71 xi 


Let us note the following feature of the method discussed in the 
book. A physical representation of the nonuniform part of the current 
is widely used in the book, but nowhere are its explicit mathematical 
expressions cited. This part of the current is generally not expressed 
in terms of well-known functions. Obviously a direct integration of 
the currents when calculating the fringing field is able to lead only 
to very complicated and immense equations. Therefore, we will find 
the fringing field created by the nonuniform part of the current on 
the basis of indirect considerations without direct integration of 
it (see especially Chapters I - IV). 


The method by which the diffraction problems are solved in this 
bock may be briefly summarized as follows. We will seek an approxi- 
mate solution of the diffraction problem for a certain body by first 
having studied diffraction by its separate geometric elements. For 
example, for a finite cylinder such elements are: the lateral surface 
as part of an infinite cylindrical surface, each base as part of a 
plane, each section of the base rim as the edge of a wedge (the curva- 
ture of the rim in the first approximation may be neglected). Having 
studied the diffraction by the separate elements of the body, we will 
obtain a representation of the nonuniform part of the current and of 
the field which is radiated by it. Then secondary, tertiary, etc. 
diffraction is stuuied — that 1s, the diffraction interaction of the 
various elements of the body is taken into account. 


This method appeals to physical considerations, not only when 
formulating the problem but also in its solution process, and in this 
way differs from the methods of the mathematical theory of diffraction. 


Therefore, such a method may be referred to as the physical theory 
of diffraction. 


A whole series of other diffraction studies which appeared in 
the last five to ten years also are able to relate to the physical 
theory of diffraction. The first work which contained the idea of 
the physical theory of diffraction 1s evidently the paper of 
Schwarzschild [15] which was published at the beginning of this 
century and was devoted to diffraction by a slit. 
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One should note that approximate solutions of diffraction 
problems would be impossible without the use of the results obtained 
in the mathematical theory of diffraction. In particular, the 
rigorous solution to the problem of diffraction by a wedge which is 
attributed to Sommerfeld [16] is widely used in this book. In 
Chapter I this solution is obtained by another method. The works of 
Fok [17, 18] served as the starting point for numerous studies on 
diffraction by smooth convex bodies. The rigorous solution of the 
problem of diffraction at the open end of a wave guide [19] revealed 
the mechanism for the formation of primary diffraction waves, and 
their shadowing by the opposite end of the wave guide. The rigorous 
theory as applied to a strip and disk allows us to examine the 
precision of the approximation theory (see Chapter V). 
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CHAPTER IT 


DIFFRACTION BY A WEDGE 


As was already said in the Introduction, the field scattered by 
a body may be investi;ated in the form of the sum of the fields being 
radiated by the uniform and nonuniform parts of the surface current. 
The uniform part of the current 1s completely determined by the geo- 
metry of the body and the magnetic fleld of the incident wave. The 
nonuniform part generally is unknown. However, one may arproximately 
assume that in the vicinity of the discontinuity of a convex surface 
{ct will be the same as on a corresponding wedge. Therefore, it is 
necessary for us to begin by studying the diffraction of a plane elec- 
tromagnetic wave by a wedre. This chapter will be devoted to this 
problem. First we will investifate the rirorous solution of this 
problem (§ 1 and 2). Then we will find its solution in the physical 
optics approach (§ 3). The difference of these solutions determines 
the field created by the nonuniform part of the current (§ 4). 


§ 1. The Rtporous Solution 


The rigorous solution to the problem of diffractton of a plane 
wave by a wedse was first obtained by Sommerfeld by the methed of 
branching wave functions [16]. Later, the diffraction cf cylindrical 
and spherical waves by a wedge also was studied. A rather extensive 
bibllopraphy on these problems may be found, for example, in the 
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Paper of Oberhettinger [20]. Since the problem of diffraction by a 
wedge lies at the base of our studies, we considered it advisable not 
only to present the results of its rigorous solution, but also to give 
them a new more graphic derivation. The idea for this derivation 
follows directly from the work of Sommerfeld. Sommerfeld found the 
solution to the problem in the form of a contour integral, and then 
he transformed it to a series. However, one may proceed in the oppo- 
site direction: first find the solution in the form of a series and 
then give its integral representation. Such a path seems to us more 
graphic, and is discussed in this section. The necessity for a 
detailed derivation 1s caused by the fact that the results of Sommer- 
feld [16] are not represented in a sufficiently clear form, which 
hinders their use. 


Let us assume there is in free space (a vacuum) an ideally con- 
ducting wedge and a cylindrical wave source Q parallel to its edge 
(Figure 1). Let us introduce the cylindrical coordinate system r, #9, 
Zz in such a way that the z axis coincides with the wedge edge, and 
the angle $¢ is measured from the irradiated surface. The external 
wedge angle will be designated by the letter a, so that O<p%a . The 
coordinates of the source Q we will designate by os o> 


Let us investigate two particular cases for the excitation of an 
electromagnetic field. In the first case, it is excited by a "fila- 
ment of electric current" 


j-== — iop,3(¢ — ry, 9 — 9%). (1.01) 
in the second case, it is excited by a "filament of magnetic current" 


[P= —iwm,5(r — rq, 9 — %)- (1.02) 


The quantities Po and mM, here designate, respectively, the electric 
and magnetic nonenes of the filament per unit length along the z axis, 
w is the cyclic frequency (o=ck =e), 8(f — fo, 9 -- Py) == 8 (F - - Pq) 5 Ir (@ —- 9)] 


is a two-dimensional delta function which satisfies the condition 
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Figure 1. The excitation of a wedge- 
shaped region by a linear source. 


- source; P = the observation point; 
L - the Integration contour in Equation (1.10). 


ffir —r. 9 — Po) rdrde =} 
with integration over the neighborhood of the point To» 9° 


Here and henceforth, we will use the absolute system of units 


(the Gauss system), and we will assume the dependence on time is tn 
the form e~™ , 


In the first case, the "electric" vector potential Ay satisfies 
the equation (see, for example, [4]) 


AA + kta — 1 jf (1.03) 


and the boundary condition 


Av = 0 with 9=0 and Puta, (1.04) 
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In the second case, the "magnetic" vector potential A satisfles the 


equation 
AAT k®A™ =a == i° (1.05) 
and the boundary condition 
a4" ‘ 
“dg 779 with g==0 and 9=a. (1.06) 


It is natural to seek the solution of the nonhomogeneous 
Equattons (1.03) and (1.05) in the form 


* ad, (&r) HO (hry) singy, gysin 4? WIth r< or 
v a,J, (kr) H'" (kr)sin ¥.99 SiN VP WIth! >So: (1.07) 
v 6,J, (kr) H" (Arq) COSY, Fg COS EP WIth PK fs 


NY 6, (Ry) Hi (kr) cos v.49) COS¥.9 WIth 1 >%e (1.08) 
* "6 e 


4.=s—. 


The products = 


J, (kr)siny,9 J, (kr) cosy, 
and 
WY (Ar) sin v,9 H"" (kr) cosy, (1.09) 


are the partial solutions of Equations (1.03) and (1.05) without the 
right-hand member which sattsfy the boundary conditions (1.04) and 
(1.06). The remaining factors entering into Equations (1.07) and (1.08) 
ensure the observance of the reciprocity principle and the continuity 


of the fleld on the are r 2 To: The Bessel functton q,, (kr) enters 


these equations when r < To because it remains finite when r + 0, 


>] 


and the Hankel function HY (ke) is taken when r > lq in order that the 


solution satisfies the radiation condition. 
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The coefficients a. and by may be determined by means of Green's 
theorem 


er =: [ auus, dS =: rdrde (2.10) 
t . 


for the contour L in the plane z = const which 1s shown in Firure 1. 
Here, the external normal to the contour L is designated by the letter 
n. Applying Equation (1.10) to the functions Ae and AD and performing 


the limiting transitions r, » r, andr, +r. in it, we obtain 


1 0 2 0 


ark 
) rile =i fae 2.) 89 — drat, 


4dr 
~“2| a! 7=4 tm, { 3¢¢— gu)rat. 


Since here the integration limits are arbitrary, it follows from the 
equality of the integrals that the integrands are equal: 


OAS oAG inkp, 
i or | pea ' tee : (? i Ye)» ( i ‘ 11 ) 
OA vat _4ckm, » 
| fet " ee : “te i (? my %)- ( 1.12 ) 


Now let us substitute Expressions (1.07) into Equality (1.11) and 
multiply both members of the latter by sinvy . Then tnterratine the 
resulting equality over ¢ in the limits from 0 to a, we find 
a,== 3 kp. (1.13) 
In a similar way, let us determine the coefficients 
by=ty i hy (1.24) 


where 


a= A eee (1.15) 
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Consequently, tne electric current filament excites, in the space 
outside the wedge, the field 


‘ on 
rece Hp, » HI (kero) J, (ke) sin v9 sin up 
am 


with r<%, 


ao 2 
é = ke “y) 4, (kr) HY (Ar) sin v.99 sin vp 
a= 
with r>fe, 
E,=£,=0, H=-1 rot E, (1.16) 


and the magnetic current filament excites, outside the wedge, the 
field 
4a? = 
i km, y oa (kr), (ht) 60S %695 C08 9 


an 
with °r<fos 
= . 


@ 
iA wm, y tad, (kta) Hy” (kr) C08 v4 605 %4P 
am@ 


with r>le ‘ 
1 
-H,=H,=0, E=—— rot. (1.17) 


Now using the asymptotic equation for the Hankel function when 
kg CO 121), we have 


WV OED mage 0.8) 


Then Expressions (1.16) and (1.17) in the region r < r. take the form 


8) 
ae 1 bYpH (br) Xx 


= =the 
x } e 7 °s , (er) sin v4.99 SiN v9, 


(equation continued on next page) 


FTD=HC--2 3-259-71 : 


iE kim! (he) X 


hed i i) 
x } ae ~ 4, (kr) cos »,9 COS v,9 
3:30 e 


or 

E, = int’ pH (kr) (0 (7, P— Pl — Ur, P+). 

H,==irk*m.H) (kro) {4 (7, 9 —%)-+ u(r. PP), (1.19) 
where 

— hed -1= % 
u(r, a=2y ee J, (kr) cos vp (1.20) 
s=4 . : 

(p=9=) > 


Let us note, furthermore, that in free space the field of the 
electric filament with a moment Po is determined by the relationship 


i] 
E, = Ik" pH (hr,), (1.21) 


and the field of the magnetic filament with the moment m, is determined 
by the relationship 


2 = izk’m,H, (kr). (1.22) 
Therefore, the expressions in front of the square brackets in Equations 
(1.19) may be regarded as the primary field of the filament — the 
cylindrical wave arriving at the wedge edge. Now removing the fila- 
ment of current to infinity (ro +m), let us proceed to the incident 
plane waves 


£,= Ee" cos (9 -— Fa) ' E,= E,=0 (1.23) 


and 
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H, == Hy,-e7 iar cos (c =te. H, Be Hy 0. i (1.24 ) 


The field arising with the diffraction of these waves ty the wedge 


will obviously have the component 
E,==E,,[alr, 9-- 9) —4(r, P-F PI (1.25) 
and 


H,=Halu(r, 9--%)4-4(r, 9+9,))- (1.26) 
Let us find the integral representation for the function u (r, 
¥). For this purpose, let us use the equation (see [16], p. 866) 


J, (er)= a fe ( Ar cos B+e, (1-3 )} df, 


(1.27) 


where the limits I - III mean that the integration contour goes from 
region I to region III (Figure 2). The cross-hatched sections in 

the plane of the complex variable B (8') shown in Fifure 2 are regions 
in which Imcosj>0 (Imcos? <0), Therefore, in the sections of the con- 
tour exterdins to infinity the integrand strives to sero, ensuring 

the convergence of the intesral. Substituting Expression (1.27) into 
Equation (1.20), we obtain 


u(r, y= 
wt 


—_! [ eitreoss * iv, P— 844) » git, (G—8-9) 
Be Pee +e |a8 


s=] 


After summing the infinite geometric progressions and replacing the 


variable 8 by 8' = 8B = m7, the function u (r, ) acauires the forn 
a(r, y= 
ue : a 
—iks cos § 
=-— le ee ef 
x | | 1 +e) <te-e | ee 
l—e —e 
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Fivure 2. The invegration contours Figure 3. The integration contour 
in the complex plane 8. in Equation (1.28). 


As a result, we obtain the well-known integral of Sommerfeld 


e— iar cas 8 


a ‘ 
ingen (1.28) 


u 


The integration contour C is shown in Figure 3 and consists of 
two infinite branches. Since the integrand expression has poles at 
the points Pxa=2em—9 (m=0, +1, + 2,...) , then for the values of $ corre- 
sponding to the space outside the wedge (0< <a} the function u (r, yp) 
may be represented (with t<a<2n, 0<9,<2 ) in the following way: 


u(r, =o, ten? with—ecics, 
u(r, oY) =o(r, ¥) with 2¥<2a—sz, 
a(r, 9) =2(r, 9) elke cos 2a—9) 

oa with 20a—"#<4<2a, (1.29) 


where 


i-o+8 
t—e 


o(r, | nat 
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or 


f are 1 elke cos ftp 
o(r, ¥):=-5— sin -— oe 
= . carrer (1.30) 


The integration contours D and Dy are shown, respectively, in Figures 
3 and 4. : 


With an arbitrary incidence of 4 plane wave on a wedge, one of 
two cases may occur: (1) the plane wave "illuminates" only one face of 
the wedge (0<%<a—«) , and (2) the plane wave "illuminates" both 
faces of the wedge (a—*<9%<27) . Let us write out in more detail the 
functions u (r, ¥) corresponding to these cases. In the case M<3—* 
(Figure 5), we have 


u(r, 9 — 9%) = ’ ’ 
=o(r, 9— elk con (=~ ze) 
pen Lene cuedee 
, 0 —_ 
main gt eytentennee | 


u(r, F— P= 
=o(r, 9— Pp) pe rere \oivinese-eeete 
Br, PAM) =I(r, P+) 
u(r, P— Po) == 9(7s P— Po) \ 


‘) <a, 
u(r, PP) ==9(F, PA Po) with §+%<7<e 


: (1.31) 


and in the case %2—-"*<%<e (Figure 6) we have 


u(r, 9--%) =. \ 
=o(r, oe. ew itrcoste— ee 
ae ine peeee: 
s=0(r, at r)pe mer : 
u(r. 9M 
oe = - ith 
s=0(r, P—- epee NE ae ae 
SC ate cone te ko ee 
u(r, ?—%)) = | - 
=0(r, 9—P~rpte ees | with 
alr, aaa [nee \ a8) 
=u(r, eb 9) bere | 3c 
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The direction $¢ =m - $5 corresponds to the ray reflected in a specular 
fashion from the first face (we will consider as the first face that 
face from which the angles are calculated), and the direction = 

=s 2a— z— 9 corresponds to the ray reflected specularly from the second 
face (Figure 5 and 6). The functions e~"™* describe plane waves of 
unit amplitude: e~""**"-W describes the incident wave, ei iirtr 
describes the wave reflected from the first face, and g~trcos(te—e—e) —— 
the wave reflected from the second face. 


Figure 5. Diffraction of a plane 
wave by a wedge. The plane 
wave irradiates only one face of 
the wedge. $4 is the angle of 


incidence. The line 9 = fT - $5 


is the boundary of the reflected 
plane wave, and the line 
ou + $5 is the boundary of 


the shadow. 


Figure 4. The integration contour 
in Equation (1.30). 


Figure 6. Diffraction by a wedre. 
The plane wave irradiates both 
faces. The line 9 = 2a - 7 = $5 


is the boundary of the plane wave 
reflected from the second face 
(6 2 a). 
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§ 2. Asymptotic Expressions 


The integral 


t a ¢ gltrcos lye 
o(r.$) =f ta f ee, 
fore meretH9 (2.01) 


in Equations (1.31) and (1.32) generally is not expressed in terms of 
well-known functions. However, when kr >> 1 it may be calculated 
approximately by the method of steepest descents [21]. In integral 
(2.01), changing for this purpose to a new integration variable 


sa fte™ sins, s*=i (1 —cosQ), 


we obtain 
sia —— 7) 
(7.9) = sae ia | N rortstiger =a 
(coe F en EEA) con (2.02) 
-—@ 
where 
n=, (2.03) 


It is not difficult to see that the point s 2 0 is a saddle point: 

as one goes further from it along the imaginary axis (Re s = 0) in 
the plane of the complex variable s, the function e-* most rapidly 
increases, and as one goes along the real axis (Im s = 0) it decreases 
most rapidly. Therefore, when kr >> 1 the main contribution to the 
integral (2.02) is given by the integrand in the section of the 
contour in the vicinity of the saddle point (s = 0). 


The method of steepest descents is carried out by expanding the 


enter’ 


integrand (except for the factor ) into a Taylor series in powers 


of s. This series is then integrated term by term. If the integrand 
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expansion converges only on part of the integration contour, the re- 
sultant series obtained after the integration will be semiconvergent 
(asymptotic). Limiting ourselves to the first term in it, we obtain: 


cos — — cos — Y 2ekr (2.04) 


The remaining terms of the asymptotic series have a value on the order 
1 
of an? and higher. 


Expression (2.04) is valid with the condition (cos 5 —cos-+.) / kr >1 
and describes that part of the diffraction field which has the charac- 
ter of cylindrical waves diverging from the wedge edge. With the 
incidence of the plane wave (1.23) on a wedge, the electric vector of 
which is parallel to the wedge edge, the cylindrical wave is determined 
in accordance with (1.25) and (2.04) by the equation 


£,=—H,=Eq-(0(r.9—%)—9(r, Poll = 


YT aes 
e 
= Bul Vie (2.05) 
where 
sin = 1 { 
j= ota con 2 — con PLES (2.06) 


When the wedge is excited by the plane wave (1.24), in which the mag- 
netic vector 1s parallel to the wedge edge, the cylindrical wave has 
the form 


H, = EL =H [0 (1.9 — me) +09 Fel = 


ie es (m ++) 
ee (2.07) 
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where 


sin — : n , 
g= ry \ot Fem cos = — con TEE (2.08) 


In the vicinity of the shadow boundary (¢ x T+ 5) and near the 
directions of the mirror-reflected rays (p=*—%, 9 2a—%—%) 
Expressions (2.04) = (2.08) are not valid, since the poles 


7 a+ 


sin “>, $= Vie sin (a= -*S 


= /Ze 


of the integrand in (2.02) are close to s = 0 and, consequently, its 
expansion in a Taylor series loses meaning. Physically, this result 
means that in the indicated region the diffraction wave does not re- 
duce to plane and cylindrical waves, but has a more complicated char- 
acter. An asymptotic representation of the function v (r, ¥) in this 
region was obtained in 1938 by Pauli [22]; here we will present the 
derivation of the first term of the asymptotic series obtained in [22]. 


Let us multiply and divide the integrand expression in Equation 
(2.02) by the quantity 


cos ¢-+ cost = i(s* —is’) (si=2cos*3-) (2.09) 


and let us expand into a Taylor series in powers of s the function 


cos $-f- cos § 


(woe Bt 


which no longer has a pole at the saddle point (s = 0) when p = ¢ + 
+ $5 =z 7. Limiting ourselves in this series to the first term, we 
obtain 


ere A seal “1 eons | en tre* 22 


(2.10) 
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The integral here may be represented in the form 


t- J Py t- J heed a 
fiat ds = ons ds J e tte) “ét. 
s*— i323 2 


Changing the order of integration here, we find 


Va thea? ° ie 
= Tea © “ os (2.11) 


and finally 


Vite |cos 4] (2.12) 


The next term of the asymptotic expansion for the function v (r, yp) 
has a value, whose order of magnitude depends on the observation 
direction: in the vicinity of the border of the plane waves (7=*=+9) 
its order of magnitude is Ti » but far from it the order of magnitude 
is 1/kr in comparison with the term written in (2.12). 


It is convenient to represent Expression (2.12) in the following 
form: 


in — cos aT 
2°" a 2 thr cose € x 
(7.4) =—— n + ye 
08 — — cos 
Var cos — 
> 4 e'“dg. 
wren $ 


(2.13) 
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Here the absolute value of the lower limit of the Fresnel integral 
always equals infinity, and its sign is determined by the sign of 
cos ~/2. Therefore, when passing through the boundary of the plane 
waves (W = $6 + $5 = 1) the lower limit changes sign and the Fresnel 
integral undergoes a finite discontinuity, ensuring at this boundary 
the continuity of the function u (r, p) and consequently of the 
diffraction field. Actually, by means of the well-known equation 


@ Sa a 
leg _¥« i< 
Js rae (2.14) 


it is not difficult to show that 


o(r, 2-+-0)=--, o(r,a—0)= ao (2.15) 


and consequently 
a(r,a0)=+e™, (2.16) 


In view of the asymptotic relationships 


eld, ee 7 id, tl" 
I= "vp fe q 375 (with p> 1) 
if a . (2.17) 


hy ¢ 
Pauli's Equation (2.13) is transformed with ¥ Br |cos 3]>1 to the 
Expression (2.04). As was already indicated above, it determines the 
cylindrical waves diverging from the wedge edge. 


By means of Equation (2.13), one may also calculate the field in 
the vicinity of the direction 6 = 2a - 7 - $6 — that is, near the 
boundary of the plane wave reflected from the face > = a; for this 
purpose, it is sufficient to replace ¢ by a - ¢ and to by a - o> 


It is also interesting to note that in the case of a half-plane 
(n = 2) Equation (2.13) gives the expression 
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~i> Vitr co < 
o(r, y) ==e aS x e“dq, 
| omm (2.18) 


which completely agrees with the rigorous solution. Actually, with 
a = 2m, when the wedge is transformed to a half-plane, integrai (1.30) 
equals 


i elds cost 
(7,9 =— ys eh 
aie a (2.19) 


and it may te reduced to a Fresnel integral. For this purpose, let 
us divide the contour Do into two parts by the point > = 0. Summing 
the integrals over these parts of the contour, we find that 


+ ie 


0(r,9) == — + J ae ene 
2 2 


con 


i 
Now changing to a new integration variable s= /2e * sin = and taking 


into account Equation (2.09), we obtain 


sf 


Ff i ie - +) ey 
v(r, p= —° % cos fi aid” (2.20) 


The integral here was already calculated by us. Turning to Equation 
(2.11), we arrive at Expression (2.18) which — together with rela- 
tionships (1.25), (1.26) and (1.31) — give us the rigorous solution 
to the problem of the diffraction of plar.e waves by in ideally 
conducting half-plane. 
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§ 3. The Physical Optics Approach 


In the physical optics approach, the fringing field is sought 
as the electro-magnetic field created by the uniform part of the sur~ 
face current 


=~ [nH,)- (3.01) 
Let us recall that here n designates the external normal to the body's 
surface, and Hy designates the magnetic vector of the incident wave. 
First let us investigate the case 0°.4,<a--7, when the incident plane 
wave irradiates only one face of the wedge (Pigure 5). 


From Equation (3.01) it follows that the density of the uniform 
part of the current being excited on the irradiated face by plane 
waves (1.23) and (1.24) has the following components, respectively 


f =~ Ey:-sin gee Meee, Pm fl a0 (3.02) 
and 
=a Hae, i, =: j}=0. (3.03) 


For the purpose of calculating the field radiated by this current, we 
will use the following integral representation of the Hankel function 
(see [16], p. 866 ) 


, Kix @ 
HM (p) = —- ef cos 12 as ahs ( e? ch tat 
—» 


iz 
-ly+io 


(0«::8< x). (3.04) 


Assuming here p = kd and changing to a new integration variable 
C = dosh t, we obtain 
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He? (haya [PS at, 
| ia (3.05) 


It is easy to show by means of Equations (3.02), (3.03) and (3.05) 
that the vector potential 


* git Wi pce 
A(x, »O)=: EB J? (2) d: j y weer ees GG 
, 7 gy Ver ee (3.06) 


has the components 
A= Eorsintel, As=A,=0, (3.07) 
if the wedge is excited by plane wave (1.23), and 
As=7-He-l, A= 4, =0, (3.08) 


if the wedge 1s excited by plane wave (1.24). Here, q, designates 
the integral 


4 ——_ 
I= as Hy (kA ye Bk (3.09) 


Let us transform it by using the relationship 


HY (& fd pai) + fe ae 
(o= /k*—w*, Imo>0, d>0), (3.10) - 


It is not difficult to establish the correctness of this relationship 
by verifying that it changes into Expression (3.04) with the substi- 
tution w= k sin t, v = k cos t and &pd?fzt=p. As a result 


ef (0 lula) 
h=7 { ty (2 cos gy, — «) ms (3.11) 
—o 
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where the integration contour passes above the pole w = k cos o> 


Let us note that integral (3.11) is a function of |y|, and let 
us change to polar coordinates according to tne equations 


X= PCOS 9, 
lyj=rsing with ¢<s, 


Ww=—rsinewith p>. (3.12) 
Furthermore, by carrying out the substitution 
w==—keost (o==ksin'), (3.13) 
we obtain 
h=m mt withecs, | 
? 
= a agape ®t (3.14) 


The integration contour F is shown in Figures 7a and 7b. In Figure 7a 
the cross-hatched areas indicate the sections in the plane of the com- 
plex variable £ in which Imcos(i—9>0 ; in Figure 7b, the cross-hatched 
areas indicate the sections where Incos(t-+.9)>0 . Now let us deform the 
contour F into the contour G, (G,) for the values ?7<t(? >), and let 
us change to a new integration variable 


C=—t—9 withe<z, 
=3§~— (2e—g)with p>. (3.15) 
As a result, we obtain the following expressions: 
Owithe>r—g, 
1 eftrcos ly \ 
=, j COS Fy + C03 (§ + 5) ar 2 gt ike cos (¢+e8) 
° koin x, 

(3.16) 


with 7 F—%, 


if ¢ < 7 and 
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Figure 7. The integration contour in 


Equations (3.14). 


Owithpceto, 


! n eftr cos lye 
1 ah la fe + Cosa — ) + 2 em ite cos(g — ey 
k sing, 
with @>a-+g,, 
(3.17) 
if¢>m. The integration contour Dy 1s shown in Figure 4, 
By means of Equation (3.06) and the equality 
&,=tkA, (3.18) 


let us find the field which is radiated by the uniform part of the 
current excited on the face @ = 0 by the plane wave (1.23) 
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OF (9 e)— EMH wIthOCpce—p,, 


OF (9, 9) with®—%<e<e, 

D; (P. %) with=<pcet a, 

oF (9, 9.) — 7 1h 608 F—60) with R+4-9< 9a, (3.19) 
21 


where 


e! 4rcos 6 dz 


ry _i.: 
Oy (# %) = 5, sin wf Cos ye; cose tq) (3.20) 


It is not dirficult to see that with a—z<9<« , when both faces of 
the wedge are illuminated, the field excited by the second face 

(¢ = a) will be described by the same equations if one replaces ¢ by 
a- $, and $5 by a- $5 in them. 


Adding the field being radiated by the uniform part of the 
current with the incident plane wave (1.23), we obtain the diffraction 
field in the physical optics approach. It equals 


of (2 9%) + en ltr cose—te)__ en ltr con st90) 
with0«p<er—y, 
OF (% Me) pen ere 80 
Oo ws with 5—%<7<8, 
OF (9. ea) -permemer 9 
with e<e- ete, 
Oo; (9. #e) with t+%<9<2 (3.22) 


if one face of the wedge (0<%&<a—*) is illuminated, and 


oF (9% P)+9, (a—p, 2—F)-+ 
+e- the cos ($ —9e_ g— thr cose + Fe) 
with 0<9<e—, 

oF (7. Me) +0, (@—9 2—F) + 

peTrese— 8) With t—- eS Pas, 
Ey | Or (ts me) +97 (3— 9, 2) + 
Fer [pee ee With aac ece, 

OF (% Go) +o (2—9, a— %) + 

ie nese ae te with r<gCla—z—H, 

DT (Ps Pe) +07 (2— 9, 2—%) + 

md cos(¢ — Se) en ite cos (20— 9 — GQ 


with 2a—"—%<pas, (3.22) 


’ 
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if both faces of the wedjte (a--r<e<cer) are illuminated. 


Now let us calculate the field arising duping diffraction by a 
wedge of a plane wave (1.24). The field scattered by the first face 
(9 = 0) is determined by the relationship 


H=—-%, } 
H,=H,=0. | ey) 


One may write the component Hy in the form 


Hes—gHad Ih, (3.24) 


* oy 
or 


H.=— s Hel, with ?<%, 


He= ZH, with @>® 


The quantity I, introduced here is the integral 
1 fvlyi—en de 
=, j Reoty, —e 
2, (3.26) 


along the infinite contour which passes above the pole w = k cos oy- 
This integral, precisely the same as integral I); is transformed to 
an integr2l1 along the contour Do: As a result, we obtain 


OF (8, Peper Th se O Gc E—% 


- = o; (%, %) with ®—%w<9<e, 
Hue U2 (? Pe) with eects, 
le (P Ge)— eth et oS pect, (3.27) 


where 


: iC rin = getter cont 
U2 (9) Pe) = = az j COS ye + cos (6 + W) 
. . 
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In the case when both faces of the wedge are illuminated, the field 
scattered by the second face also ts determined by Equations (3.27) 


and (3.28) in which one need only replace $y by a - o9, and ¢ by 


a- 9. 


Then adding the field radiated by the nonuniform part of the 
current with the incident wave (1.24), we find the diffraction field 


in the physical optics approach. This equals 


5 (2, %) enim (9 ~ toh em iar conte +00 


? | of (7. %) + ev ite cos (¢-%)) 
7.= With #—%<9<e 


Hes 
oF (9. Go) ent cee — 9) 
; with =P EM, 
0, (9. Fo) . with s+9<9<4 


if one face of the wedge 1s illuminated, and 

07 (7, %)-+0, (a—~9, a — %)-+ 

4+ e far cos (5 ~ 184 gottr ccs (9 + 0 

with Vcgqt—h, 

OF (% %)-+ oF (a— 9, a~9,)-+ 

er mtd with *=—-%<Fqa—z, 
He 122 (% Pe) boy (a—9, a—g) + 
ferme vith aac, 
0; (9. %) bos (a—9, a—g) + 
ferret Withec ec 2a—n—9,, 


oF (?. Fe) +o (a— 9, &— 9%) 
us eit tos (¢ ~ 9) ae ete cos (2a — oe — 9) 


with 2a—r—9,< 9a, 


if both of its faces are illuminated. 
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(3.29) 


(3.30) 


The interrals Vie Vy menerally are not expressed itn terms of 
welieknown functions. However, by using the method of steepest 
descents, it is not difficult to obtaln thelr asymptotle expansion 
when kr >> 1, Far from the directions ¢ = 7 + 4 and @ = 2a - Tf - %> 
the first term of the asymptotic expansion gives us the cylindrical 
wave diverging from the wedre edre. In the case of wedge excitation 


by a plane wave (1.23), these cylindrical waves are determined by 
the equation 


] yg 
E,=—H, = E,,-fe. © a od, 
£,=H,=0, (3.31) 


and with the excitation of the wedge by a plane wave (1.24) they are 
determined by the equation 


(we >) 
=F =-H,,.go.$  ___. 
H,=E,==Hu-@ ) lak ° 


H,=£,=0, (3.32) 
The functions f° and 2° have the form 


___ Bin % 
CORY + OS vy” 


oe. See 
i COS p $ Cosy, (3.33) 


if one face of the wedre W<cy<a—) , and 


fozx- Si Y ie Soe sin (a ~~ yo) 
COS y cus yy cos(a— yp vs (t— ye) , 
= _sing " sin(a ak 
Sia cose + cosy cos(a— ,) + cos(a—yo) ' (3.34) 


if both of its faces (a@-1<q¢0<%) are tlluminated. The index "0" for 
the functions f° and n° means that the cylindrical waves (3.31) and 
(3.32) are radiated by the uniform part of the surface current. (3°). 
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§ 4. The Field Radiated by the Nonuniform 


Part of the Current 


In § 1 and 3 we represented the rigorous and approximate expres- 
sions for the diffraction field by integrals along the same contour 
in the complex variable plane. By subtracting the approximate expres- 
sion from the rigorous expression, we find the field created by the 
nonuniform part of the current. It is determined by integrals of the 


type 


{ ¢. ? Pe 6) elt oe Sag, (4.01) 


Ir 
which, with the replacement of the variable € by s=//2e * sins » are 
transformed to the form 
a 


e™ j 9 (2, % a. se" ds (4.02) 


~—o 


and may be approximately calculated by the method of steepest 
descents. 


For this purpose, let us expand the function q(s) into a Taylor 
series 


9 (2, F Pee S)=Fo + 4S +9-5*+... (4.03) 


Let us note that expansion (4.03) does not have meaning only in the 
particular case 


P=2+9, with 9,=0; =, 


g=22—2x—9, withp,=a—R, (4.04) 


when the observation direction ($6) coincides with the direction of 
propagation of the incident wave glancing along one of the wedge faces. 
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Substituting serfes (4.03) Into Equation (4.02) and then perform- 
ing; a term by term Integration, we find the asymptotic expansion for 
the field radiated by the nonuniform part of the current. We limit 
ourselves to the first term of the asymptotic expansion, omitting 
terms of the order (kr)7"* and higher. As a result, the required field 
from the nonuniform part of the current will equal 


o t(# +4) 
Eu= —H,=Eul' ee 
E,=H,=0 (4.05) 


with wedge excitation by plane wave (1.23), and 


(a+) 
ee |, aoa ies 

i aE 
H,= E,=0 (4.06) 


with wedge excitation by plane wave (1.24). 


By calculating, with the help of Equations (4.05) and (4.06), 
the nonuniform part of the current, it is not difficult to see that 
it is concentrated mainly in the vicinity of the wedge edge. But the 
field created in the region kr >> 1 by this part of the current has 
the form of cylindrical waves, the angular functions of which are 
determined by the relationships ‘2) 


P=/—f* g'=g—e@' 


where in accordance with § 1 and 3 we have 
2 at ee 
a . (ae cat atts] 
ere.) 
(*=5) 
(4.08) 


(1) 
Footnote appears on page 42. 
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and 


jo — ne __ 
~~ €d8 g + COS Go 


al f 
e=- TET (4.09) 


if one face of the wedge is illuminated (that is, when %<%<e—= ), 
and 


eae. sin % sin(a —+,} 
I — C08 p+ 508 gy % COS 6 — 5) + CaS (a—y,) ” 
g=— ain sis sin (a — 0) 
“COSY -++COSG_  COS(a— 9) + cos(a—y,) ” (4.10) 


if both faces of the wedge are illuminated (that is, when a—™<qo<a ). 
Let us recall that the functions f and g describe the cylindrical 
waves radiated by the total current — that is, the sum of the uniform 
and nonuniform parts, and the functions fo and 2? refer to the 
cylindrical waves radiated only by the uniform part of the current 
(3°). 

Let us note certain properties of the functions ci and gl, The 
function ri = rica, >, $4) is continuous, whereas the function gi = ge 
(a, 9, 5) undergoes a finite discontinuity when $, = a - 1. The 
reason for this discontinuity is that the uniform part of the current 
differs from zero on the face along which plane wave (1.24) is propa~ 
gated (with $) *a- 7). In the case of radar, when the direction to 
the observation point coincides with the direction to the source 
(o = 95), both functions et and gi are continuous. There is no dis- 
continuity of the function gt with $ = $5 = @-=- 1, because the current 


1 


element does not radiate in the longitudinal direction. 


On the boundary of the plane waves (that is, when @="x+qo and 
q@=2a—n~—Go ) the functions f, f° and g, &? become infinite, whereas 
the functions rt and gt remain finite. In accordance with Equations 
(4.07) = (4.10), they take the following values 
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aa t a 
=——_—_ — +7 n= Zete=. 


4 -— — 
g cos A cos 


(4.11) 


if @2tT - %o» and $9 <a-T, 


ta } 
=—sin— - | 
P =—— +5 ete nt pets $— 
cos = — cos ——— ; 
sin («— 9,) 


~~ cos(a— 9) + cos(e—q) ° 


a= +e, — pets E+ 


cos — — cos 2 


sin(a — 9) 
Sea emeo © (4.12) 


ifo=t- %% andaewm < 9 < 7, and 


P (4.13) 


ifgog=znt %q» and o> <a-mn. The value ¢ = 7 + %% with a-m < $5 
< mt corresponds to the angle inside the wedge, and therefore is not 
of interest. In the direction of the mirror-reflected ray ¢ = 2a - 7 
- %o» the functions ct and gi are determined (with a = 7 < ) <7 ) 
by the following, equation: 


+tetg(a—e) +n tes. (4.14) 
. s 

5 at oe 
ETS sag Oe ag 


ir 
cos n — 0s n 


7 t z 
+7 clg(2—9)—getg= 
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The functions rt and gt have a finite value everywhere, except 
for the particular values ¢ and % enumerated in Equation (4.04). The 
graphs of the functions ri and gi (Figures 8 - 13) drawn in polar 
coordinates give a visual representation of the effect of the nonuni- 
form part of the current which is concentrated near the wedge edge. 
In particular, they show that this effect may be usbstantial for the 
fringing field not only in the shadow region (1+¢o<9Sa,), but also in 
the region of light (O0Se<a+q) . The continuous lines in the figures 
correspond to the functions ft (f+ <0). The dashed and dash-dot lines 
correspond to the functions gi — the dash lines refer to the case 
gi < 0, and the dash-dot lines refer to the case gi > 0. 


apg 


‘a ape 


6-0? 2-30" 


aust 
Figure 8. The diagram of the =4 


field from the nonuniform part 
of the current excited by a 
plane wave_on a half~plane. The 
function fl (or g+) corresponds 
to the case when the electric 
(or magnetic) vector of the 
incident wave is parallel to the 
wedge edge. 


Figure 9. The same as Figure 8 
for the case $ = %o- 


Let us turn our attention to the next important aspect. As is 
seen from § 1 and 3, the nonuniform part of the current on the wedge 
is described by a contour integral which is generally not expressed 
in terms of well-known functions. But in order to calculate the field 
scattered by some convex, ideally conducting surface with discontinuli- 
ties. (edges), the indicated expression still must be integrated over 
the given surface. Obviously, such a path is able to lead only to 
very cumbersome equations. Therefore, henceforth, when calculating 
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= 
a: 


Figure 12. The same as Figure 10 
for the case $ 2 9° 


Figure 13. The functions rl and 
gt for a wedge (¢ # o> a = 210°). 


the field scattered by composite bodies, we will not integrate the 
explicit expressions for the nonuniform part of the current, but we 
will endeavour to express these integrals directly in terms of the 
functions rt and gi which have been found. 


§ 5. The Oblique Incidence of a Plane Wave on a Wedge 


Above, the diffraction was studied of a plane wave incident on 
a wedge perpendicular to its edge. Now let us investigate the case 
when the plane wave 


Ean E e208 4 6508 + cond) (5.01) 


falls on the wedge at an oblique angle 1 (0<1<3) to the wedge edge 
(Figure 14). 


From the geometry of the problem, it follows that the diffraction 


field must have that same dependence on the z coordinate as the field 
of the incident wave, that is 
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E=E(«x, yyelttert, ° 
ee (ore) 


Using Maxwell's equations 


rot H=2 —ikE, rotE=ikH, §=(5.03) 


Pigure 14 Diffraction by a one is able to obtain the follow- 


wedge with oblique incidence of ing expressions for the radial and 
a plane wave. y is the angle ‘ 
between the normal to the incident azimuthal components of the field: 


wave front and the z axis. 


! 1 OM, o£ 
E-=— waar (F op beet). 
1 toe dHe\ 
y= ra (> Gy — SST ~G 
| on cosy OF 
§,= may (ar PH): 


wa! _ (28s 4 097 Me 
H + Sm): 


¢ asin®y\ or "FP (5.04) 
The functions E, and H, in turn satisfy the wave equations 
SE, +e E,=0,4H, +8} H,=0, (5.05) 
where 
beste and k,= esiny. (5.06) 


In § 4 we found the fields (4.05) and (4.06) which satisfy the 
equations 


AE,+4E,==0, 4H,-+4'H,=0 (5.07) 


and which are created by the nonuniform part of the current excited 
on the wedge by the plane wave 
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E = E,e7!* (seen te + ine (5.08) 


Representing Expression (5.01) in the form 


E am Eg! 7 ites conte + 7 0in ve (5.09) 
and comparing Equations (5.05) and (5.07), we easily find the field 
created by the nonuniform part of the current with the irradiation of 
the wedge by plane wave (5.01). For this purpose, it is sufficient to 
replace in Equations (4.05) and (4.06) k by k,, and E,, and Hoe by 
Eye" and Hye" . As a result, we obtain 


e 
i ayee 


E,= —H, = | a (9. *) TS eltt cory, | 


‘(4 +7) 
= _ ' eg ttt cosy 
H,=EF,= Hyg (?, Po) ¥ icky e ° | 


(5.10) 


The angle > introduced here is determined by the condition 


eittremes COS %) __ qth (x con ay + y sing.) P (5.11) 
hence 
__. cos B 
1S Pe cee (5.12) 


The remaining components of the field created by the nonuniform 
part of the current witn the oblique incidence of a plane wave are 
found from relationships (5.04), and when kr >> 1 they equal 


E,=—ctgy£, H,=—ctgyH, 
Ay Hee, (5.13) 


The equiphase surfaces for these waves have the form 


rsiny +2 cosy =const (5.14) 
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and are conical surfaces, the reneratrices of which form the angle 

n/2 + y with the positive direction of the z axis. Thus, with oblique 
irradiation of the wedge by a plane wave, the field created by the 
nonuniform part of the current is a set of conical waves diverging 
from the wedge edge. The normals to the phase surfaces of these waves 
form an angle y with the positive direction cf the z axis and are 
shown in Figure 15. These waves may be represented in a more graphic 
form if one introduces the components (see Figure 15): 


E, = £, cosy — E,siny, 
Then the final expressions for the fringing field in the far zone will 
have the form 


p= Ea te (5.16) 


Now we are able to proceed to the application of the results 
which have been obtained for the solution of specific diffraction 
problems. The simplest of them is the problem of diffraction by an 
infinitely long strip which has a rigorous solution [23] in the form 
of Mathieu function series. However, in the quasi-optical region 
when the width of the strip is large in comparison with the wavelength, 
these series have a poor convergence and are not suitable for numeri- 
cal calculations. Therefore, the requirement arises for approximation 
equations which are useful in the quasi-optical region. The derivation 
of such equations for a field scattered by a strinv will be given in 
the following section. 


§ 6. Diffraction by a Strip 
Let us investigate diffraction by an infinitely thin, ideally 


conducting strip which has a width of 2a and an unjimited length. 
The orientation of the strip in space is shown in Figure 16. 
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Figure 16. Diffraction of a 
plane wave by an infinitely long 
Figure 15. The cone of diffracted strip. The section of the axis 
rays. y (-a < y < a) shows the trans- 
verse cross section of the strip 
with the plane z = 0; a is the 
angle of incidence. 


Let a plane, electromagnetic wave strike the strip perpendicular 
to the edges. Let the direction of propagation of this wave form an 
angle «(la]<) with the plane y = 0. The field of this wave is 


represented in the form 


E= E,e*“ coset gains) H=H,e"* (scone + ysine) (6.01) 


The uniform part of the current excited by the plane wave on the 
strip has the components 


p=o, 


i —_ go Hos elt” ela * 


thy singe 
s 


I= ge Foscosae ; (6.02) 


Substituting these values into the equation for the vector potential 
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a 2 =-_ 
= j Piddy (eee 
Acer P49) yaya” 

J Jy: +u—v' +6 (6.03) 


and taking into account relationships (3.05) and (1.18), we obtain 
the following expressions in the region r »>> ka’; 


A, = 0, 
s(ee3) 
A, = 74, sin [ta (s'n a — sin 5)} © 
~ & slas—sing Y2:ar ° 


i (aco— 
sta [ta (sine — sing] es ) 


2 
A= ry Ey, cosa sina— sing Y 2-80 


(6.04) 


The components of the fringing field in the cylindrical coordinate 
system equal 


E,=—H,=ikA, Hi=E,=ikA,, (6.05) 
where 
A, = A, cos 9 — A, sing. (6.06) 


Substituting Expressions (6.04) here, let us determine the field 
radiated by the uniform part of the current 


, (aro) 
as sin (4a (sina — sin 1 @ 
Pm =—H, 22, 68 ae ¥ EE’ 


angen ; (t+) 
= = sin a (sina — ala }} e 
E, = Hi= 2H, cose Sa sina—eing = “Yasukr (6.07) 


This field may be represented in the form of cylindrical waves diverg- 
ing from the strip edges 
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F, vam ly: Eu: Exe (f(eors: - sine) 3 
lors +) 
“YW va 

Eva M, 2 Mea (gr(lpemrame— ry 


seed) 


View 5 (p08) 


+ fe(2ye da atne-- a fs : 


+¢° (je -(taqeine-- ans 


Here the first term: correspond to the waves) from edge lo (yo 2 oa), and 
the seceomd terms correspond toe the waves from edge 2 (y # <a). The 
funetions rooand gooare determined In the right: half-plane (i71<3) 

by the equat Pons 


PP (Wer -- f9(2) Pres 


sf cone 
ef (1) == ~- g*(2)- *pine—s.ag’ (@.aa) 


Now fet ous f nd the ffLeldo radiated by the nonuniform part of the 
current. Assuming the strip ta suftfefent ly whde (kas? LD), ome ts 
able to approxfmately consider that the current near fie upper edee 
Is the same ascon the fdenily eonduct tag half-plane --esyea , and 
near the lower edpe ft fe the same as oon the dni tf-plane --a es vec, 
Thevefore, En accordance with § Fy the tleld from Che nonuniform port 


of the eurrent flowlog on the strip may be represemted fm the form or 


the sum of the edee eytindrfeadl waves. 


E,: Wye _ [prayeeers- sing) - 


(rs) 
ta ~ (ta (sine. - sing) 
saye sl rea 
Eye: Hy-= Hee le'(\) eltiisine- abla 
(ees -) 


7 tha (ai e--sing) af 
-fa'(-e ars 


(os ray 


where the Punet Pens r! fund yr! are determined dno the rfiht tend tp dane 


. & 
(Ivis 2) by the eqpiunrt Pons 


SoU | hematoma rey rea 


PY=FY— 1). PRA=1A—-FE). 


‘ 6 e 
a (Y=g(I)—e(1), g'(2)=£(2)—6°(2) (6.41) 
in connection with which 
cus : — sin = 
)== ~ “naam my 
cus “+? +sIn —" 
M2)=5 tas eng 
coe ote -- sla —t 
8()= ~~ asiny 
— cos" $9 + sin tS? 
6Q)=--— igs aey 
(6.12) 


The functions r° and ge are described by the relationships (6.09). 


\ 
As a result, the fringing field (the sum of the fields radiated 
by the uniform and nonuniform parts of the current) will equal 


(=—i.= Ex-[f (leone | 


+ f (2) ew ite (sin ee sin | is) 
binke ” 


Ey = Ha Hes lg (Iyer + 


. Qy en Mtasine—stagy et 
+4 (2) ya (6.13) 


Consequently, the resulting field is expressed only In terme of 
the functions f and -¢ which determine the cylindrical wave in the 
rigorous solution (see § 2). The field Lis the superposition of two 


such waves which diverge from the edges 1 (y = a) and 2 (y = en). 


Substituting into Equations (6.12) the explicit Expresstons (6.1.°) 
for the functions f and pf, we obtain 
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E,=—H, =, cos [ta (sin a — sin ¢) 
eet 


a( (14>) 
4 i sin[ta (sta 2 — sin ¢)} 


sia +52 " p 2.20 
—_ — a— 
E,= Hi=He ‘a {ka wit y) + 


afeed) 
+i isis (4a (sin 2 — sin a—sin»)) 


¥ 2a 
a = | (6.14) 


These equations are valid when r >> ka? and lel<z. Moreover, it is 
assumed that ka >> 1, since only under this condition is one able to 
consider the nonuniform part of the current in the vicinity of the 
strip's edge to be approximately the same as on the corresponding 
half-plane. In the case of normal incidence of a plane wave (a = 0), 
Equations (6.14) change into expressions corresponding to the first 
approximation of Schwarzschild [15]. 


From relationships (6.03) and (6.05), it follows that the elec- 
tric field is an even function, and the magnetic field an odd function, 
of the x coordinate measured perpendicular to the plane x = 0 (in 
which the current flows) 


E, (x)= E,(— x), H, (x)= — H,(—x). (6.15) 


Therefore, on the basis of Equations (6.14) and (6.15) one is 
able to write the expressions for the fringing field in the region 
2 
x <0 (where 7Sl?/<#) 


i. | 


sis 3 


j sinha (stn a —~ sin 7 ‘(+a) 


cos Y tbr 
E,=Ha= Ha: eet + 
ala 2 


(6.16) 
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4 Malietlos = sinh (ers) 


atypia 


(6.16) 


Here one must select the upper sign in front of the braces when 
¢ > 0, and one must select the lower sign when @ < Q. 


The resulting Equations (6.14) and (6.16), 1n contrast to 
Equations (6.07), satisfy the reciprocity principle. It is not diffi- 
cult to establish this by verifying that Equation (6.14) is not 
changed with the simultaneous replacement of a by ¢ and of ¢ by a, 
and Equation (6.16) is not changed with the replacement of a by t + 6 
and of ¢ by a- 7 (if —<9<~-7 ) and with the replacement of a by 


wn - 9 and of ¢ by n- a-(if 3<9<e). 


However, the indicated equations lead to a discontinuity of the 
magnetic vector tangential component H, on the plane x = 0. This is 
connected with the fact that, by considering the nonuniform part of 
the current in the vicinity of the strip's edge to be the same as on 
the corresponding half-plane, we actually assume the presence of 
currents on the entire plane containing the strip. In order to refine 
the resulting expressions, it is necessary to solve the problem of 
secondary diffraction — that is, diffraction of the wave travelling 
from one edge of the strip to its other edge. In other words, it is 
necessary to take into acccunt the diffraction interaction of the 
strip's edges. As we see, it is also necessary to take into account 
the secondary diffraction in the case a=% when the H, component 
of the fringing field must equal zero. 


In Chapter V, we will return to the problem of diffraction by a 
strip, and together with the investigation of the secondary diffrac- 
tion, we will present the results of the numerical calculation based 
on Equations (6.07), (6.14) and (6.16). 
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FOOTNOTE 


Footnote (1) on page 27 The designations used here differ 
slightly from those used in the 
papers [7 - 11]. The functions f 
and f+ there were designated by 
fy and f, respectively. 
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CHAPTER IT 


DIFFRACTION BY A DISK 


The problem of diffraction by a disk has a rigorous solution 
[24-26]; however, it is not suitable for numerical calculations in 
the quasi-optical region when the dimensions of the disk are large 
in comparison with the wavelength. The physical optics approach 
used in such cases sometimes gives erroneous results. In particular, 
the fringing field calculated in this approach does not satisfy the 
reciprocity principle. 


In this Chapter a refinement of the physical optics approach is 
carried out. First the diffraction of a plane electromagnetic wave 
by a disk with normal incidence (§ 7-9) is investigated, and then 
(§ 10-12) diffraction by a disk with oblique incidence of a plane 
electromagnetic wave is investigated. 


Normal Irradiation 


§ /. The Physical Optics Approach 


Let an ideally conducting, infinitely thin disk of radius a 
(Figure 17) be irradiated by plane wave 
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£,=—H,=—H,,e™", 
£,=H,=0. (7.01) 


The uniform part of the current 
excited on the disk by wave (7.01) 
is determined by Equation (3.01) 
and has the components 


h=—2Ms. 2=2=0. (7.02) 


Figure 17. Diffraction by a Let us find the field created by 
disk of a plane wave propagated 
along the z axis. this current. 


Since the diffraction fieldin the far zone (R >> ka“) is of 
interest to us, the vector potential 


e au 
A(x, ys z= [dp (1, 9) — 4 
9 y e J J 7 r (7.03) 
may be simplified by using the relationship 
r= K+ 9° — 2pR cosO = R —pcosO, (7.04) 
where 2 is the angle between p and R, and 
cos OQ =sin § cos () — 9). (7.05) 
As a result, we obtain the simpler equation 
- @ 22 
1 ft? be ~ikp cos ® 
Mey =a [och | ee, ern a (7.06) 
0 
Continuing by using the equations 
H=rot A, rotH =—ikE, (7.07) 


it is easy to show that in the spherical coordinate system the 
fringing field components with R >> ka? equal 
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E,=H,=ikA,, 
£,=H,=0, 


(7.08) 


where 


A, = A, cos@— A, sin 9, 
\ } (7.09) 


A, = (A, cos ¢ +- A, sin 9) cos 8—A, sin 9, 


Substituting here the values 


. “4 . 
4,2 -- Hea pe A(hasindy (7.10) 


A, = A,=9, 


which result from Equations (7.02) and (7.06), let us find the field 
radiated by the uniform part of the current in the form 


E,=H4,= ~ iaHes 4 = ay cos 0/, (ka sin no . 


E,=—H, seco: 2 thesia) | eit? 


The function J, (ka sin #) is a first order Bessel function. By 
using its asymptotic expression 


Atha sin) =|/ —2 cos (ka sin — *), (7.12) 


nko sind 
which 1s applicable when ka sin >> 1, one is able to rewrite Fqua- 
tions (7.11) in the following form: 


sin 
E,==H,=— ith: V aaa a i; 2 cos 0 NOLX 


feito, ome 


Pes et | (7.13) 
SH Ty SS oT oe | 


¢ Fak sind sind * 
—1 (+ na errr) s 
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The resulting equations show that in. the region R >a’, kasin§ >! 
the fringing field may be investigated as the sum of the spherical 
waves from two "luminous" points on the rim of the disk, the polar 
angles of which respectively equal y = ¢ and » = 7 + 9. It is not 
difficult to see that thesewaves satisfy the Fermi principle. 
Actually, of all the points on the disk's surface, the point po # a, 
¥ = 9 is the closest to the observation point (R.% 9) , and the point 
90 @a, » = 1 + ¢ is the furthest from it. 


However, Equations (7.13) describe the radiation not only from 
the two “luminous” points, but they determine the field radiated by 
the entire "luminous" region which is adjacent to the line connecting 
the points p =a, » = ¢ andp =a, punt o. 


Let us show that the luminous region actually makes the main 
contribution to the fringing field. For this purpose, let us cal- 
culate the field radiated by the currents which flow inside the 
sector encompassing the line » = » (Figure 18). Let us take the 
angular dimensions of the sector in such a way that its arc, which 
equals 2a%o,» would occupy the first Fresnel zone. When this is done, 
the angle 9% will satisfy the equation 


a(! cos @) sind =, (7.14) 


In the case being investiga- 
ted by us, when the condition ka 
sin >>! is fulfilled, we have 
from Equation (7.14) 


ey er rere es (7.15) 
hence 
Figure 18. Calculation of the 
field radiated by the "luminous" = (7.16) 
region of the disk. ia po . 
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The vector potential of the currents flowing In the Ind!teated 
sevctcr is determined by the equatton 


es, 

pee pee H, meni d —ihpain 8 cos 2, di 

ste AR jaf av (7.17) 
A, = A,=O 


Taking into account the condition ka sin ®>!, one may show that the 
field created by the currents of this sector will equal 


a 
Eo, ~He hens 89 X 


tap elt® than d He 1 Yy 
slaze ao Hs o( ) 


sae R V kasint , 
E,w—Ny~Ha sina X (7.18) 
cose ef#% = ~itastn O44 1 
Xaab a ° +0( sean): 


The amplitude of the expressions which have been found is approximate- 
ly 2 times larger than the amplitude of the first terms in Fquation 
(7.13). Moreover, expressions (7.18) and the corresponding terms 

in Equation (7.13) differ slightly in their phases: the first have 


the factor As oe and the latter — the factor —e* , The result obtained 
is similar to the well-known thesis in optics that the effect of a 
wave 1s equal tc the effect of half of the first Fresnel zone (see, 
for example [27], p. 132). 


In the vicinity of the directions 0=0 and O=2, when the 
azimuthal components lose their meaning, for the purpscse of studvine 
the fringing field it is more convenient to use the Cartesian 
components . 

E, = (EF, cos$-+ £, sind) cos 9—E, sing, 
E, =(E,cos9 + E, sin) sing 4-E, cose. (7.19) 


Turning to Equations (7.11), we find that when §=—OQand d=e 


E,=0, E,=—iH San Oey (7.20) 
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Conseyuently, in the physical optics approach the field scattered 
in the directions 0=0Qand %-=" preserves the polarizaticn of the 
incident wave. 


§ 8. The Field From the Uniform 
Part of the Current 


Let us proceed to calculation of the field created by the non- 
uniform part of the current with normal irradiation of the disk. 
Since the latter is concentrated mainly in the vicinity of the disk's 
edge (9 = a), the vector potential corresponding to it will equal, 
in accordance with Equation (7.06), 


ao Rn : 
Anti. 5 dp J F(p, geTem bens ag, (8.01) 


The inner integral is calculated with ka sin® >I based on the 
stationary phase method (see, for example [21], p. 256), and Equation 
(8.01) is transformed to the form 


PD oe eee 
A(x. y a= t is sia’ ze ‘ x 


' ; e (8.02) 
xlf P(e, Ys) eeeeraytTy (p. %) el ea] ; 


which allows one to interpret the fringing field as the field from 
a luminous line on the disk. This line is a diameter, the polar 
angle y of the points on which equals 


1s =? and $,==2-+9. (8.03) 


Assuming the diameter of the disk is sufficiently large in com- 
parison with the wavelength (ka >> 1), one may approximately assume 
that the nonuniform part of the current near the disk's edge will be 
the same as on the corresponding half-plane (Figure 19). On the basis 
of § 4, the field from the nonuniform part of the current flowing 
on the half-plane ~°Sy¥,“a@ may be representedinthe form 
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e ’ 
i alte 


—ika ind 


E, =IkA, (N= E,, (ly aaa : 


H, (1) —IkA, (1)cos 8 = 


- aney —iha da 6 
(1) Vink e ‘ 


(8.04) 


= H,,,@° 


and similarly the field from the current flowing on the half-plane 
—@Sy,S, may be represented in the form 


dd 


——_—_gita sin 9, 


EQ —=URA, (2) By 1S 


: : ana?) (8.05) 
Hy, Aya — kA, Osten Ke g'(2) ree ot giaa sine , 
Here 
i (nat. 
A(t) ty 2 ( 7) f sige, 
(8.06) 


‘ (seez 
A(2) = 7 '( ‘)t f Pipette dn, 


and the functions rt and gi are determined for the right half-space 


(o<t<+) by the equations 


cosy + sla 
PM=K)+o%, as 
5 (8.07) 
. cos —~.—sto 
P\)=(%)—sy. Pcl aia 3h 
cosy — ala g- 
a' (1) = 9 (1) + Sey »@()=—-—_— . | 
shia: a (8.08) 
a'(2)=g(2)— Sh, (2) —__;~—. | 


From relationships (8.04) = (&.06), it follows that 
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f ii, (4) em aa sin® dy= | 


tle —tha sin ® 

=— 75 Ef (Ie ’ 

- (8, 09) 
fA perry 


—-@ 
=— 5 E.,I' (2) e!*° and 


and 


@ 

i, (n) eit y, sind dy= 
Figure 19. Diffractien by a disk. a 
The half-plane L lies in the 


=, 4 3 sika sing 
plane of the disk. The edge of ~Teneond /Yo,,8' (Ie , 


(8.10) 


the half-plane is tangent to the a ide gin ® 
circumference of the disk at the fi,@e i aS 
point y, = a, x, 2 0 (a is the ra ae 
asin 
radius of the disk). “Tarcost /t,,'8' (2)€ , 


In accordance with the assumption of equal currents on the disk and 
on the half-plane, one may consider the following equalities to be 
valid: 


fro. 4) ea trein® 4, _ i J (» e!™ 198 on, 


e e ; 8.11 
JO sel Mdm J Fiperen'dy oo 
=a 


Therefore, the field from the nonuniform part of the current flowing 
on the disk will equal 


4 (ta vin 7) - 


£,=-—H, (2)e 


— 4% oe ’ 
oF V 2:.24 sin 6 [i 


aye (eon F a 


R.* 

(8.12) 
_ yy laMe, (ta un o—7) 
E,=H, = [2 @e = 
3a 
mete rt) eae 
where in view of (7.01) 

E,, = — Hz cos9, H, =— Aes sing. (8.13) 
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For the direction§=0, we have according to Equation (8.01) 


(8.14) 


2s a 
Ant. fa fare. )a>, 


but in accordance with equalities (8.09) = (8.11) 


i 9) dp= if jag Hen-C08 4, 
(8.15) 


fe I (. 9) do= 25, Hye sing. | 


Consequently, 


P e 
=e. Jeosvae Sf. yde+ 


ye a i (p. 9) a|- 
(8.16) 
A= e|- -fensafe (p. ue 


+ f sin $y fi Pe. nal =0, 


that is, in the direction of the main fringe (@=0) the field from the 
nonuniform part of the current equals zero. 


By using the Bessel functions Jy and Jo for the field from the 
nonuniform part of the current, one may write the equations 


Ey =— Hy = (2) —F* (4 (asia 8) $- 
+i Q)+P (LY, (eosin Oy} , 
E=H,="7 = {[a'(2)— a" ee 
"+ Ha! (2)-+ g (I)], (ta sin SS , 


(8.17) 


which with ka sin ®>1I change to Expressions (8.12), which were already 
found. In the direction0=0, these equations give a field which 
equals, in accordance with (8.16), zero, and with intermediate value:.: 
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Subst fitat ine here the explfielt expresstons Cor the funet fone 


foand yy we aredve at the final expresstons for the Cringtny pled 


E, me H, = -— iallye yy (a stne) 
oe 
Ja: J, (tas te 8) pe sens, 
cos + | 
E, =H, = — ‘then [nGese) 4 (9.03) 

sla-5- 

+i J,(ka somes a) Fe ~ sin ?- 
cos > 


These equations are valld in the ripht half-space (0<0< +) - Inthe 
left halfespace (4 <?<s) » the rringing rleld ts eastly found by 
assuming that Its eleetrle field fs an even function, and Its marnet te 
fleld an odd function of the « coordinate: 

E, (z2)= E(-— 2). 

Hy (2) =—H,(—2). (Q,04) 


a 
Consequently, In the region z < 0 (that fs, when ~7 S88) 


laHes [ 1 (ta stad) 

E,=—H,=- aac I, (a ~~" 
cot 5° 

pres |e cos, 


sla 5 t 
. away 
, Ey=H, = lattes | Jutta sia ne 
a 
ig pea = sine: 
sia 


Assumtne that tn Equations (9.03) and (9.05)@==0 and dea, renpest ively, 


we obtain 


jka® eft? 
| a = 7 ‘Hes: “a . E,=0, chow 


whieh Ls equivalent to the physteal opt tes appreach [see Fquat fon 
(7.20) ]. 


MDH a se Oe Ty A 


Expressions (9.03) and (9.05) agree with the result obtained by 
Braunbek [29] for the scalar fringing field in the far zone. It is 
also interesting to compare these expressions with the precise numeri- 
cal results obtained by Belkina [34] by the separation of variables 
method in the spheroidal coordinate system. It turns out that even 
with ka = 5 a satisfactory agreement is observed between our approxi- 
mation method and the rigorous theory. In Figures 20 and 21, graphs 
of the functions V"(0) and V' (8), are presented which allow one to 
calculate the fringing field on the basis of the equations 

E,=—H, = ou (8) o= cos 9, ) 
Ey=h, = He, vO sing. | 
The continuous curve corresponds to the rigorous theory [34]. The 
dash-dot curve corresponds to the field from the uniform part of the 
current, and the dashed curve corresponds to the field calculated 
according to Equation (9.03) and (9.05). 


(9.07) 


Oblique Irradiation 


§ 10. The Physical Optics Approach 
Let us investigate the general case when the plane wave 


Em £E, elt (wistn'y +2 cos (10.01) 


falls on the disk at an arbitrary angle to its nxis. Let us take 
the spherical coordinate system in such a way that the normal to the 
incident wave front, n, vould lie in the half-plane p=-> and form 
an angle y (o<1<-) with the z axis (Figure 22). Adhering te the 
investigation procedure used in the previous sections, let us first 
calculate the fringing field in the physical optics approach. 


The uniform part of the current excited on the disk by wave 
(10.01) is determined by Equation (3.01) and has the components 


— ¢ Ly iky stn : 
R=— Hy, f=—sH,,e™ t 1 f=0. (10.0.) 
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Figure 20. The function ¥@@ for a disk 
with a normal incidence of the wave. The 


varicus curves correspond to different 
approximations. 


The field radiated by it is found, as was done in § 7, by interrating 


(with the condition R >> ka’), In the case of E=polarizaticn of the 
incident wave (E,lyoz), this field equals 


ne 
Ee = H =. iaE,,-cos 1 cos 8 cos 9 Lee oy) - 

2 7 
J,(ka pO Epy alt (10.03) 


E,= — Ay=—iak,,-cosysin¢- ii 
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Figure 21. The function V@ for a disk 
with normal incidence of a plane wave. 


The various curves correspond to different 


approximations. 


and in the case of H-polarization (HN, Lyoz) 


ikR 
E,=H =~ iH, cos ® sing MG2V Et a , | 


Yeni RE 
oar eee Nea VIF EaH otk 
E, == Hy = ~iall,, cos p+ Vii | 


=sin 4 sing —siny, 


4==sin 9) cos 9, | 
VY Pep =0. 
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C1004 5 


The quantities A’ and yw in Equations (10.03) and (10.04) are deter- 
mined in the following way: 


(19.95) 


Assuming 92-5 and b=2—1(3 <0<*), in the resulting expressions, 
let us find the field scattered by the disk in the direct!on toward 
its source. With E-polarization of the incident wave, it equals 


. ine 
E mH, = ME 8p aoa sin ye 


E,=H, =0 (10.06) 
and with H-polarization 
sae 
Ey, = She ct’ jitkasiny’s, | Ean 
£,=H,=0. 


Using the asymptotic expressions for the Bessel functions, one 
is able to show that when R >> ka” and ka Ep > 1 the fringing field 
is radiated from a luminous region on the disk. In the case when 
oY P+ pt, the luminous region is increased and in the limit (when 
A = » = 0) the entire surface of the disk starts to "shine". 


§ 11. The Field Radiated by the 
Nonuniform Part of the Current 


Let us calculate the field in the nonuniform part of the current. 
‘J (p. Jee, gpeite nts, (11.01) 


Its corresponding vector potential 


A= fats, gear ret Ne (11.02) 
by means of the stationary phase method is transformed with 
ka/P+e>1 to the form 
Fy AE Pe 
=F wrireae ‘[Jiete ia 
; ete (11.93) 


a 
- ifse. Zaye Pap | a 
1d 


Here 
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= 3, Y=E-s (11.04) 


are the stationary phase points 
and the quantity 6 is determined 
by the equalities 


nee rie » cos8=—~—---—, (11,05) 


Figure 22. The oblique inci- From Equation (11.03) it 


dence of a plane wave on a disk. 2 
nis the normal to the incident llows that with R >> ka” and 
wave front. ka /#EpS 1 the main contribution 


to the fringing field is given 
by the luminous region adjacent to the line y = Vis v = o> Thus, 
the stationary phase points vy» V5 physically correspond to the 
luminous line on the disk surface. 


In order to calculate the vector potential (11.03), it is 
necessary for us to first express the nonuniform part of the current 
on the half-plane in terms of its field in the far zone. For this 
purpose, let us introduce the auxiliary coordinate systems Xi» Vy 
and X55 Yo (see Figure 23), and let us take the following designations: 


Qi» B, (a5; Bo) are the angles between the normal to the incident 
wave front and the coordinate axes Xp» Vy (x5, Yo)3 


oled = -99) is the angle between the z axis and the projection 
of the indicated normal on the plane x, = QO; 


$3 (6) = -$5) is the angle between the z axis and the direction 
from the coordinate origin to the point P(yy> z) which lies in the 


plane x, = 0 and is the projection of the observation point P(x, v, 2); 


rm is the distance from the origin to the point PCY) > om 


The quantities introduced here are determined by the equations: 
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CUS 4, === SiN YCOS},, cos 9, ==sin 7sin4,, 
a 3, =t—$,, 
sing, = ¢! dy “us W cuny 


sina,’ PAT sine, 
tin ters( — 9) 
sing = Viton isinin—@ (11.06) 


= coed 
Vira dentG—n 
iF, ar 


r,=3RV 1 —sin*6 sin®(), —9) - 


COS 7, == 


Furthermore, let us write the expressions for the field from the 
nonuniform part of the current excited by wave (10.01) on an ideally 
conducting half-plane —co<y,<a. In accordance with § 5, they have 
the form 


aay 
E, seas efMt Ce & ent (r» ? ) se ene 9°~ein ¢,) : (1 1.07 ) 


(ues 3) iha(sta g®—sta +) 


H,, en et CBR or (Fis )° Vas e ° 
i] 


where 


k= ksina,, (11.08) 
k, (sin 9, —sing) =k Yi? Ep, 


sin non we nt ul e 
2 wm 


e cocoa «Ee! | 

P(t 9, sin y} — sin y, sin ¢f — sing, ® r 
) 

eas " ree eas | 

] 

7 

, 


2 (A> j= ae a ee ee eres 
sin 9) — sing, sing, —siag, ” 


(11.09) 


p de 
(-#<n <4) 
On the other hand, this field may be expressed in terms of 
the vector potential 


‘y ib s,cos § int cove, COE FPR 
Aap] inet hy | eM ame ® = (11-10) 


—@ —2— 
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By means of equation 


which follows from Equation (3.10), if one substitutes 2 = t, w= ¢, 
ad = -ip, k = -iD, in it we find that 


ts iy , ‘2A 4ly—wte' 
Aa Ems (Hy Hey FFG eM May (11.12) 


Taking the fact into account that 
the nonuniform part of the current 
is concentrated mainly in the 
vicinity of the half-plane edge 
and using the asymptotic repre- 
sentation of tne Hankel function, 
we obtain 


1 tn" iA(x, 608 abr, a1 eHet ° ihyd, 
Ami jee fsme 4% (17.13) 
-@ 
Figure 23. Diffraction by a disk where , = cos 4, -— sin 3, Sin p, = 
with oblique incidence of a plane 
wave. The half-plane L lies in =sin {sin ¥, — sin 9 cos(p, — 7) 
the plane of the disk. Its edge hs etait ee ALE 
is tangent to the circumference V eee. (11.14) 
1 — syn? d sin* (}, — 9) 


of the disk at the point x, = 0, 
y, =a (a is the radius of the 
disk). In the case when v= 6 [see 
Equations (11.04) and (11.05)) 
the function 1 takes the value 
O,= — f/f (11.15) 


Starting from expression (11.13), it is not difficult to shew 
that the fringing field in the far zone is described by the following 
equations: 
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E, ==— iksina, cosa, sing,A, +(e sin®s,A,, (11.16) 
H, =s— ik sina, cos9,A,, : 
where 
A,, =, (¥,). A, =O, (bs 
nal oe (12.17) 
=ty =e (&(e, cos Sey 
é 
1,03 f 4,,(0 oF Ay 
“— (11.18) 


I a (9,) = j J, (y) eam 


Then by equating expressions (11.07) and (11.16), we find the desired 
connection between the nonuniform part of the current on the half- 
plane ~©<¥,“@ and its field in the far zone 


H# vig 
1, @=— Ea Tine con py (Pe 1 em a : 


pera (En. f* us ~ 


(11.19) 
—cosa, tg eH, g'(9,. 9° ener ry, 


One may show in a completely similar way that the nonuniform 
part of the current excited by wave (10.01) on the half-plane 
—co<y,<a creates, in the far zone, the fringing field 


w= tk bes 2, cos2,sin?,A, -ék sin’a,A,, (11.20) 
. == —iksin 2, cos p, A,, 
where 
A,,= 8,1, (ta) A,, = 941, (42) 
(11.21) 


1 = ikte, «Ina, = 0, cop ai > 
= 7k inn, e ’ 


ene Be a 
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a . ws 
FA 9) f4,(et Fae, | 
-@ 


. ‘ieee (11.22) 
Iv) = f 4,(we' dae sae 
-= 


On the other hand, in accordance with § 5 this field equals 


eee +) thya(sin "sin ¢y) 
E,=e™ as * E.,J' (Ps, ? ) aS 4 5) 


(11.23) 


than ¢4—sln Gy) 


Hee ibevconn H,,8'(e 7, ) Be ae 
he 
Here 
k, (sin 2, ~ sin 9) )= — Yap yp", (11.24) 


and the functions /*(7, 9) and g'(%9%) are determined by the equations: 


ata con tH ? 
2 2 cus 9) 
PQ %2)=— a ace 
#1 — 810 & sin ¢; — 91 9, 
0 0 
n-th at 9 
“ne PPM het ish Lie Soe) (11.25) 
oe sing —sing, sin 79 — sing, 
‘3a a 
Vrs <}). 
Equating the quantities (11.20) and (11.23), we find 
c etait Te ; 
A (C= = Fas “sine, con 9 Hg! (P., Vo )» 
: ikaV RF wt io a8 
1,0) = Gage ata, — Enel" (Pas PI (11.26) 


C082, te Po, 8! (Par Pe II- 


In this way we established the relationship between the nonuniform 
part of the current on the half-plane and its field in the far zone. 
Now let us return to a calculation of the field from the nonuniform 
part of the current flowing on the disk. 
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Since the disk is assumed to be large in comparison with the 
wavelength, the nonuniform part of the current in the vicinity of its 
edge may be approximately considered the same as on a corresponding 
half-plane. Consequently, the integrals in Equation (11.03) will 
approximately equal the corresponding integrals from the current on 
the half-plane: 


j I, (eo. $) ev Vee wap <y !,, (¥.). 
J+. (, yer P+ Map = 7 (,), 


: es (11.27) 
f4.00 0 YF Map, Mi 


Sint ge FR ay I, (dy). 


As a result, the vector components of (11.03) may be represented in 
the following form: 


I = eltR 1 i 
A= V Wire re “tl, GFE, Cdl | 


' (11.28) 
I Uy ia 
A,=+Y SS ; a aid Hil, Gal | 
Then substituting these values into the equations 
E, = ikA, =k [A,, sin (7, — 9) — A, cos(}, —9)I. (11.29): 
E, =ikA, =ik (A, cos (¥, —9) + A, sin (4,— 9)] cos ® ; 


we find the field from the nonuniform part of the current flowing 
on the disk 


i= 
eS __4 ale sin(t; —9) 
E,= H,= VY ie | i, (a 
4 
i te *) cos 3, tgr.|X 
Xla' (es 9; enateing ig' (5 ¥, ) eel en] (11.30) 
En, a UP a Demi RH — 


side’ if' (Pes ?3 ) ae ; 
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E,= H—-— ae * cosa 


: sh BE EOS 2 ish f 


— cus sy - — —#) 
YV2rka Vit tpt “Ee \ Ho,, E SIN dy COBY + 
+ — sin 1%) — 9) ~ 9) 


sinta, O80 tg 9 Jie" (7. ? 1 penta” uO 


— gt ( ?; yeltal’ Ors + +E, sin(4, 


- (11.31) 
hate 
x<if (?,. ?, ) enter __ if* (Pas ?, ) elisa ives ] 


The resulting expressions are valid when ka/i*ty" >1. 
may be slightly simplified to 


They 
ata eine 
E =H. = * 
OV tnt pe ® “RX 
X{— Has, cos 4 Set Ig! (Pas F ° ye ieaV ite 
— igh (oy of oY] — B,,, =O 
XIF (9, of emteVE_. ift(g,, gt )eitaVPow] \ 
: 2 (11.32) 
¢ fan 
E Sie ay 
= a! } Inka V + p® “R- 
X {— Has, San Let (ae of pees 
— ley ge | + E,,,cos9 29 2s 
x (f° (A ?; ) em iha Vi +e __ if (Pas a) eitaV tee aaa | } 
if we use the identities 
sta (J, —9)__cos(t, — 9) __cos # sin (2,—9) 
ein Be Sinta, > COS% t¢9,= agate 8 
cos (+, —¢) 1 sin (9 — ¢os(?;—9) (11.33) 
sin sina, €O5 94 at aia cos 2, (B= ein? ea e- 


The operations carried out above may be briefly summarized in 
the following way. 


The field from the nonuniform part of the current 
on the disk 


aAnast 
Se 


| “fal (. He d 


is found (without direct calculation o 


9 the current) in terms of the 
known field of an auxiliary haif-plane 
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baja, e108 a, +9, Ole aes = 8 


ant, : [sene?*as 
-@e 
by a replacement of 


$3e. nerapny Jsenedy 


in those cases when $ 8 1° The functions ¢ and o, are determined 
by the equations 


@=sin sin» — sin B cos (¥— 9), 
== sho y sing — sin 8 cos (¢—p)y/ Taster ae (11.34) 


Solution (11.32) was determined exactly in this way with kay’a*+y"> 1, 
when, for auxiliary half-plane whose edge touches the rim of the 
disk at the pcints y = 6, » = 7 + 6&6, the phase 4 was equal to ¢. 


A solution to the problem using this method also is possible 
in the case 
\ 


%=1 9= 7: (11.35) 


when 7 Se 0. The direction 6=y, => corresponds to the princi- 
pal maximum of the scattering diagram, and therefore is of special 
interest. Substituting the relationships 


arr Ta ‘ 
nn ikde sin 1 cost " 


I, ae. (E,, 4 -L cus a, tz 9} M,,.) 


ny tkin as a 


sin has 


| (11.36) 
} 


which follow in this case from (11.19) into the equations 


fe 
A, = = al ((1,,c0s y-+/,, sin ¢) dy, 
, (11.37) 
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detail] the fringing fleld in the ineldent plane (=0. gz r) where 
the expressions (11.32) take the form 


Aq a 
By ayo age ale * bipaye™l Se 
with °=>5: &>y1. 
ane 
ikap —dtap, ¢ 
B= Hya Bal faye +p ey (12.01) 
withe=y. 8<T. 
ae ; ian 
ihop oes —hepy CO 
Ea Hy = ag E ual pe" — if Meg 
with p=—-y; 
rs ie —itap ; sot ihap e 
Ey Hy = Fa [— Bt (le ig! (2)e a 
withe=+, &>1. 
Le 
aie ee ar ihap - —t e/* 
Bet, =) Garant ne tet ae a (12.02) 
with p=, b<Y, ies 
tt 
Hor iad ° —ihd Pied 
Ev=H,= i ‘taig (2) ¢*" — ig"(I)e TS 


with @=— >: 


The funet toms rhcry and rh) eorrespend te the is or the auxtlbary 
half-plane -OC#VGa , and the Cunet fons rlcn) and 4 Vay cerrespenad 
to the field of the halfeplane m@S"Soe Tn aecerdanve with 


squat bone Crp.) and (11.29), they ere deterritned ty the oxpress ten: 


oy hr 
hMW=1)—,<"L j= iat a 
siny—sind ° siny —sin® : 
Ty noe EY 
— sin—5-"--c08 —5—- 
eos 8 2 2 
a’ (1) == RU = ae aaa g (1) ==- sing ocane ‘ Cte ney 


(Maat Dem ont imued: on oadext perce.) 


aue sin el —tos a 
hQ)=lO+s5 =a 12)=—-; 


sing — sind 


con 8 goa etek (12.03) 
8'(2)== 8 (2+ ane 8(2=— -—ainqg@a an 
if =F and ®<+ ; and 
+ —y } 
cong sin TT ogg 
Pay=fy— say + sind’ iY=— ara 
Saab sin o+1 tT cos = 
a'(lj=g 80 aay ena a(ly= "Suny sio® , 
se aes (12.04) 
eee sin —! +cos 4 
POQ=(+ spp enb (OW ~ ayy ies 
+7 6—y 
sin ZT 80S a 


td 
g'(2) = a(2)+- ame > £Q2)= sing+sne 


if g=— FZ and 6<>. 
It was mentioned above that when ka >> 1 in the direction 

8=7 (0°<7<55°), P= the field from the nonuniform part of the current 

is negligibly small in comparison with the field from the uniform 

part. Therefore, for the field from the nanuniform part of the 

current one may write, with the help of Bessel functions, the 

following interpolation formulas: with eso 


2 


E,=—H,= fore {i (yy — 2) 48) — 


; eit 
HP O+MOMAO SR 


(12.05) 
Ey Hy = 28" (Let (1)— 2! AO 
— ita (1) +" QO» 
and with p=—> 
Eps Hy a2 FONG) + C1 .a8) 


het ; ela 
HHP +P QO Z-- 
(Pquation cont bidesd ono mest or res, 


etiC-2 3-05 0-71 a) 


£y=H, = 2h 1g) #240 


ft (12.06) 
ie’ (1+ g'(2)1 4, @)}—>z-- 
where 
&= ka (sin 6 — sin 7). 
t= ka (sin © +- sin 7). (12.07) 


These expressions are valid in the region 0<0<—; when 9>1 andi >! 
they change to Equations (12.01) and (12.02), and in the direction 
O=y, c= > they give a field equal to zero. 

Using specific expressions for the functions rl and gi, it is 
not difficult to establish that the total field scattered by the 
disk in view of Equations (10.03) and (10.04) may be represented in 
the following form: 


with @=4y 
E,=--H, = 2171) —f21400— | 


@ft# 
— HM N+INL OE. 


- (12.08) 
Ey =H, =" (le) — 2 2, O— 
: ft 
— Hg (1)+ e240} 
and with I= > 
yey SU 12140) + 
+ EIA OS. 
{ (12.09) 
Epa =F le (2 M404 | 
+ ile +2 (24, @) e~- } 
It is convenient to write these expressions as follows 
A= —- Hy = “52 50, ges | 
(12.39) 


E, a H l= = — (4, yon. | 
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where the functions E and = are determined in the region Osts sy. by 
the equations: 


5, 1) |e — Og AO ten 
. i w= 
sin 5 con tt 


£(8, 7) 
20, 1) a2 j;_ 10 vith pat, (12.11) 
Os ee ar : 


and in the region + <f<s 


he aC 
cosy ae : 
Fa, 2) + 40 _ se =with g=—s. 
x8, 7) 


(12.12) 


~~ 


Here assuming y = 0, we obtain the previous relationships (9.03) and 
(9.05). 


In the directions 8=y7 and ®8=s—7 (with gus), where the 
Scattering diagram has a principal maximum, it follows from Equations 
(12.11) and (12.12) that 


Si = (1) = — ka cosy (12.13) 
and 


E(e— 1) —E(e--- 7) = — ka cosy. Cet) 


In the direction toward the source (@=~y. 15) the functions 
¥(0) and 35(8) take the values 


£8 
® pea AO (t). (12.15) 


58) 


Here considering fA=:2, we obtain 


E(x): — —5 (x)= ka, (12.16) 
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which corresponds to the physical optics approach [Equation (7.20)]. 
The functions £ and £ allow one to calculate 


e,=earlElt, s,=e0°(E/ (12.17) 


which are the effective scattering surfaces with the E- and H-polar- 
izations of the incident wave. Let us recall that, by definition, 
the effective scattering surface is a quantity equal to 


om dent! Sh (12.18) 


where 


S= 5, RefEW) (12.19) 


which is the energy flux density averaged over one oscillation cycle 
{the Poynting vector) in the scattered wave, and So is a similar 
quantity for the incident wave. 


In this way, we obtained the expressions for the fringing field 
which approximately take into account the nonuniform part of the 
current. In the incident plane (9== 3) , they have a form which is 
rather simple and convenient for calculations. It is also interesting 
that in this case they satisfy the reciprocity principle as distinct 
from expressions (10.03) and (10.04) which correspond to the uniform 
part of the current. It is not difficult to prove this by verifying 
that Equations (12.11) are not changed with the simultaneous replace- 
ment of y by 9 and of ® by y, and Equations (12.12) are nct chanred 
with the replacement of 3 by 7m - y and of y by s—@0 [in the case 


p=—-y |. 


However, Equations (12.11) and (12.12) lead to a discontinuity 


of the magnetic fleld tangential component Hy on the plane z = 0 in 
which the disk lies. As in the case of diffraction hy a strip, the 
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reason for this is that we did not consider the interaction of the 
edges. It is also necessary to take into account this interaction 
in the case of glancing incidence of the plane wave (T=), when 

the fringing field components E, and yy must be equal to zero. 

Let us point out once again in conclusion to this section that 
expressions (12.11) and (12.12) near directions 0=7, 0= e~—1 (with 
225) have an interpolation character, but in return they allow 
one to represent the fringing field in the incident plane x 2 0 in 
a convenient (uniform) form which frequently is of greatest importance 
(compare § 24). 
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CHAPTER III 


DIFFRACTION BY A FINITE LENGTH CYLINDER 


The distinctive feature of this probiem is that, in addition 
to the nonuniform part of the current on the cylinder's surface 
which is caused by the discontinuity, there also exists a nonuniform 
part of the current arising as a consequence of the smooth curve of 
the surface. This part of the current has the character of waves 
travelling over the cylindrical surface along geodesic lines [36] — 
that is, along spirals on the cylinder. These waves, which as they 
move strixe the edge of the cylinder, undergo diffraction and 
excite secondary surface currents. In turn, the nonuniform part of 
the current resulting from the discontinuity undergoes diffraction 
while being propagated over the cylindrical surface. It is clear 
that specific consideration of all these effects is a verv comriicated 
problem. 


However, if all the linear dimensions of the cylinder are 
sufficiently large in comparison with the wavelenrth, these effects 
may be neglected when calculating the fringing field in manv cases 
which are of practical interest. In particular, thev may be neprlected 
when calculating the field scattered in the direction toward the 
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source [5, 37]. In this case it is sufficient to consider only the 
nonuniform part of the current which is caused by the discontinuity 
of the surface, and we will do this in this Chapter. The equations 
obtained in this way are generalized to the case when the observation 
direction does not colncide with the direction of the source. 


§ 13. The Physical Optics Approach 
Let us investigate the diffraction of plane electromagnetic wave 


E= E, eitivete T+evoe 7 ( 1 3 .01 ) 


on a finite, ideally conducting cylinder of radius a and length Zz. Let 
us position the spherical coordinate system in such a way that its 
origin is at the center of the cylinder, and the normal n to the 
incident wave front lies in the half-plane 9= 5 and forms an angle 

"7 (0<1<+) with the z axis (Figure 24). 


An incident wave having an 
arbitrary linear polarizaton always 
may be represented as the sum of 
two waves with mutually perpen- 
dicular polarizations. Therefore, 
for a complete solution of the 
problem, it is sufficient to in- 


Figure 24. Diffraction of a plane vestigate two particular tases of 
wave by a finite cylinder. n is ; 
the normal to the incident wave 
front. 


incident wave polarization: 


(1) E-polarization, when the incident wave electric vector ts 


perpendicular to the plane (Eo _L yoz)and 
(2) H-polarization, when HoLyoz . 


The uniform part of the current excited on the evlindrtcal sur- 
face by wave (13.01) has, with the F-polarization, the comrene it: 
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f=— = E,.'sinysinge™, 
f= £ Evasinzcosge™, (13.02) 
\ f =< E,,-c0s 1 cos ge” ,” 
and with the H-polarization it has the components 


=f =0, | 
i =—Z- asingel’®, | (13.03) 
where 


‘b= asinysiny + € cosy. (13.04) 


Let us calculate the fieia created by these currents in the region 
| 
‘3-7. 


The vector potential of the fringing field is determined by 
the equations 


os r d¥ 
A=ildy { PG. 0) with y=0 (13.05) 
and 
“ 
Par mer elitr 
ae | dy Me % y) Sal with y>0, (13.06) 
7 
where 
r=VP FOV FeRe ee 


Since the field in the far zone (R> kat, R>kM) Ais of interest to us, 
these expressions may be simplified by uSing the relationship 


r=R+ua sin) sins ~€ cos 8. (13.08) 
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As a result, we obtain a simpler equation 


U 
r ‘ 
an 
Am ES [ eteuneains yy J 2, weer ar, (13.09) 
t 
a 


Since the current components are described by the functions 
i@e"* | then the problem of finding the field reduces 
to a calculation of integrals of the type 


essentially 
a 
] 
f eSeerr—e Mag j F(pyel? %dy= 
ca bd 
2 
t ‘ # (cos T—Cos @) —_ 
= TeimH | (13.10) 
it 


a 
a FE Tier dy 


The integral 


8 Z 
fice tuy, p= ka (sing sin) 


(13.11) 
when p >> 1 is easily calculated by the Stationary nhase method. The 
stationary phase point is determined from the condition isin y ==0 
and equals 

eg 
bs z° (13.12) 
Then assuming g=— +s, we find 
a 
r r 
J raper raya rye’ [eoF a= 
= -$ (13.13) 
: | I~ —ipsice 
; 2n ee 
=P Trier | rat=V Fre 


—-@ 


As a result, we obtain the following expressions for the 
vector potential: 
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with E~polarization 
@ . git 
Ar zy Eusiny 3 -/, A, =A, =0 (13.14) 


and with H-polarization 


Ava He end, Ae =A, = 0s (13.15) 


where 


e 
1 eos J—toe 8) “i # (cus J—cos 6) 
e —® 


i= ik (coa,y — von 0) x 


——+— athe (in pool 8) +d > (13.16) 
xV ka (sta y + in 8) . ‘ 


The fringing field in the region 93> is determined by the 
relationships 


E,=— Hy =ikA,, 
Ey =H, = —ikA,sin 8. (13.17) 
Therefore, with the E-polarization it equals 
ika ,, ..  @ft® 
E,=—H,= 7 Enssint G1] (13.18) 


£,=H,=0, I. 
and with H-polarization 


ine 
Ey =H, = — Fo Mae sind 1, 
B= H,=0. 


The resulting equations show that the field scattered by the cylin- 
drical surface is created mainly by a luminous band adfacent to the 
cylinder's generatrix with q=4=—> . The radiation from this 

band may be represented [see Equation (13.16)] in the form of soheri- 


cal waves diverging from its ends (points 2 and 3 in FPirure 24). 


Now let us write Expressions (13.18) and (13.19) in a ferm which 
is most convenient for calculating the effective scatterinr area 
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E=— Not oe, (8, %). 


(13.20) 
E,=H, =F Me 7), (8, 1). 
Here 
Yi =Gsiny, Y= —Gsin’ 
2. =Gsiay, 2 = —Gsind, (13.21) 
and 


Y ‘ 2 
g =? V aarp 6) x 
sia (omr con | tha (ein p+3in 4d (13.22) 
x cos ~ — cos : e 


The index "0" on x and x. means that the field was calculated in 
the physical optics approach (based on the uniform part of the 
current), and the index "c" shows that this fringing field is created 
by a cylindrical surface. The effective scattering area, in accord= 
ance with (12.17), is determined for a cylindrical surface by the 
relationships 


8 
of, 2= ra") S| = 2e*sin’y |G, 
of, =a" [nat sin’ (Gt. (13.23) 
In the direction of the mirror-reflected ray (®=y), we have 
of pet, = halt sin 8 =? sin 8, (1327) 


In the direction toward the source (9 =<—y) , the functions x. 
and > equal ; 
-tika sin 646m 


7° sin} sin(e/ cus 4) . 4 as 
y Tyr =f ist. as : (13.25) 


These expressions are valid if ka sin @>!. It is not difficult 
to see, by means of equations (13.02) - (13.05), that the: frinimini- 
field equals zero if y = 0 and @=a. Thus in the case cf radar (that 
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is, in the direction toward the source) we find an expression for 

the fringing field in the region ka sin #3%Jj anc in the direction 
@=-a1, Naturally the desire arises to write interpolation equations — 
that 1s, equations which would provide a continuous transition from 
the region ka sin8®1 to the direction @2x. Now let us note thit 
the field scattered by a cylinder is comprised of the fields scattered 
by the lateral (cylindrical) surface and the base (end) of the 
cylinder. In the physical optics approach, the fleld scattered by 

the end of the cylinder 1s equivalent to the field scattered by a 
disk. But the field scattered by a disk is described by Bessel 
functions. Therefore, it 1s also advisable to express the field 
scattered by the cylindrical surface in terms of Bessel functions. 

As s result, the desired interpolation equations for the field 
scattered by the cylindrical surface may be represented in the form 


x =-J'=- ain ind (eitteest e MEMO EY (Cy id, (6)], 
(13.26) 
C = 2kasin 8. 


From this it follows that og So = 0 in the direction @=«, and with 


the conditions ka sin 81 we obtain Equations (13.25). 


The field being scattered by the cylinder's end (by the disk), 
in accordance with equalities (10.06) and (10.07), is descrited in 
the physical optics approach ty the equations 


P+ 

2 = =e tat con® | 
?|- 5A je 

2, 

Consequently, the field scattered by the entire surface of the 


cylinder will be determined in the plane e=— > by the equations: 


(13.27) 


i eR oe 

Bay FE EDO | 
7 eth? rye 

£,= oop oe ry (9), | 


~~ 
jo 
wy) 
e 
rt 
™~ 
~ 


where 
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(el ae el cus * lJ, (€) _ iJ, (6)] ~ 


° 2 cosh 
5” (13.29) 


we SOE J, cc eltteov e_ Qrasin 9). 


Tnese equations allow one to determine in the physical optics 
approach the effective scattering area of a finite cylinder. 


§ 14. The Field Created by the Nonuniform 
Part of the Current 


Let us find the field from the nonuniform part of the current 
caused by the surface's discontinuity. Figuratively speaking, the 
field scattered by the cylinder is created by the "luminous" regions 
cn its end and lateral surface. Mathematically this field is 
described by the sum of spherical waves from the "luminous" points 1, 
2 and 3 (see Figure 24). Obviously the field from the nonuniform 
part of the current also will have the form of spherical waves diverg- 
ing from these same points. 


In the case when the length and diameter of the cylinder are 
sufficiently large in comparison with the wavelength, one may 
approximately consider that the nonuniform part of the current near 
the discontinuity is the same as that on a corresponding wedge. The 
field radiated by this part of the current in principle may be found 
in the same way as in the case of the disk. However, such a method 
is rather complicated. We will find the desired field by a simpler 
and more graphic method, starting from a physical analysis of the 
solution obtained for the disk. 


For this purpose, let us investigate the structure of waves 
(12.01) and (12.02) which are radiated by the disk. These equations 
include the factor 

—— ia {aR eth rae 
V2zka(siny + sia») 


e 
Rae VR appease (14.01) 
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Here = Is the unfoldin;s coefficient of the wave. It shows how 
the fleld is formed with increasing distance from the disk: the 
agiffracted wave which is cylindrical near the disk unfolds [tntn a 
spherical wave as the distance from it increases. The coeffictent 
(sing -{-sin 8)—"* is proportional to the width of the luminous resion on 
the disk or, in other words, to the width of the first Fresnel zone. 
Thus, in Equaticns (12.01) and (12.12) the functions rt and at depend 
only on the body's geometry — more precisely, on the character of 
the discontinuity. 


Therefore, it is entirely natural to assume that the similar 
waves which are being scattered by a cylinder have the same structure 
and differ only in the functions ri and gt which correspond in this 
case to a rectangular wedge. Consequently, in the direction toward 
the source, the field from a nonuniform part of the current flowing 
on the cylinder may be represcnted when ka sin 4%1 in the following 
way: 


E, = —}f ait Eutraye | (:-*) + (tl cos 8 


—{f? (2) efter 1 i Bye eo (tj heer Paid : 

] 

ie t—<) siatcose (14.02) 

E,=H,= oe He | ia" sae 8 z)¢ 
—{e' (2) efter Ot (3) eto Og et) ‘ ee 


In accordance with § 4, the functions rt and gl 


the equations 


are determined by 


ia) =" ( , 
¢ 


= + 
g(t) ‘ oo 1 
es, } 4. coae 
[eS ee — ind? BY 
ee — coe” =) bass (14.03) 


f'(2) ‘ee ' - 
== —~—S i are ae ina 
g* (2) sen — 1 cos — — cos 
(Equation eontinued on next rare.) 
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_. cos8 __ sion & 
“sind * 2cos8 * 


aS a (a 2 


e 
. 
| SSP ES CD 


oe cost 
Sere ete ae ) = ast (14.03) 
ee eng abet | Fens * 
cos — 7 — eos n 
_where 
3 
ae (14.04) 


In Chapter IV, we will show [see Equation (17.25)] that in the 
direction #=2—y=e# one may neglect the field from the nonuniform 
part of the current flowing on the cylinder in compartson with the 
field from the uniform part, tfka>l. Therefore, for the field 
from the nonuniform part of the current, one may write with the help 
of Bessel functions the following interpolation equations: 


TT 
E,=—H=7 bey) (5), 
(14.05) 
E,=H = Ne cay (8). 


F (0) = 141, (+i, eo? — 
—f(3){4, is, (et ® 

3 SBE O Mile = (14,08) 
— g' (3), —1, (Qype ®, 


and the functions mM, v2 and wm ui respectively equal 


M 


myst =f" (2), ae = (1) =F g' @, (14.07) 


= =| 


or 
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® 

mL sia-y 7 ' 

M (moa Fae a = 
a 


cos > — cos 


cos 67 sin? 


2 x a) sind — Zcose ° 


v y= 2(—t—= (14.08) 
n s 


The resulting Equations (14.05) change when ka sin®>I into Equations 
(14.02), and in the direction §=* they give a value equal to zero 
for the field. \ 


§ 15. The Total Fringing Field 


Summing Expressions (13.28) and (14.05), it is not difficult to 
see that the total field scattered by a cylinder will equal 


amo 
E,= —H, = $ Eu GY. 

th OL 
E,=H. a Yo. (15 ) 


where 


fai cos t 
J) = (41, © + iN Oe 


’ 
a (3) [4, (f) —f. (9)| en itces @ r) | 
y, (0) = [MJ, (8) + INS, (8)] eiticos . | 

r 


— g (3) J, (5) — iJ, (6)} en iticos 8 
C= 2kasind ; 


(15.02) 


and the functions M, N and M, N are expressed only in terms of the 


functions f and g which correspond to the asymptotic solution for a 
rectangular wedge 
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=! =) 


p=/M=1e) yprevee (15.03) 


or 


| sin ry =m } 
+ | ae a (= fn a—2— Py a6) ° 
coe 7 n coo cos” 


(15.04) 


The functions f(3) and ¢g(3) in turn are determined by the equation 


‘}(3) eae ——- ‘ , 
gQ3)J 8 ( ee, =a) (15.05) 
. a 


~ =} cos ——— 
cosZ—) com : 


Thus only the functions f and g are included in the final expression: 
for the scattering characteristic of a plane wave by a cylinder. 


In the direction 8=2, the functions 2(6) and s(%) take, as in 
the case of a disk, the values 


£ (2)=: — 5 (2) == koe, (15.06) 
and with %=:35 they respectively equal 
2 2 

sok a Sia ; ' 4 ; 
V(4)= - ~~ —+ pele: tn!) 0) iSO), 

cos — I ° 

4 2 

~~ sin — (15.07) 
mH=(2 ——" 1 cig £ in 

cos— — 1 


+(4 ctg + itl) 4,0, 


where ¢ 2 2ka. The terms in this equation which contain the factor 
kiZ refer to the field from the uniform part of the current, and the 
remaining terms refer to the field from the nonuniform part cf the 
current. 
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In accordance with (12.17), the effective seattering area of 
the cylinder is determined with the E-polarization of the incident 
wave by the function 


3, == na" |B (8)|", (15.08) 
and with the H-polarization of the incident wave by the function. 
3, = #a'| &(8)|*. (15.09) 


Let us note that Fxpressfons (15.02) for the seattering fleld 
may be obtained directly on the basis of an analogy with the 
Equations (12.06), omitting the calculation of the fields from the 
uniform and nonuniform parts of the current. In the same way, one 
may obtain the expressions 


_ al 
i= 


F(®, 0) = (MJ, (+ INS, (he? 


1 Fi teen 0 + con Oy 


—f@(AO—MOle 


(cos © + cos &) 


a (15.10) 
£(0, 0) =[M/,()-+iN, Jer 


. i ws 
— 214, — Wye Fr 


which are suitable for calculating the fringing field in the region 
2 a 
Sm 7 <h OSs 8, =2—-1), The quantities here equal 


t= ka (sin 8 + sin 9,) (152142) 
2 ) es 
wo s 
= t 
a _y n—0—0, -- ® za 
cos ~~ cos Py cos -~—"Ccos a 
fa n 
N ae 2 a (15.12) 
nya" (oactece 
, cos = — cos — 
cos = —co8 aaah os — —cas eh): 1 
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—~ 
ra 


bel) 


s 
f(3) a ' = t 
=— |.) UO ee 
a(3) Ja ( re ane Ee ae) 
Expressions (15.10) satisfy tne reciprocity principle — tiat 
ls, they do not change tneir values if one interchanres ® and 4%. 
When 6=8, they change into tne previous Expressions (15.02). 


Equations (15.02) and (15.10) describe the radiation from the 
currents flowing’ only on part of the cylinder's surface: on whe ane 
end (when z 8 a and on half of the lateral surface (-" < y» < 0). 
Moreover, these expressions do not take into account the nonuniform 
part of the current caused by the curvature of the cylindrical sur- 
face. Therefore, they must be refined with values of 8 and 4 


which are close to = and n. However, in the case 4 =f — that is, in 
the direction towards the source — these corrections may be nerlected 
if the parameters ka and kl are sufficiently large. ‘Numerical cal- 


culations performed by us on the basis cf Equations (15.02) show 
tnat this evidently may te done already when ka = 1 and k? = 10 7. 
The graphs of the functions A, =, E(v)}? and oh aj Ei) * constructed for 
this case in Firures 25 and 26 agree with the experimental curve (1) 
(the dashed line): the position of the maxima and minima basically 
agree, and the number of diffraction fringes is the same. For the 
purpose of illustrating the effect of the ends, we constructed a eraph 
of the effective scattering area for those same values of ka and kl 
taking into account only the uniform part of the current on the ‘ 
cylindrical surface (Figure 27). A comparison of Firures 25, @c and 
27 shows that tne effect of the ends begins to aprear when & = 120°, 


Footnote (1) appears on page 89. 


PID-HC-2 3-259-71 86 


© wo 
pe) 
a tt Ge 
whe SC 6 & EG 
GO VU ey O Uet 
Yu ® puyn 
gx t} Fp ay 
Taga Gd cad d 
& oO © 4% oO oO 
Bo ke : = a teow & 
«gue gc Oosg 
—= EO £0 HOM ca 
eco a ~_A> Ba 
oO AY 
YY o ou ial 
“U8 en om «nN 
IN vey ed a a ke Mee 
Nee oh 1H & w fF 
belie IS Al asdd 
Un EA ooned 
= nae fy “a O 
30 &&eo. 3 oa Wed Gc, 
Vo Sh) worocnme 
oad OL) Apa OT 
gk, fe 


+h 
eet LLL ooo eo 


5 ees 
ESanaem 


a 
ci 
CC 
» 
c 
@ 
E: 
ad 
f. 
ue 
ie 
at 
a 


fH 


FTD-HC~23-259-71 


eee 


Figure 27. The effective scatterins area 
of the lateral surface of a cylinder in 
the physical cptics approach [see (13.26)]. 
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l. on pare 86. The experimental curves shown in Figures 25 and 26 
and also those in Figures 31, 32, 65 and 71 were 
obtained by Ye. Il. Mayzels and L. &. Chugunova. 
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CHAPTER” EY: 


in this Chapte: we will refine the physical crties arrreach 
for certain ctheriodies of rotation, whese surfaces have ciroular 
discontinuities. We will limit curselves tc the case when a plane 


electremarnetic wave falls on the todfes alens thety symmetry axis. 


az before, we will assume that the linesr dimensions cf the 


=] 
ft 
a 
b+ 
%, 
Q 
a) 
f 
ob 
t 
7] 


bedies are large in ccemparison with the wavelenrth. I: 
r 


ity cf a circular disecntinut 


*3 
i) 
>) 
ct 
Us 
b* 
3 
ct 
a 
< 
ie 
iP) 
poe 
a 


2 


nea co 
face of rctaticn may be aprreximately considered tc be the sare 
as that on a ccrresponding conical tedy. Coensequentiv, it 
ficlent to study the field from the nonuniform part of the current 
f 


caused ty the circular disecntinulty of the sur 


§ 16. The Field Created ty the Nonuntfern 


Part of the Current 


Let a plane electromagnetic wave fall on a eonical bedy in the 
resitive direction cf the z axis (Fieure 28). Srem the relotienshins 


E=—-\ (graddivA+#A), | 


‘ (16.01) 
==rotA | 
we find the following expression. fer the fringing field in the wave 


Zone: 


== Hy =ikA,, ee 8—0 : 
Emo Hesibhs oO (16.02) 
and 
£E,=—H, =ikA,, 1 a= 
E,=H,=ikA, Juttn > (16.03) 


The vector potential is determined by the equation 
Asst aa j[fas yer™ cus “(a —Ssinw) do 
(1€.94) 
Dita aval 


Here ris the distance from the discontinuity to the observation 
point, J405) is the surface current density flowing on the trradiated 
side of the body, and Jos) is the current density on the shadewed 
side. The upper sign in the exponents refers to the case 4=e, and 
the lower sign refers to the case 68=0. Cince the nonuniform rart cf 
the current is concentrated mainly in the vicinity of the discontin- 
uity, the vectcr potential ccrrespcending to it may be represented in 
the form 

ce @ 

A= sei ((( here rat 
694 


Ld 


(16.05) 


=a 0) eo? it cos “dt) dy. 


Fhe Z 


Obviously the nonuniform part of the current near the discontin- 
ulty of a conical surface may be considered to be approximately the 
same aS on a corresponding wedye (Figure 29). In the local cvlindri- 


eal coordinate system Ti» o> Zi» the field from the nonuniform part 
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Figure 29, The dthedral ante 
correrpondings to the diseont ine 


Fieure 28. Piteractton of a uity ef a confteal surface. 


plane wave by a contecal body. 
The plane wave Is propayated 


of the rrent flo - oon sue 
alone the z axts, f current Plowing ¢ teh 


a wedre fs determined In the 
far vone (kr) >> 1) by the frollowtny equations: 


(2 2 


Moy Beg itd vn | (16.06) 


where 


a 
A= a on Sate ‘ PRL cutee teh pp 
t eae [| J, ints i+ 

: . (16.07) 

= ’ : ote. bey | 
+f Sere us}e . 

a 
Here the upper aten in the exponents: refers to the case g) som tow, 
and the lower styn — to the case $) 3 wu. On the other hand, In 


§ 4 it wan shown that this fleld equals 


s 
~ 
es 
~ 
_ 
So ee ee ee 


where &,. (2). H,,1) aro the values of the tnetdent wave amp dttude oat 
a Ade 
wd 1 * : 
the wedte edge, and ff and et are angular funet fons eharnetertodas 


the seattertoy diagram. 


Let us introduce the designation 


ir : * it cos 
a Heyer rae +f Hy Ge** br (16.09) 


Equating Expressions (16.06) and (16.08), we find 


J w= foe) 4 


__ Hg, W) 
he tka jaa 


= Tada 8" (16.10) 
The components Jz) and Joy are mutually perpendicular, and when 

®=0 and ®=s they are parallel to the plane xOy (Figure 30). The 
different orientation of the unit vector 4, when 8=0 and 8=2 is 
connected with the fact that the angle $4 is measured from the irrad- 
lated face of the wedge. In the original x, y, z coordinate system 
the vector J has the components 


J,=/, siny—J, cos¢, | a 
Jy=—J, cosy—J, siny | e 


and 


J,= J, sing+J, cosy, 


=e. 
J,=— J, cos¢+- J, sing just . (16.12) 


Substituting Expressions (16.10) here, we obtain 


Ja = 7EGg (P Ee, (9) si0 } — Bg, (#)c08 4), 


—f€ , with $=0 
Sy = Fqe UE ne, (9) 003 9 +.8'H,, (9) si0 9} : 


(16.13) 


and 


I= Tag Eg, (Hsin? + a'H,,, (9) cos 4], 


J, = ie TPE ae, So aH, ip)sin 4} 


t 


en Doak. (16.14) 
! 


Now identifying the current near the conical surface discontinulty 
with the current on the wedre, we find the components of vector 
potential (16.05) 
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A, = rn ah “ue, fP)siny — 
My, (cos sabi, 


] 

‘ 

7 |. =0 (16.15) 
Ase a “P(E, (oreosp | 
t 
d 


He, (9) sin a] chp 
Figure 30. The relative 


orientation of the and 

unit vectors C6) and 

ez) in the cases @=0 

and dO=m., A, =: aie wfure ve, (7) sind + 


+ BH, (9) cos spats, 
with ®==® (1616) 


Aye. —: 98S" fre. cos: - 
c) 
—2'H,, (sin ade, 


Furthermore, let the Plane wave be polarized in such a way that 
Ey! |ox. Then 


Eu, (9) = Eye sing, Hy, (¢) = — Ey, cosy. (16.17) 


Considering these relationships and Substituting Expressions (16.15) 
and (16.16) into Equations (16.02) and (16, 03), we find the fteld 
from the nonuniform part of the current which 1s caused by the 
circular discontinuity of the conical Surface 


Fam Ham 23 (pf + gy 


with ®=0 (16.18) 
F,=H,=0 | ; 
and 
F,-=-H, .- the gt yf 
b nd / R r with Hs:: ¢, (16.19) 
Ey >H,=0 a . 


Equation (16.18) is applicable for the values OSes Fac OrSe , and 
Equation (16.19) for the values 0Se<t O<e". In the case of 4 


94 
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ad 


disk(»-: +. 8 ), the fleld from tie nonuniform part of the current 


equals zero on the 7 axte, cinee e! = --} = -l1/2 when a=, and 
re wt me al/o whee @ ae fe ocpure (F.le)). 

Using the resulting cf. oto ot nti: fs the followine seetions,we 
will calculate tue Pf te ‘torte. area (in the direction f=ae) 
for specifie bidfes. oa ‘hut they are lrradfated by 


the plane wive 
Foe We Ee (16.20) 


and their linear dimenctons are lvree fn comparison with the wave- 
‘length. 


§ 17. A Cone 


Let a cone (Figure 28) be trradiated by plane electromagnetic 
wave (16.20). The uniform part of the current which is excited on 
the cone's surface has the components 

f= x £,,sinwe'™, ' 


(17.01) 


= ae Een conmcenget 
and creates in the direction % =x (with Ra’, R>RP) the fleld 
’ i A et 
£,>—H,= — Ewa tt e-a+ 


+ Eva (4 tate S tga) rem, (17.02) 
F,=-H,=0. 


Here the first term describes the spherical wave divereing from the 
vertex of the cone, and the remaining terms describe the spherteal 
wave from its base. 


Tne field caused by the discontinuity of the surface at tne 


cone base is a soherical wave, and is determined tn accordance wit! 
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(16.19) by the expression 


Sot aR 

a : e" ' 

2=—-H,=— $F teot Ta ic ne (17.03) 
] 
| 
J 


cos a cos 7. 


£,=H,=0, 


where 
nai4et8, (17.04) 


An asymptotic calculaticn of the rigorous diffraction series 
for a semi-infinite cone [38-40] shows that in the direction@#=a 
one may neglect the effect of the nonuniform part of the current 
caused by the conical point. Therefore, summing (17.02) and (17.03), 
we obtain the following expression for the fringing field: 


E,=—H,=—Eu [: tgte(1 — ean 


ge (17.05) 


Let us point out the following important feature of the resulting 
equation. In the problems which were investirated in the previous 
chapters, the edge waves of the fringing field were expressed only 
in terms of the functions f and -. But now inthe equatlon for the 
spherical wave from the cone's base, in addition to the term which 
depends on f and ¢g [the last term in the bracket of Fquation (17.05)], 
there is an additional term [term «-é/stg’me™”’ {n Equation (17.05) ] 
which does not depend on these functions and ts determined by the 
uniform part of the current. Therefore, it is impossible to rerresent 
the resulting spherical wave from the cone's base only in terms of 
the functions f and g which characterize the total edre wave diarram 
from the corresponding wedge edge. This important fact was not 
considered in [41, 44], as a consequence of which their authors did 
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not succeed in obtaining correct results for a cone with an arbitrary 
aperture angle w(0 < w « #/2). 


The effective scattering area in accordance with (12.18) is 
determined by the equation 


a= ea E', (17.06) 


where the function £ is connected with the fringing field by the 
relationship 


Se ee ae es 
= y= — g ex RB (17.07) 
and equals 
| eS 
Ea Ligtesin kf emp —__*_ uns, (17.08) 
cos C08 


The analogous function in the physical optics approach may be written 
in accordance with (17.02) in the form 


EP as gt tytesin kel! — tgaet™, (17.09) 


With the deforming of the top part of the cone into a disk 
2 
(o~ F. (-0) Equations (17.08) and (17.09) are transformed, respec- 
tively, to the form 
f= —tka—— cig, | 


DP? =s — iha, 


(17.10) 


Furthermore, it follows from (17.08) and (17.09) that for larre 
values of the parameter ka(ka >> teow) the functicns © and a may 
be represented in the form 


z sin (17.17) 
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2° = tg w e7""!, (17.12) 


Thus even in the case of short waves ( ka >> tes but R >> K2*), 
our Expression (17.08) does not change into the physical optics 
equation, but substantially differs from it because 


(17.13) 


and 
3° == eat tg’. (17.14) 


With this 


(«co Fem) ge | (17.15) 


that is, for sufficiently short waves (or for sufficiently large 
dimensions of the cone) the function o is proportional to o°, The 
coefficient of proportionality here does not depend on the cone 


dimensions, but. 1s determined only by its shape. 


This result is graphically illustrated by the curves giving 
the effective scattering area of a cone (w 2 10°25', k = nr, 2 = 90°) 
as a function of its length (Figure 31). Whereas our equation (the 
continuous line) is in satisfactory agreement with the results of 
measurements (the small crosses)‘1) , the physical optics approach 
(the dashed line) gives values which are smaller than the experimental 
values by 13-15 dB. For sharply pointed cones, the nonuniform part of 
the current has an especially large values. In Figure 32, a curve is 


constructed for the effective surface of a cone (ka = 2.75 1, 2 = 90°) 


Footnote (1) appears on page 113. 
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with its deformation into a disk 
(w + 90°). The discrepancy be- 
tween our curve and the physical 
optics approach here reaches 
almost 30 dB when w = 2°. 


Expression (17.08) obtained 
by us also allows one, in con- 
trast to the physical aptics 
approach (17.09), to evaluate 
the role of the shape of the 
shadowed part of the body and 
Shows that the reflected signal 


Figure 31. The effective scatter- Will be larger, the closer this 


ing area of a finite cone as a Shape is to a funnel-shaped form 
function of its length. The 

function 3 (the continuous line) (Q =n - w). Thus, for example, 
was calculated on the basis of in the case w 2 10°, ki = 10 r 
equation (17.06) which con- 

siders the nonuniform part of (k = 7) the sipnal reflected by 
the current in the vicinity of the cone may exceed by 15 dE the 
the circular disccntinuity. 4 : 

The function 0° (the dashed value corresponding to physical 
line) correspcnds to the physi- optics (see Figure 33) if 2 = 


cal optics approach. 170°. 


Let us note that our Expression (17.13) 1s equivalent to the 
expression presented in the above-mentioned papers (41, 44]. HBowever, 
the latter expression 1s applicable only for sharply pointed cones, 
whereas we have, in addition to (17.13), Equation (17.08) which is 
suitable for cones with any aperture angle uw (0 Fes re 

The calculation method discussed may be peneralized in the cnse 
of asymmetric irradiation of the cone. However, with asymmetric 
irradiation, generally speaking, it is necessary to take inte account 
the nonuniform part of the current caused by the point of the cere. 


In concluding this section, let us calculate the effective 
scattering area for a body which is formed by rotation arcuni the 


PTD=HC $2 3-259-71 99 


Figure 32. The effective scattering area of a finite cone 
as a function of the vertex angle. 


z axis ef the plane figure shown 
in Figure 34. Integrating the 
uniform part of the current, it is 
not difficult to show that the 
field scattered in the direction 
d=" by the lateral surface of the 
truncated cone (Figure 35) is 
determined by the equation 


TT ENELLSE 
Fy IN TT Tt es 
a 4 
| 
Saqgun sean 
rE EL TT RT 
PET TTT ft 
BRRmRGA ee 


MTTTiltititiiiii tf 
TIT TT Try Ty Try rer 
*“FTTTLITIILeT Tt LL OL 


E,=—H,=E,:: alas tg? @, -}- 3 By) enh as 


Figure 33. The effective scatter~ 


‘i A it R 
ing area of a finite cone as a +( deters A Pig o,) ie qt, bh) Sgt 
function of the shape of the ; 
shaded part. (17.16) 
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Figure 34. The generatrix of Figure 35. The generatrix of 
the surface of rotation. of a truncated conical surface. 


Summing this expressing with (17.02), where the quantities 7 and a 
must be replaced by Ly and a,» we find the field from the uniform 
part of the current flowing on the entire illuminated side of the 
body 


Ly) 
E,=—H,=—% 5" fe 


+[jt tg*e, sin ate 4(1 — eM lig oer \s ‘ 


ty? wsin ki, e'** — tue 4. 
(17.17) 


The field radiated by the nonuniform part of the current is 
determined in accordance with § 16 by the equation 


2 a 
“= sin 
£,= —H,=— el i +tg w—tgee™"+ 
all” Pema 
2 (17.18) 
a, ml 2b (t, +f,) eftt 
— | —-———_--— 4 tzu, Jer orn | —_., 
+2/ con = — eon 5! ie - 
where 
a =1P 2, nyt pts, (17.19) 


Now summing (17.17) and (17.18), we obtain a refined expression 
for the field scattere1 in the direction 8-s 


poate ig? sin kl,eM 4 (17.79) 


(Equation continued on next pase.) 
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. aikt ! 2 i gible + 2:40, 
Pees get "F gat ty?w, sin &l,e"* ‘.f- 
cos a, 98 m 


(L520): 


Z 
— sin— ike 
a, Ae eS, ox gail sn) en: 


Consequently, the effective scattering area will equal 


} ; ial, a a it, 
ia; tg’ wsin ke + ay 


C8 en, 


? 
a== FQ, 


» (17.21) 


dai asikh 4 ae My Me gti thts 


+ . tg?u, sit &/,e 
i] 


In the physical optics approach, the analogous quantity equals 
o? = 2a; | ty2osin kie™” — tg wet 4 
t] 
! . i A Hy A 3 (17.22) 
+[ ia; tetas ki,e* + (1 sar ev" tgw, |e | : 


\ 
When the top part of the cone is deformed into a disk (e+3. +0), 
Equations (17.21) and (17.22) take the form 


2 . i] . i 
¢ ea | —~ika, -- x tgs tye tg, sin kien} 


2 2 . 
~- sin — 17.2 
a oe eS tht, (17.23) 
WN sg heaped 
cos fs — cos My ’ 


¢ = ra; | = ika, tue tg? o, sin kie™ + 


+ (1g) tga [. (17.24) 


In these expressions assuming w, * 0, we find the effective 
‘scattering area for a finite cylinier 
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ai | -- ita, boty 74-2" lth 
3 == -- te ic — ® 
io En ae (17.25) 
3.2 naj (ka,)’, (17.26) 
in connection with which 
Pees 1 
A= A= +e (17.27) 


Equation (17.25) 1s more precise than Equation (13.06) which 
was derived in § 15, where the value of the field in the direction 
Q=2 was taken in the physical optics approach. 


§ 18. A Paraboloid of Rotation 


Let us calculate the effective scattering area of a paraboloid 
of rotation r? = Q2pz (Figure 36) which is irradiated by plane wave 
(16.20). The untform part of the current excited on the parabcloid's 
surface has the components 


f= = E,, sin ae’, 
i= (18.01) 
=F Fas oer J 


Integrating this current, it is not difficult to show that in the 
direction ?=w it radiates the field 


E,= =? Hy, == cs 5a— : re) ty fea 18 o7) 
ie | (107024 


Here a is tne radius of the base of the paraboloid; /==5. = J ctgw is 
its length; a is the angle between the z axis and the tangent to the 
generatrix of the paraboloid (r? = 2pz). At the point z= 7, the 


anzle a takes the value a=e(tge= 2) . 
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The (feld from the nenuntfern part 


by Equation (17.02). The ffeld from the 


to 
aad 


nonuniform part cf the current which 


eaused by the smeotn curve of tne parnae- 
holetdi’s surface eyuals zero in the case 


os 


of symmetric radiation [45]. ‘Therefore, 


Fisure 36. Ciffraction 
of a plane wave by 4 summins (189.02) and (17.03) we find the 
4c Ee ae a ete tier e . 
Parebeioy dO: bern expression for the resu:ting fringing 
field, 
eae 
a&,, ava =gug pet 
sao or tga-- oa emma al Rr e 
. cos 7 — cos 
E,=4,=0 (18.03) 


Consequentiv, the effective scattering area cf the parsbololid 


will be determined by the relationship 


2 1 
= sin -- 
a a 2 ’ ” 


cos 7 — COs 


which, when the paraboloid is deformed int. a disk (o+ 7, 4-0, 9 


const) » + ~ ansfcormed ad the .OPn 
32278 ika -|- ctg -- . (18.05) 
| ‘ 


Comparing Excressi:n (18.54) with the equation 


a*-= ra’ tg*w J} —e |, (16.006) 
which physical optics sives for the effective scatterins area, we 
see that they differ signiflcantly from one another. First of all, 


the oscillating character of the funetion 09 draws our attenticn: 


*| 
ae 
i 
' 
Q 
tv 
w 
i) 
So 
pa 
es) 
1 
~ 
—~ 
— 
oS 
Bi 


of the current eagused by the discentinuity 


ef the parabolciad'ts surface is determined 


ges equals vers §f a whele number of hal fevavec 


1 3! 
t=: A a. ? 3. é 5 . 
+ Ma: sal - --) ts fitted inte the leneti: of the caratc aa and 
{ct takes 2 maximum value If a half-fnte,sral number of hal fewayes 
A 1 
(: sy fat yj 123...) 1: ecntained in this lenrth. 


A calculation performed by us on the basis of Sausticon (1£.04) 
for parabcloids with the parareters 2 = 90°, teu = 9.1 (% = ©) shows 
(Ficure 37) that, although the ccetllatinz character cf the effective 
Scattering area 1s preserved, the amplitude of the oscillations ts 
only 2 dB, and the maximum values of the function o exceed the corres- 
ponding values tn the physical optics approach by almost 13 4B. A 
Still stronger diverger:ze between the results of our theory ind 
Dhysical optics is detected when the paraboloid is deformed into a 
disk (Figure 36, ka = jn , k 2 7, 2 2 90°, w+» 90°). 


Figure 37. The effective scatterine area of a finite 
paraboloid as a function of its length with a constant 
value of the angle w (tgw = 0.1). The diameter of the 
base varies. 
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—— 5s, a lesd 
---43! ie: A 


é 


Pieure 38. The effective seatterinre area cfia 
finite parabolo!ad as a function of ‘ts lenrth 
with a eonstant radius of the Lasse. 


As in the case of a cone, the shape of the shadowed purt turns 
out to have a substantial influence on the retlected sienal. For 
example, for a paraboloid with the paramcters ka = on, KP = 19 a, 
tru = O11 (k 


Increase of Q(w < 8 << mew) (Pieuire 39). 


mn), the reflected sienal trereases by WE di with an 


HYCe 2959-7) pa 


Figure 39. The effective scatter- 
ing area of a finite paraboloid 
as a function of the shape of 
tne shaded part. 


In concluding this section, 
let us dwell for a moment on 
the question of calculating the 
effective scattering area for 
bodies of rotation of a complex 
shape, whose elements arc the 
lateral surfaces of truncated 
The field from the 
nonuniform part of the current 
arising in the vicinity of 
circular discontinuities may be 
determined without difficulty 
from Equation (17.03). The field 
from the uniform part of the 


paraboloids. 


current is found by quadratures. 
Thus, the field being created in 
the direction #=<= by the uniform 


part of the current which flows on the lateral surface of the truncated 
paraboloid r? = 2pz(p = a,tgu, = astgus; see Figure 40) ts determined 


by the equation 


u Peer les 
Ea= — Hy = — F Eye laity oy -— arty oy oe y gti — 


£,=H,=0 
Here 


! 
l= 7 (as ctg wo, — 2, ctg »,) 


- (18.07) 


(18.08) 


is the height of the truncated paraboloid (the distance between its 


bases). 


corollary of Expression (18.02): 


Let us ncte that Equation (18.07) is a simple alrebrnaic 
it is the difference of the fields 


scattered, respectively, by the paraboloid of heirht hth= 5 ant 


2 
by the paraboloid of heirht b= 35° 
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§ 19. A Spherical Surface 


Incident wave (16.20) excites a 
surface current on the surface of an 
ideally conducting sphere (a radius of 

—t , p and a center on the z axis at the 
point z = 9). The uniform part of 


Figure 40. The generatrix of this current has the components 


the lateral surface of a : jus 
truncated paraboloid of p=—s 92 cos Ge, | 
rotation. 
i= (19.01) 
i? =: 5- Eys sin @ cos ye**. 


The currents flowing on a spherical ring cut from the sphere's 
surface by the planes z 2 Ly and z = ty + ls (Figure 41) create, in 
the direction 9==2, the fleld 


s==7, £5 (-(3 tga, =u) eth | 


i\ nj ike 
+(Fte%.— 4) etme, (19.02) 
£,=H,=0, 
where 

f,=p(1 —sine,), } (19.03) 

ft, == (sin @, — sin w,); 

—— _ 41. _ 

P= osu, cosa, * (19.04) 


Here a, is the radius of the first cross section; a, is the radius 
of the second cross section; w, (os) ids the angle between the @ axts 


and the tangent to the meridian at the potnt s = Lj (2,5 zy + fy). 
Furthermore, assuming in Equation (19.02) w, =: 5 (, == const) , we 


obtain in the physical optics approach an expression for the fPleld 


scattered by the spherical segment (Firure 42) 
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Figure 42. A spherical segment 


with a conically shaped base. 
Figure 41. A ring cut from the 


surface of its sphere by the , ; 
planes z 3 Ly and z = Ly + Lo. E,=—He En | —ytctat 
Pr dal 


+(Stg0—f) eM]. (19.05) 


Here we used the new designations 


\ 
@=4,, Sy 


= 1, =9(1—sine). (19.06) 


Equations (19.02) and (19.05) are simplified if ka, >> 1 and 
ka, >> 1. Thus, the field from the spherical ring will equal 


y ry (aR 
E, =—H, = —£,.(% tg e, ~—3 tga, em) on , 


(19.07) 
and the field from the spherical segment will equal 
arp elt 
E,=— Hy =: — "7 (sec a—tg we). (19.98) 
If here one assumes w 2 0, then equation 
aF,, e* .69 


gives us the field scattered by a hemisphere. The value of the 
effective scattering area corresponding to it will equal, in 
accordance with (17.06), 


of = za", 


(19.10) 
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Now let us find the field scattered by the spherical segment 
considering the discontinuity of the surface; one mav neglect the 
perturbation of the current. as a consequence of the smooth curve of 
the surface if ka >> 1 [74]. The nonuniform part of the current 
which is caused by the discontinuity creates in the direction 0= 
the field (17.03). Summing the latter with the field (19.08), we 
find the desired field 


2 2 
= aE yaf 1 Tar) lth 
Bama taste r° (19.12) 
Os — cos a. 
Consequently, the effective scattering area of a spherical segment 
will equal 
= in * ; 
t a 7 2iht 
o=2q' tos @ ee oe te . 
cos-, — cos a (19.12) 


n=142t8, 


In the physical optics approach, a similar quantity is determined 
by field (19.08) and equals 


#— xa*|-!  —tgwe", 
(19.13) 


cos @ 


With the deforming of the spherical surface into a disk 


a? 


(2+ F..!-+0, Q=sconst ) Equations (19.12) and (19.13) are transformed, 
respectively, to the form 


) 


1 r 3 
=: ra? |ika + — ctg — 
a + ct 5 | (19.14) 


s¢=- ra’ (ka)’. 


Tt follows from Equations (19.12) and (19.13) that the effective 
scattering area of a spherical segment is an oscillating function of 
its length. The oscillation period equals 3. Numerical calculations 
rerformed on the basis of these equations showed (Ficure 43) that, 
with small angles of the discontinuity (2 = 15°), ene may still 
neglect the field from the nonuniform part of the current. In Fisure 
44, graphs are constructed for the effective scatterinr area of a 
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Figure 43. The effective scatter- Figure 44. A comparison of the 


ing area of a spherical segment effective scattering area of 

as a function of its length with a spherical segment (continuous 
a constant radius of the base. line) and a finite cone (dashed 
The function o (the continuous line) which have the same bases. 


line) is calculated on the basis 
of Equation (19.12) which con- 
siders the nonuniform part of 
the current near the discontin- spherical segment and a fintte 
uity. The function o®% (the 


dashed line) is calculated from cone (the dashed curve) which 


Equation (19.13), and corres- have the same diameter and tase 
ponds to the physical optics ehawe 
approach. Pes 


The results obtained in this Chapter show that the reflected 
signal depends suostantially on the shape of the shaded part of the 
body, and increases with an increase of the concavity. However, 
since the nonuniform part of the current is concentrated mainly near 
the discontinuity, that part of the shaded surface which is several 
wavelengths away from the discontinuity evidently will not have a 
noticeable effect on the reflected signal and may be an arbitrary 
shape. 
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It is interesting that our expressions, which agree satisfactor- 
ily with experiments, even with large (in comparison with the wiave- 
lengths) dimensions of the bodies, do not change into the physical 
optics equations, but differ from them substantially. At the same 
time, physical optics, contrary to the widely held opinion concerning 
its reliability in such cases, leads to a significant discrepancy with 
experiments. 


The method used ‘n tis Chapter allows one to calculate the 
effective scattering area associated with the symmetric irradiation 
of any convex body of rotation, the surface of which has circular 
discontinuities. It may also be generalized to the case of asymmetric 
irradiation. However, when doing this it is necessary to take into 
account the nonuniform part of the current caused by the point and 
the smooth curve of the surface. 
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FOOTNOTES 


1. on page 98. See footnote on page 86. 
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CHAPTER V 


SECONDARY DIFFRACTION 


In the previous chapters, an approximate solution of diffraction 
problems was carried out which was based on the representation of the 
fringing fleld in the form of the sum of the fields from the uniform 
and nonuniform parts of the surface current. The first field was 
found by quadratures, and the second field by approximation; it was 
assumed that the nonuniform part of the current near the discontinuity 
(edge) of a surface is the same as on a corresponding wedge. 


However, the fields found by such a method are actually the 
fields from the currents flowing, not only on the flat and curved 
parts of the body's surface, but also to some extent on the geometric 
extension of these sections. The error in the expressions for the 
fringing field which is thus introduced is most significant with a 
glancing incident wave, when the edge zone occupied by the nonuniform 
part of the current is noticeably broadened, and also with a glancing 
radiation, when the direction to the observation point forms a small 
angle with the given section of the surface. In these cases, the 
results obtained earlier are in need of substantial corrections. We 
already talked about this briefly in § 6 and § 12. 
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Por the purpose of refining the solutions which were found 
previously, it is necessary to assume that in actuality the currents 
flow only on the body's surface, and that a wave travelling from one 
edge to the other will undergo a perturbation at the latter. The pro- 
cess of forming the fringing field when this occurs may be investigated 
in the following way. The edge wave propagated from one of the edges 
is diffracted by the other edges; the waves arising with this in turn 
are diffracted by adjacent edges, etc. In this chapter, we will 
investigate the case when the dimensions of the surface faces are so 
large in comparison with the wavelength that it is sufficient to limit 
oneself to considering the diffraction of only the primary edge waves. 
This phenomenon we shall call secondary diffraction. 


In this chapter, secondary diffraction by an infinitely long 
strip (§ 20 - § 23) and by a circular disk (§ 24) is studied. The 
solution of these problems may be obtained by means of the principle 
of duality from the solution of the diffraction problems for an 
infinite slit and a circular hole in a flat, ideally conducting screen. 
In the latter case, the physical treatment of diffraction of edge - 
waves is significantly simpler; it is exactly for this reason, therefore, 
that almost all diffraction studies of edge waves are related to holes 
in a plane screen. However, we will not take such a path, but we 
shall investigate a strip and a disk directly. This approach has the 
advantage that it is easily generalized to the case of three- 
dimensional bodies. 


§ 20. Secondary Diffraction by a Strip. 
Formulation of the Problem. 


Let an infinitely thin, ideally conducting strip of width 2a and 
unlimited length be orientated in space as shown in Figure 45. A 


Plane electromagnetic wave incident normal to the strip'’s edges is 
directed at an angle a to the plane xoz and has the following form: 


E=E,e?" cme + onee H=H,e*" tose + vain) (20. 01) 
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In § 6 approximation expres-~- 
Sions were found for the fringing 
field in the far zone which did 
not consider the interaction of 
the edges. In the case of F- 
polarization of the incident wave 
(E.| 02) , these expressions may be 
represented in the form 


Er=— Hi, = Eu lf (lee one— any 


i{arae 
“PL ape Menke mses e ( ) 
Figure 45. The transverse cross V 2abr 
section of a strip with the E.=H,=0, 
Plane xoy, x = 0, y = a and - 
x = 0, y # -a are the coordi- (20.02) 
nates of the strip's edge; n is 
the normal to the incident plane 
wave front. and in the case of H-polarization 
(H, {| 02) 
Hy= E, = Hyalg (1) eine ~ 2 9) 
2 wi (e+) 
— dae (sine ~ 4 
+a (2)e ‘ | Y dxkr ® 
HM =£,=0. ' (20.03) 


Let us recall that the functions f and @ included here are determined 
in the region (elgg (when lal<>) by the following relationships: 


cos Ft aot? ; cos tty sin 
1] ae (Oe away (20.04) 
a(l)=—f (2), g(2)= —J (I). (20.05) 


The first terms in Equations (20.02) and (20.03) describe cylin- 
drical waves diverging from edge 1 (y = a), and the second terms 
describe the cylindrical waves diverging from edge 2 (y = -a). The 
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nonuniform part of the current on each side of the strip also has the 
form of waves which diverge from edges 1 and 2, and are an "analytical 
extension" of the corresponding terms in Equations (20.02) and (20.03). 
The current wave encountering the opposite edge is reflected from it. 
Or else one may say that each of the cylindrical waves propagated from 


edge 1 or 2 undergoes diffraction by the opposite edge (secondary 
diffraction). 


If the strip's width is sufficiently large in comparison with 
the wavelength, then one may appr«ximately assume that the oncoming 
current wave near the strip's edg: will be the same as on a corre- 
sponding half-plane excited by a linear source, the moment of which 
is selected in a definite way. It is also obvious that the current 
waves reflected from the edge will also coincide. Consequently, the 
problem of secondary diffraction by a strip may be reduced to the 
problem of the diffraction of a cylindrical wave by a half-plane. 

> The field created at the point P by a current filament parallel 
to the half-plane's edge and passing through the point Q (Figure 46) 
may be found by means of the reciprocity principle. In the case of 
E-polarization, it is determined by the relationship 


E,= 2 £.(Q). (20.06) 


and in the case of H-polarization 


Mm, 
Here p, (m,) is the electric (magnetic) moment of the current fila- 
ment passing through the Qq; Pes (m2) is the moment of the auxiliary 
current filament passing through the point P with the coordinates 


(o", R), and H, (Q) or E. (Q) is the field created by the auxiliary 
filament at the point Qq. 


Now let us remove the auxiliary current filament to such a 
distance that the cylindrical wave arriving from it may be considered 
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to be a plane wave on the section 
from the edge of the half-plane 
to the point Q. In this case, in 
accordance with § 1 and § 2 the 
field created by it at the point 
Q will equal 


E,(Q)= Ey2(0) (ald, ¥ — 9") 4 +9"), 
H, (Q)= Her (0) [4 (d. ¥ — 9") +4 ld. 7+P))- 


(20.08) 


Figure 46. Diffraction of a 
cylindrical wave by a half-plane. The functions u introduced here 


Q is the source, Q* is the 
mirror image of the source , are determined (for the values 


and P is the observation point. O%<#"<*) by the equations 


AE Vindeos tS 
u(d4,¢ —?)= eee oe) = x elda+ 
7 ae , 20 cos = 
+ een tno Fete 

9 with ®+9° <% <2, 


e ne Cad +57" 
15 Vikd cos = 


ud y teem fag 


y* 


+{ ender ee thoy <t—F 
0 with r—P <P cls, 


¢! +o" 
2 


(20.09) 


and the quantities E, (0) and He (0) are the values of the primary 
field created by the auxiliary filament at points corresponding to 

the half-plane's edge. In accordance with Equations (1.21) and (1.22), 
this field may be represented when kR >> 1 in the form 


9 


Eu. (0)= ik’ pa YW ee ; | 
! 
Ha (W)= ik’ Gee } 


as (ae) 


(20.10) 
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Consequently, an electric current filament located above an 
ideally conducting half-plane excites, at the point F, the field 


see FRR 
: tee ph ; 2: 
E, = ik" p,|utd. 2 -- 7")-- ald, ¢ + 9")} V xe \ a) (20.11) | 
and a magnetic current filament excites, at the point FP, the field 


, 4 ae 2x : uy 
H,= ik?m,|a(d. 2 -- 9) 18d. TF iy ape ( oo (20.12) 


-It 18 easy to see that the exponent e*-4** =" 4n these expressions 
corresponds to the primary cylindrical wave arriving at the observa- 
tion point P, and the cxponent ¢@itl&-¢-™@' +e" corresponds to the reflec- 
ted cylindrical wave. 


The moments m, and P. must be selected in such a way that in 
the dlrection 9" #2 (Figure 446) the filament would create a field 
equal to the field of the primary edze wave above an infinite, ideally 
conducting plane. We will conclude tnese calculations in the follow- 
ing sections, but for now let us make still one other ecmment on tne 
formulation of the problem. 


In the previous chapters it was shown that the scattering object 
may be approximated by a series cf sources — “lumincus" lines and 
points. Therefore, the problem of seccndary diffraction may be formu- 
lated as a rrotlem of searchinc for functions which describe the 
continuous change of the field of each such s:urce during tl.e passage 
through the boundary of the lisht and shadcw torresponding to the 
source. 


§ 23. Secerncary Diffraction by a Strip (H-Pclarization) 


A current fliament with the moment m_ which ‘ts positioned abeve 
an ideally conducting piane (h = 0, Figure 46) creates in space the 
field 
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Ay aban 2H (RR). (21.01) 


Rar from the fllament (when kh, >> L), thio tleld fa desertbed by the 


Aoympoot le express don 


(ams <) 
Hy Ack Soe, eas 
¥ 2:aR, 
Rut the primary ede wave in the direetion @" = on takes the value 


é aR +) 


Hy £40, (Q) 7 (Q) “Yk, ’ (0b. Od) 


whrevre Ho 6) Is the ffledd of the thetdent plane wave at othe potut oy 
r(Q) ts the value of the angular funet fon of the primary ovdl Podifecst 
wive In the dfreetton towards the opposite edge of the strip.  Equort- 
faye Expressfons (21.02) and (beg), we Cin the (i bament ts moment , 
the CLleld of whleh we use to approximate the petinry ede wave, fu 


the form 


my = ay HaniQH(Q). Crem) 


Aso a opesult, the tleld ereated by the (hlaments leeated ahove 
the hatfeplane —@erynk and correspond ity: toe edie 1 (eM ture yy aay 


be represented Cor proton fap r fro the form 
HA) AF 4-1). hate 


The Punet fon 


Wise (5 - 2) is 
+ | 7. @ ey being. singe 
fin ee 8 ee ye ag ; gl Woe 


. t : 
deseribes: the wave padtated by Che souree om mbt thre Mune Pari 


bt 


PTPoH Ce. f60-7 pia 


“1, ! = ie aie aa (ele amoin ont 
H~(\)= Hug (1) X | ea i ae 
(a oy 


+ Hua(l) ee eltetein e- sia er (21.07) 


describes the wave radiated by 
the source mo The sum of these 
waves equals 


Figure 47. The problem of 
secondary diffraction by a 

strip. ¥ tacoe( 3-4) 

mM» and mM), are the sources, 


the fields of which are used 

when approximating the primary 7 
edge wave being propagated from “ (ee wae : 

edge 1 (y = a); +H g(t: — itu grin a - sin ed 
+ - Y 2zke 

M5, and m5, are the sources, 


the fields from which are used 
when approximating the primary 
wave from edge 2 (y #® <a). (21.08) 


e as in 
~ —— the (si o— sin 5) 
i ae 


H, (== Hag (WX j e'N'dq Y 


The first term in this expression is the Gesired secondary wave from 
edge 2, and the second term is the field radiated by the filament 
which is loacted above the ideally conducting plane x = '0 and has the 
moment 


eas { we tha sina 
mo ce He (eens, (21.09) 


es seat" 


where 


g(ty= eid 


een 
¢a- 9 


(21.10) 


FTD-HC-2 3-259-71 120 


Summing the secondary wave which has been found with the unper- 
turbed primary wave from edge 1; we obtain 


H,(1—-2)= 
Wien (7-5) 
oo ihe ‘ 
= Hy 2 8) x { e'dq veer 
4 
i{ars = 


Ve a ina (oln o— sing) 
Ae Va (21.11) 


This expression reduces to zero if one assumes $ = = 1/2; consequently, 
the secondary diffraction eliminates the field discontinuities which 
occurred in the previous approximation when $¢ # = 1/2. However, in 
the direction $ # 1/2 the field (21.11) is different from zero. Since 
H, is an odd function of the x coordinate, the relationship 

Al = nO means that the fringing field components H, and E9 will 
undergo a discontinuity with a transition through the direction 

¢ = t/2. The reason for such a jump, as before, is that in our calcu- 
lations the plane x = 0 is a plane of currents. By finding the 
secondary wave from edge 2, we actually considered that the diffrac- 
tion takes place not on the edge of a finite width strip, but on the 
edge of an ideally conducting half-plane —acyqm. 


Again the resulting discontinuity has an order of magnitude of 
vaye - It is clear that one may completely eliminate the field 
discontinuities only with consideration of multiple diffraction. 
However, the calculation of fields arising with multiple diffraction 
requires specific consideration of the following terms in order of 
smallness in the expansion of the primary edge in inverse powers of 
ir (see, for example, [46]). All this greatly complicates the cal- 
qilations. Therefore, we, using the condition ka >> 1, will limit 
ourselves to an investigation of secondary diffraction, and in order 
to eliminate the discontinuities in the plane x = 0, we will proceed 
in the following way. 
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Let us consider the quantity #(1) in the Expression (21.11) to 
be a function of the angle ¢ [see Equation (20.05)], that is, let us 
replace @(1) by the function g(1). In this case the equation 


H,(t—2)= 
Wie coe(7—5) 


2 ¥ let ad ike (stn e—sin 9) 
=H,,- = e(t) x \ e dq Yiu ° + 
ea e 
_ (wes 
+ Hae (a ght tine-sing) (21.12) 


will give qualitatively correct results not only when 7 y> but 
also with all other values of ¢. Actually, the Fresnel integral is 
Slose to zero if Vkacos(+--+)>1 » and in Equation (21.12) only the 
second term remains, as must be the case. Therefore, Equation (21.12) 
may be investigated as an interpolation equation, and it may be 
applied with any values of a(iPi<z,. It is easy to establish that 
now the fringing field does not undergo a discontinuity with the 
passage through plane x # 0, since Expression (21.12) becomes zero 
when tS ° 


It is interesting to note that Equation (21.12) automatically 
follows from Equation (21.08) if in the latter equation one replaces 
g(1) by g(1). Essentially, this substitution is equivalent to the 
assumption that the moments of the filaments, the fields of which are 
used for approximating the primary edge waves, depend on the radiation 
direction (that is, on the azimuth ¢ of the observation point) 


er gar Hag else’ (21.13) 


Such a determination of the moments of the auxiliary linear sources 
is used, for example, in the work of Millar [47]. 


Precisely in the same way that Equations (21.06) and (21.07) 
were obtained, we find (when x > 0) 
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ofan t,s6 
Wa cos (++3) 


Ht (2)-= us Hug (2) J eft dy 5 eiterln emsin gy 


2 (21.14) 
Wisco (F+-$) 
°: ef Oe. ieee ) 
Ha (2) = > Hugi?) J ef dg Fae an 
i(e+z) (21.15) 


+ Away (2)° WE ea ite (sin e—sin of 


These expressions give the field created by the filaments which are 
located above the ideally conducting half-plane ~*<y<e@ and have 
the moments 


mo = -- mt = os Hag (2)e7 re, (21.26) 


In accordance with Equation (21.04), here 


S ad 


@(2)=¢(2)f to (21.17) 


Furthermore, summing (21.14) ana (21.15), we obtain 


2V Ra cos (F +f) : 
ear 
Hi, (2)= = Hut z (2) \ “dq rz ew ite (sin e—sin m4 
@ 
i(ar ps 


am e '  miky qata e—<ine) 
t Hug (2)-Fee (21.18) 


Here the first term is the desired secondary wave from edge 1, and 

the second term is the field radiated by the filament which is located 
above the ideally conducting plane x = 0 and has the moment Ms 5° 

Summing the secondary wave which has been found with the unperturbed 
primary wave from edge 2, we have 
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Wiacos (+3) 
i 


= .Y ay : 7 
H,(2— = My, R12) \ e'“dq ee ee 


@e 
: waa —the (vin o—sin gt : 
+ Aug (2) vie ° . (21.19) 


It is not difficult to see that the resulting expression becomes zero 
if one assumes 9= 5 in it. Consequently, the secondary diffraction 
eliminates the field discontinuity which we had earlier (§ 6) when 
=F, but at the same time it leads to a field discontinuity when 
=—-5- Again the resulting field discontinuity may be eliminated 
by the above indicated method, replacing the quantity ¢(2) by g(2) 
— that is, by assuming the moments mM» and n>, depend on the obser- 
vation angle ¢. Actually as a result of such a substitution, we 
obtain from (21.19) the expression 


We coe (T+) 
2 % ie aa —ika (sin e-sin gp 
H,(2—1)= Ha > €(2) \. ede Vime + 


eo : , 
t (a+ > 


+Hug(yaa e 


—tha (sin e~—se ¢) 
. 


(21.20) 


which vanishes when gute - This expression may be investigated as 
an interpolation equation which describes the field created in the 
region lel<-5 by the primary wave of edge 2 with consideration of 
its diffraction at edge 1. 


Now summing (21.12) and (21.20), we obtain the following expres- 
sion for the total field scattered by the strip: 


H, = Hq[G(1, pg (ier + 
1(ae + +) 
+G(2 gg (apeMrene—eny er . paces 


Here 
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. Wises (Z-4) 
Cie * etd 
O= 72° x J : (21.22) 


is the shading function of the primary wave travelling from edge l, 
and 


G(2, y= zee 
(21.23) 


is the shading function of the primary wave travelling from edge 2. 
These functions show that the primary wave from edge 1 undergoes the 
greatest perturbation when 92> » and the wave from edge 2 under- 
goes the greatest perturbation when p=. 


An important property of Equation (21.21) is that it becomes 


zero when a — that is, the field discontinuity which we had 
earlier at the plane x # 0 is completely eliminated. 


In concluding this section, let us return to Expressions (21.11) 
and (21.19) which lead to discontinuities of the fringing field in the 
plane of the strip (x = 0). One may show that the sum cf these 
expressions 


WV ka cos (7-4) 


” u = ; ika fein o—sins 
Hi=Hat vag] a0 \ edge 


i(*+7) 
+ Herlag (I) eter Pg (ape anne 


(21..24) 
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agrees, when Vkacos(+-+-$)>1 , with the asymptotic solution obtained 

in the book [50] by means of integral equations. The solution found 

in [50] has the greatest precision when a=0,"9=0 , and it is completely 
useless 1f e=xt-p or patsy 


§ 22. Secondary Diffraction by a Strip (E-Polarization) 
It is known that a current filament with an electric moment Po 
which is found at a distance h from an infinite, ideally conducting 
Plane (see Figure 46) creates in space the field 


E, = ik*p,e [Hi '(kR,) — HO (kR,)). (22.01) 


With small values of h (and Ry 9? kh?), this expression is trans- 
’ 
formed to the form 


i (amet) 
£.= ~itpwasing Woe), (22.02) 
The primary edge wave is determined by the relationship 
sna) 
E.=Eul(9)l (ae (22.03) 


where E42 (a) 1s the value of the incident plane wave field at the 

point q (Ry = 0). Consequently, the primary edge wave in the direction 
vw * 0 may be investigated as the wave from a current filament located 
above an ideally conducting plane if one assumes the filament moment 

to be equal to 


i f(¢) 
Pe=Fayig Eu l9) sing . (22.04) 


The fie'd, created at the point P by the current filament with 
a moment Po which is parallel to the half-plane's edge and passes 
through the point Q, is determined by Expression (20.11). Expanding 
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the right-hand member of this expression into a series in terms of the 
small quantit, h(h + 0) and limiting ourselves to the first term which 
is different from zero, we obtain 


E.= ik'p. lald, ¥—¥)— 


Fe t(ee-= 
Fe ), (22.05) 


—aartenls V 
a 
By means of relatior.ships (22.04) and (22.05), one may show that 
current filaments with moments Piz and Pie which are located on the 
ideally conducting half-plane —a@<yeoo and correspond to edge 1 (see 
Pigure 47) create in the region lei<> the field 


Wien (5-3) 
Buen Eat cos¢ j ef dg - 


+ asia (tft Fc teeinem sey, 


2 ‘ wes 
+ Eu (1) cos 9 * VEE genes ee (22.06) 


The current filaments with the moments P5, and a which are located 
above the iceally conducting nalf-plane —eceyqa and correspond to 
edge 2 create in the same region the field 


Vive (4-5) 
Ear Eu is cos@ j efde+- 


= sip (+ + Henn wT (2) ema tets e—s18 9) is 


+75 ha 


ifere 7) 
—ihe (sia oie ¢) 


i ig e T 


(22.07) 
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The first terms is. Expressions (22.06) and (22.07) ere the desired 
secondary waves, and the last *erms in the expressicns are the fields 
from the current filaments located abcve the ideally conductinz piane 
Zz © 0 and having the moments 


Po oe aang Eu Tite’ Saute 5 va gles Bua Haye ine oon (5.8 
where 
H+ cece. 3 Ji2v~2 Lo wgee (22.09) 


Summing the secondary waves which have veen found with the 
unperturbed primary waves, we obtain the total field scattered by the 
strip 


1VGen (5-$) 

Fysn =H, = Bar vin{ie i[eose x j efdet 

tg (t-4) 
2Vitew (> +4) 

+1) [cosyx J eag + 


pe, sin 4+ He eubett- sing) nee 


Ret! > sin ‘| ¢ inet ome 


2y/ he 
fe Fea Uf (1) et 8M 4 7 (2) emtawin omar ot }x 
s(97) 
e 
xe = (22.10) 


Now assuming, as in the case of the si-polarization, that the 
moments Pie and Po. depend on the angle ¢, by replacing 


ry Tu) Le, 
1) by soe and jt2) by a (22.11) 
we obtain 
B,m— Hm Eu lPU, 9)f (1) Oe Oy 
e4) 


—jdapeiae—singh @ 
HFG. HIGHS u (22.12) 
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where 


aie wag cos (+-+) 
{ eda 


HMhaitsin gy — 
$5 “aay 
4V ia 


cos( 7-4) 


Pay pod Vine ce (4+ 5) 
eae i J e'“dg + 


etitatl =—sing) + 


WEE +H) | 


F(2, 9)= 


(22.13) 
are the shading functions. They show that the primary wave from edge 
1 undergoes the greatest perturbation near I=—> » and the wave from 
edge 2 undergoes the greatest perturbation in the vicinity of => 


§ 23. The Scattering Characteristics of a 
ale Plane Wave by a Strip 


Expressions (21.21) and (22.12) which were obtained above for 
the field scattered by a strip approximately take into account the 
interaction of the .dges and are valid when |p|<-3 . However, they 
are not applicable with a glancing incidence of a plane wave on a 
strip (when a= ). 


In order to find equations which are applicable in this case, 
let us proceed in the following way. Let us write the expressions 
for the field radiated by the strip in the direction a with the inci- 
dence of a plane wave in the direction » (Figure 45) 

Ey = — H,= Ea F(R, 2) j (In + 
i (™e +) 
+F(I, a) f (2) e7i8s Siam 9] ee, 
HE = Has {G (2. a)g (1) eftaineetneh 
i fare <) 


F e 
+6 (1. a\g (2) enimies a—sin s) ] 7 ra 


. 
Se Seem 2 a 


(234-05) 
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Here lel< 5 » but cannot approximate 21/2. Now let us note that 
the ses ree sions for the fringing field must satisfy the reciprocity 
principle — that is, they must not change with the simultaneous re- 
Placement of a by $ and ¢ by a. Comparing Equations (21.21), (20.12) 
and (23.01), it 1s not difficult to obtain the expressions 


Ey = —'H) = Ed F (1, 9) F (2, a) fryers 
ne) 
HF (2, gyF (1, apf (Qype rerin® eas ; 
Mie E, = Hy[Q(1, 9)G (2, nspcyememenoe 
ieez) 


View” 


+ G2, 9G (1, a) g (2) enim 9 
(23.02) 


which satisfy the reciprocity principle, have no discontinuities any- 
where, and are suitable for making calculations with any values of a 
and e(es. lei<z) - From the second equation of (23.02), it 
follows that Hq, = Ey * Q when e=+> — that is, the fringing field 
does not experience discontinuities in the plane x = 0. Moreover, 
H, * Ey * 0 with any values of 9 if a=+- + — that is, a plane wave 


polarized perpendicularly to the strip does not undergo diffraction 
with a glancing incidence. 


The resulting Equations (23.02) may be investigated as interpo- 
lation equations. Actually, with Jal€3 when Vkacos(7.= 5 )>1 the 
functions Fil. 2), F(2, 2) G(l,a) and G (2, a) are close to one, and 
Equations (23.02) change into the previous Expressions (21.21) and 
(22.12). But if {!el[€4 and Viacos(--~-t)>1, then the functions 
Fil, =) F(2,9) , G (1, 6) and G (2, ¢) are close to one, and Equation: 
(23.02) change into Equations (23.01). Let us recall that the fun:- 


tions F and G are determined by relationships (21.22), (21.23) and 
(22.13). 


In the direction of the principal maximum of the scattering 
diagram (¢ = a), Equations (23.02) take the following form: 
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“, = H,. [ (2ite COS 2 +. NG ],°2) G (2, at 

at 7) 
Vine ' 

F.= Ey] (2k cox 2. ‘gan oO a)F(2, a)-+- 


+G(1, y PEO _ G19, 4) Ge] 


; e 
ane —- 
‘(are +) 


FL 9) SERED (2, ay EM | se (23.03) 


a 
Hence when 3=—>7F we have 


H,=0, 


3 Re Wha 
ema) 6 YR Porat et 
\ A . 


4V ies 
. wie i 
oe eee 
vee tiv § £4 Vye. ear (23.04) 


It is interesting to observe that Expressions (23.02) to some 
extent take into account, in addition to Secondary diffraction, also 
tertiary diffraction. Actually, for the values Isle and lle F 
we have 


G(1, 9) G2 a) eftaisin a—sin g) eltaisine—sing) 


‘ s i s 
efter2+ 10 afsin¢) 7 efta(2—sin a~sin oy T 4 
_————— 9 _ 
2 xka oon (7-4) 2 xka cos (+4) 
i eftait—sin atsin 9) 


eos (F—F)ea StS) (23.05) 


— 


The physical meaning of the four terms in the right-hand mexber of 
this equation is illustrated in Figure 48 (Figure 48a corresponds to 
the first term; Figure 48b corresponds to the second term, etc.). 


Taking into account condition (6.15), one may write the equations 


for the fringing field in the lert half-space ( $<l#<= put jal< ) 
in the same form as (23.02). Thus, the functions G(1, a), G(2; a) and 
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‘a? by a? 


Figure 48. The schematic diagram of the 
waves corresponding to the various terms 
in Equation (22.05). 


F(l, a), F(2, a) will, as before, be described by the relationships 
(21.22), (21.23) and (22.13). The remaining functions in Equat‘ons 
(22.02) will be determined when > * <|9l<e by the following equations: 


: 22¥ie sin (+ +#) 
Gl, Er e . J e'“ug, 


a2 a 
G(2, 9) = = ee" e"'dq, 
0 (23.06} 
a 220 fa sin| | +}) ’ 
F(1, 9) = Jee J eda 
° Ve _— 
0 
mn { ayceraik 
Via sn(Z+ iy 
- Wha stn ( 777) 
F (2, = ee ‘| ae 
2 ily, a (23.07) 
WW ka Ain (3-5) = 
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Calculations of the scattering characteristics were carried out 
based on the equations derived above. These scattering characteris- 


tics are the functions h(a, ¢) and e(a, 9) determining the fringing 
field by means of the relationships 


i (s+ 7) 


eae 


E,= E,,kae(a, of = 


j : 
H,= H,,kah(a. nf ze (23.13) 


The calculations were performed for the values kKa=Y28 and £g=:V80 .~ 
In Figures 49 ~ 62, the following designations were used: 1) the 
functions h and e correspond to the rigorous theory; 2) the functions 
hy and ey correspond to the field from the uniform part of the current 
(the physical optics approach); 3) the functions h, and e, correspond 
to the field from the uniform and nonuniform parts of the current, 
but without consideration of the interaction of the edges; 4) the 
functions hy and e5 correspond to the fringing field with considera- 
tion of secondary diffraction calculated on the basis of equations 
(23.13), (23.02) and (23.10). Thus, in accordance with § 6, 

Ce CON 1 ee, } 


h,=cos et 2 [ka (sina — #10 9) 


~~ ha(sins — sing) (23.14) 


and 


a slofka(sin 2—sin 9)] | jf03 {ka(sins — sia 5)] 
i ne 


s— ‘ a+ is 
h, sin >t cos 4 (23.15) 
where ales, lel<z , 


The results obtained show that our approximation equations arree 
satisfactorily with the rigorous theory already when tu~-1 28 , 
although in the given case approximately one and one-half wavelengths 
are fitted into the width of the strip. In the direction toward the 
source (9=—-e+2. C<ac3 » and also with glancing irradiation of the 
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The various 


The scattering diagram of a field by a strip as a function of 
the incident angle of a plane wave (a) and the parameter vka. 
curves correspond to various approximations. 


Figure 49, 
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Figure 51. The same as Figure 49 when 
@9 = -7 +a, 


strip (a = =-1/2), when the functions €, and Nps and e, and hy lead to 
qualitatively incorrect results, the functions e5 and hs give, as in 
the remaining cases, fully satisfactory results. Actually, the curve 
In,| coincides almost everywhere with the curve |h| (Figure 49-54) 
within the limits of graphical precision. But the calculated values 
of the function le] differ from the corresponding values of the 
function |e] only by hundredths of a percent (Figure 55 - 62). The 
better agreement with the rigorous theory associated with the E- 
polarization is explained by the weaker interaction of the edges in 
this case. A certain discrepancy of the curves |hs | and |h] in the 
vicinity of the principal scattering maximum is explained by the 
interpolation character of our equations. 


As a consequence of the interpolation character of Equations 
(23.02), the integral scattering diameter obtained from Expressions 
(23.03) when a = 0 does not coincide with the interral diameter found 
by Clemmow [46] in the form of the first terms of an asymptotic ex- 
pansion in inverse powers of “ka. However, our equations,as disttnet 
from the similar equations obtained by othcr authors, allow one to 
calculate the scattering characteristics with any Incident anrles of 
the plane wave. 
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The function h (a, 9) for a strip (a = 0, ka = 


0). 


Figure 52. 


137 


BSR wie 
Seal 


td 


T 


The same as Figure 52 when a = -7. 


Figure 53. 
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Figure 54. The same as Figure 52 when @ ® -7 + a. 


Let us note that the functions e(u, $) and h(a, 9) for Figures 
49 = 52 were calculated on the basis of rigorous series which were 
obtained by the separation of variables in the elliptic coordinate 
system (compare C23) ")), 


4 
§ 24. Secondary Diffraction by a Disk 


Let us refine the approximate solution of the diffraction pro- 
blem for a disk which was found in Chapter II. 


Let an infinitely thin, ideally conducting disk of radius a be 
found in free space. Let us orientate the spherical coordinate system 
in such a way that the normal n to the incident wave front would lie 
in the half-plane 9 = 1/2, and form an angle 1 (0<1<4) with the z 
axis (Figure 63). Let us prescribe the incident plane wave field in 


E=E,e"” eID RECeO TY: H= He" slay +e con a) (24 : 01) 


()pootnote appears on page 162. 
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The same as Figure 55 when a = -n, 


Figure 57. 
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Figure 62. The same as Figure 59 when 4 = -1 + a. 


In accordance with § 12, the fringing field in the plane 
¢ = tn/2 is described (when R >> ka°) by the equations 


E,=— A =™ (02.9) — 10.9140 + 
+112 9-10, M40) 5. 


Hy = Ey 5 {let2. )— (1, D+ 


jak 


$ile(2 Y+g(l. AAG) a (24.02) 


These expressions are valid when <5 and <F - The quantities 
included in them are determined by the relationships: 


cos “t. nie 
iQ. J= — aoa 
cosy sta as 
(2, \=—-— gains * 
a(t, Jo —/2, 8), (2, J=— (1%, | 
(24.03) 
C== ka(sin 8 — sin 8), (24.04) 
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Twith e=-> 


|-1 with 9=—>- (24.05) 
Let us note that here 
FQ) with e=+> 
£(1.)= eee 
f (2) withe=—+> 
jie) with e=> 
[24= «’ 
[F(Q) withe=—> 
g() withe=+> 
g(l, »)= a’ 
18 (2)with P= —- = 
jee with p= 
g (2, %) = 2’ 
[aQ)witha=— > 
(24.06) 


an@ the functions f(1}, f(2), g(1) and g(2) are determined by the 
Equations (12.03) and (12.04). 


When c >> 1, Expressions (24.02) take the form 


ae i( 3a } 
E,=—H,=F2[1@, ae = 


Rr (24.07) 


They show that the fringing field in this region may be investigated 
as the sum of spherical waves from two luminous points on the rim of 
the disk with the polar angle yw = t7/2. The diffraction by a disk of 
each of these waves may be studied as was done in the case of a strip, 
but we shall proceed differently. 
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Figure 63. The cross section of igure 64. 
a disk with the plane yoz; n 
is the normal to the incident 
wave front. 


Excitation of a half- 
plane by an elementary dipole 
which is located at the point Q. 


Let us compare the shading of spherical and cylindrical waves 
by a half-plane. Let an ideally conducting half-plane be found in 
free space, and let there be an elementary dipole at the point Q 
(Pigure 64). Let us find the field in the plane perpendicular to 
the half-plane'’s edge and passing through the point Q. 


In accordance with the reciprocity principle, it is determined 
for the electric dipole by the relationship 


= Pt 
Ee = 5 Ee (Q), (24.08) 
and for the magnetic dipole by the relationship 
ms, 
= 2 HQ. (24.09) 


Here p,(m,) is the electric (magnetic) dipole moment found at the 
point Q; Poe and Mj, are the moments of the auxiliary dipoles which 
are placed at the point P; E,(@) and H,(Q) are the fields created by 


the auxiliary dipoles at the point Q. 
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Now let us remove the auxiliary dipoles to such a distance that 
the spherical wave arriving from them may be considered to be a plane 
wave on the section from the half-plane's edge to the point Q. In 
this case, in accordance with Equation (20.08), the field created by 
the wave at the point Q will equal 


E, (Q) = E.. (0) (uld, o — 7’) — u(d,¢ +-9")], 


Hs (Q) = Hys (0) [u (d, 9° —9”) 4-4 (d, ¢ +-9”)}. (24.10) 
The expressions 
é 
Ex) =k pa G, Hye 0) = ery (24.12) 


determine the fields created by the auxiliary dipoles in free space 
(with the absence of the half-plane) at the point 0. 


Consequently, the fields excited at the point P by the electric 
and magnetic dipoles which are found at the point Q above the half- 
Plane equal respectively 


E,=p.[u(d,9 —9") ald, toytn, } 
Hiab [u (9 —9") Lady + ey. J (24.12) 


With the absence of the half-plane, these dipoles create at the point 
P the field 


oa iar 
— "@ ~lede get 
E,=k' * Re i la di 


H,= k'*m, - ev itd con (6 — pre) (24.13) 


Comparing Expressions (24.12) and (24.13) we find the shading 
functions 


F= fu (d, P’ -- ?”) —2 (d, 9’ + ?”)] etd cos (6° —~<? ZR 
[ (24.14) 


G=[u(d. 9 — 9") 4u(d, 9 + gry ettrorer—w-, 
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In the case wher the current filament passes through the point 
Q parallel to the half-plane's edge, the field at the point P is 
determined — in accordance with (20.11) and (20.12) — by the 
equations 


E,=ik"p,(u(d,¢ — 9") —a(d, 9 + 9")] Viz )] 
—- t{ag—= 
Haidar —o) 40.9 te we ‘| (24.15) 


With the absence of a half-plane, these sources create at the point 
P the field 


(re 


s 
E,=ik'p, Vz e 4 ) en itd coe (997) 


ifaae > 
H,=ik'm, i= e ( 7 ) eid cos (5°97), 


(24.16) 


Comparing Equations (24.15) and (24.16), we obtain the same Expressions 
(24.14) for the shading functions. Consequently, a spherical wave in 
the direction perpendicular to an ideally conducting half-plane is 
shaded by it the same as a cylindrical wave. 


Let us note, however, that Expressions (24.14) are not equivalent 
to Expressions (21.22), (21.23) and (22.13), since the first represent 
the shading function by a half-plane of a wave from a sirgle source, 
and the latter represent the shading function of an edge wave which 
we approximate by waves from two sources located on both sides of the 
corresponding half-plane. Since the shading functions of spherical 
and cylindrical waves are the same, the edge wave shading functicns 
of a strip and a disk also will coincide. 


Therefore, the approximation expressions for a field scattered 
by a disk which take account of secondary diffraction may be repre- 
sented in the region p=---~,0<9<-+ (with ¢ >> 1) in the following 
form: 
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Starting from Expressicns (24.18) and (24.19), it is not diffi- 
cult to write interpolation equations for the fringing field which 
are suitable for any values of y and 6 in the interval (0; 772), 
but when 9 = 2/2; 


jaE,, 


Ey — Hy {LF (2,8) F(1,2) (2,8) — 
—F (1.9) F (2,4) £(142)]44 ©) ELF (2.9) F (1, 2)/ (2, 894 
FF (19 F (2,2) (1,91 4,0)} > (24.20) 


£,=H, =o {Ig (2, 0) G(1,2)g¢2, a= 
— G (1,8) G (2,4) a1, 2) J, (©) + 1G (2,0) G1, 2) g(2, 8) + 
£6(1,9)6(2,2) 2(1, 21,0) (24.21) 


Let us note that when y = 0 these expressions will be valid for 
any values of the azimuth $¢ , since then any point of space may be 
considered to be located in the incident plane. \ 


In the direction of the scattering diagram's principal maximum 
— that is, when 9=7,.9?=+ — the fringing field (24.20) and (24.21) 
takes the form 


2 jan 
E,=—H, ait By -F nF (bn cost, | 


iar 
=H, = BE HG (2,16 (1.1) ZF 008 1- (24.22) 


However, these expressions have an interpolation character,and with 
small values of the angle y it is impossible to consider them to be 
more precise than the simple equations of § 9 and § 12. In particular, 
with y = 0, when the fringing field must not depend on the incident 
wave polarization, they give values which are different for the E- 
polarization and H-polarization by small yguantities of the order of 
= . Therefore, in this case (when y = 0) it makes sense to use 
Expressions (24.20) and (24.21) only far from the z axis, switching 

to Equations (24.02) near the z axis. 
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Equations (24.20) and (24.21) have the following important pro- 
perties, They do not have discontinuities, they include the case of 
glancing incidence of a plane wave, and they satisfy the reciprocity 
principle. From them it follows that £,=H,=0 when §8=+ — that 
is, the fringing field does not experience a discontinuity on the 
plane z = 0. Moreover, E,=H,=0 with any values of ®& if Y=-p 
— that is, a plane wave polarized perpendicular to the disk's plane 
does not experience diffraction with glancing irradiation of the 


disk. 


As in the case of diffraction by a strip, the new approximation 
expressions consider to some extent tertiary diffraction [see 
Equations (23.05) and Figure 48]. 


i 


Using Condition (9.04), it is not difficult to write equations 
*£ 
for the fringing field in the left half-space (7 <O<5 e=t +) 


Ey — Hy = LF (2. #0) F (1,8) (2,2) — 
— F(1,a— 0) F (2,3) 7 (1,2) 4, (+ 
PALF (2,—8) F (1,3) (2,8)-+- 


Z eit? 
HF (1e— F(Z (1L NAO. - (24,23) 
Ey=H,= <P {6@.2—HG0, 2) g (2,8) - 
—G(1,2— )G (2,8) 2(1, 2) 4,4 
Fi 1G (2,2 —8)G(1, 2) g(2,8)+ 
$G(1,2—2)G (2,3) 9(1.2)4.0 } 2, (24.23) 


where the functions f and g are determined by the equations 


. ‘ } 
g (1,8) =f (2, 8) == ——_ ————-—_ 


sind—sinte i 
cos ai —~sin os i 
f(1.8)=g(2,=——*2,___?_ | 
aoe : (24,24) 
In the direction towards the source (=*—¥, 9=—-5) , the 


fringing field equals 
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aE, 
E,=—-U,= = (LF? (1.8) f (2.8) — 


— FP(2,8) f (1.2915, @)-E AEE, 8) (258) 
+ F271 OL 

Ey=H,= 9” fic*(1, 2) g(2,8)— 

— G* (2,8) 9 (1,291, +4 (64(1.7) 902.8) 4 


+6" (2,2) 2(1.3)14, 0} 5 
DaLMAOpy-. | (24.25) 


where 


2siny 


_ s-t 
(1,8) = 9 (2,8) = —L tear, | 


Lo _ t—sinyg 
12.8)= 801.8 ag (24.26) 


As was already noted, it makes sense to use Equations (24.25) 
only far from the z axis, changing to Expression (12.15) of the pre- 
vious approximation in the vicinity of the z axis. A calculation of 
functions E(e—y) and S(s—y) (Figures 65 and 66) which determine 
the effective scattering surface [see Expressions (12.17)] was per- 
formed on the basis of these equations when ka = 5. A comparison was 
carried out of this calculation with the results of measurements. 

The two experimental diagrams (the dashed Lines) (2) depicted in 
Figure 65 characterize the experimental precision. As distinct from 
the previous approximations, which lead in this case to qualitatively 
incorrect results [see Equations (10.06), (10.07) and (12.15)], we 
observe a satisfactory agreement of theory with experiment. 


For verifying the results obtained, a calculation was also 
carried ou* of the functions video ) and vi2Ic4 ) [see Equations 
(9.07)] when ka = 5 (Figure 67 and 68) with normal irradiation of a 
disk by a plane wave. Curve 1 corresponds to the field calculated 
from the rigorous theory [34]; curve 2 corresponds to the field from 


(2) 


Footnote appears on page 162. 
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Figure 65. The diagram of a disk's effective 
scattering surface when the plane wave's 
magnetic vector is perpendicular to the inci- 
dent plane. 
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Figure 66. The calculated diagram of a 
disk's effective scattering surface when 
the plane wave's electric vector is per- 
pendicular to the incident plane. 


the uniform part of the current (the physical optics approach). Curve 
3 corresponds to the field from the uniform and nonuniform parts of 
the current, but without the interaction of the edges. Curve 4 corre- 
sponds to the field with consideration of secondary diffraction. 43 
is seen from these graphs, consideration of the edge interaction re- 
fines the previous approximation and ensures better agreement with 

the rigorous theory results. 


The problem of secondary diffraction by a cylinder may be solved 
by a similar method. However, considering that the corrections which 
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(1) 
Figure 67. The function Vig for a 
disk with normal incidence of a plane 
wave (curve 4). Curves 1, 2 and 3 from 
Figure 20 are drawn for comparison. 


depend on the secondary diffraction here are small (on the order of 
1 dB) when ka = m, kl = 10m, and the equations are substantially more 
complicated, we shall not cite them here. 


In the problems investigated above, the edge waves have the 
character of cylindrical or spherical waves — that is, they decrease 
rather rapidly with the distance from the edge. Therefore, in the 
ease when the linear dimensions of the faces are approximately two 
wavelengths, it is sufficient to limit ourselves to a consideration 
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Figure 68. The function viz) for a 
disk with the normal incidence of a 
plane wave (curve 4). Curves 1, 2 and 
3 from Figure 21 are drawn for 
comparison 


of only secondary waves. In Chapter VII we will investigate the 
problem of a dipole in which the edge waves decrease so slowly that 
it is necessary to consider multiple diffraction. 


§ 25. A Brief Review of the Literature 


In this and previous chapters, approximation expressions were 
obtained for the scattering characteristics of a plane wave by various 
bodies. These expressions were derived with the help of physical 
considerations which do not pretend to be mathematically rigorous, 
and they are adequate for sufficiently short waves. In the literature, 
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there are a number of works in which similar results were obtained. 

A majority of these works also are not characterized by mathematical 
rigor, and they are based on certain physical assumptions. Therefore, 
one may relate them to the physical theory of diffraction. Only in 

a few works (related to the simpler diffraction problems) did they 
succeed in obtaining specific results at a higher level of mathematical 
rigor ~~ more precisely, while developing asymptotic methods of mathe- 
matical diffraction theory. 


We will briefly list the most important results obtained in a 
number of papers and books, grouping the material in the following 
sequence: 


l. Diffraction by plane, infinitely thin plates (an infinite 
strip, a circular disk) and diffraction by auxiliary apertures in a 
flat screen (an infinite slit, a circular hole). 


2. Diffraction by three-dimensional bodies with edges (a finite 
cylinder, a finite cone, etc.). 


3. Other diffraction problems. 


When investigating the first group of diffraction problems, it 
is necessary to keep in mind the principle of duality [4] which enables 
one to easily change from a strip to a slit, from a disk to a circular 
hole, etc. In the literature as a rule, they preferred to investigate 
apertures in an infinite flat screen, whereas in our book, diffraction 
by a strip and a disk was studied. This approach facilitates the 
transition to three-dimensional bodies (see the remarks at the 
beginning of this chapter). 


Based on the time of appearance (if we do not consider the works 
of Schwarzschild [25] which we talked about in the Introduction), one 
should first of all mention the works of Braunbek [28 = 30] which 
were devoted to the diffraction of a scalar wave by a circular hole 
in a flat screen. Assuming that the plane wave is incident normal to 
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the screen, the author obtained an approximation solution in the form 
of a surface integral. The boundary values of the integrand were 
taken from the rigorous solution to the problem of diffraction by a 
half-plane which was found by Sommerfeld. The field was calculated 
in the far zone on the axis of the hole and far from it, and also on 
the axis near the screen. Using this approach, Braunbek recently 
solved the problem of scalar wave diffraction by an aperture in a 
concially shaped screen [31]. 


In the papers of Frahn [32, 33], this method was used for the 
diffraction of electromagnetic waves. Diffraction of a plane wave 
incident normal to an ideally conducting screen with a circular hole 
was investigated. The field was calculated in the hole and on the 
axis, and also the field in the far zone and the transmission coeffi- 
cient (the ratio of the energy passing through the hole to the energy 
falling on it) were calculated. 


In these works of Braunbek and Frahn, secondary diffraction was 
not considered. The expressions obtained by them for the fringing 
field intensity in the far zone agree with similar expressions 
following from our equations (§ 9). 


Karp and Russek [51] studied diffraction by a slit in the case 
when the incident wave's electric vector is parallel to the slit edge. 
They investigated each semi-infinite p.urt of the screen as a half- 
plane excited by the incident wave field and a "virtual" source 
localized on the edge of the opposite half-plane. The moments of 
these sources were determined from a system of two algebraic equations 
which were obtained by using the asymptotic expressions resulting 
from the rigorous solution for the half-plane. Secondary diffraction 
was considered, and partially the general interaction. Special 
attention was allotted to calculating the transmission coefficient, 
but equations for the scattering characteristics which would be 
suitable with all directions of incident wave propagation were absent. 
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Clemmow [46] and Millar [47 ~- 49] in their works calculated the 
transmission coefficients with normal irradiation of a slit and a 
hole, and also the field in the hole. The solution was sought by 
means of curvilinear integrals of the fictitious linear currents on 
the aperture edges. The interaction of the edges was considered. 

The case of inclined irradiation was not investigated, since it turned 
out to be too complicated for investigation by this method. 


The "geometric theory of diffraction" of Keller [42 - 44] which 
deals with diffraction rays is of special interest. The phase and 
amplitutude corresponding to each diffraction ray are determined at 
each ray point on the basis of geometric considerations and the law 
of the conservation of energy. The initial diffraction ray amplitude 
is assumed to be proportional to the incident ray amplitude at the 
point of its diffraction. The unknown proportionality constant be- 
tween the amplitudes and the initial phase difference -s determined 
from a comparison with the results of well-known solutions of diffrac- 
tion problems. In this way, the fields scattered with the normal 
incidence of a plane wave onaslit and hole in a flat screen are 
found. These fields are obtained with consideration of multiple dif- 
fractions, but they are not precise wave equation solutions, since 
their calculation was started from approximation relations. Moreover, 
geometric diffraction theory is not applicable near caustics, and also 
in the vicinity of the scattering diagram principal maximum. 


In a recently published paper of Buchal and Keller [52], a new 
method for the solution of diffraction problems for holes in a flat 
screen was proposed. The caustics and shadow boundaries here are in- 
vestigated as thin boundary layers, inside of which a rapid field 
change takes place. This method supplements geometric diffraction 
theory, and in particular enables one to find the field at caustics 
and on the shadow boundary. 


Recently, the method of integral equations has been applied to 


the solution of diffraction problems of holes in a flat screcn. In 
particular, Greenberg [53, 54] reduced the solution of this problem 
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to an integral equation for a "shadow" current which 1s, in our termi- 
nology, half the nonuniform part of the current. The resulting 
integral equations may be solved (with any ratio between the dimensions 
of the hole and the wavelength) by the method of successive approxi- 
mations. Moreover, they allow one to obtain asymptotic expressions 
which are suitable for short waves. In Reference [55] Greenberg found 
an asymptotic expression for the current on a strip with ka >> 1 (2a 
is the strip's width). Greenberg and Pimenov [56] obtained a similar 
solution in the case of normal incidence of a plane wave on a circular 
hole. Using the same method, an asymptotic expression was found for 
the current on a flat ring [57], the width and inner diameter of which 
are a great deal larger than the wavelength. 


The above listed works [53 - 57] already relate to the mathemati- 
cal theory of diffraction: in them the first terms of the asymptotic 
expansions for the current were obtained with the desire evidently to 
also be able to calculate the following terms. Unfortunately, the 
asymptotic expressions which have been found up to now refer only to 
currents, and one is obliged to calculate the scattering characteris- 
tics by means of numerical quadratures [56]. As a consequence of the 
rapid oscillation of the integrands, such a method leads to rather 
unwieldy calculations and does not enable one to formulate a clear 
representation of the fringing field formation, and also does not 
allow one to study this field properly. 


Millar [58] investigated the problems of electromagnetic wave 
diffraction by slits in a flat screen. The system of integral equa- 
tions obtained by him for the current is solved by the method of 
successive approximations. The field in the hole is calculated from 
the currents which are found, and then on the basis of the field in 
the hole the field in the far zone and the transmission coefficient 
are calculated. All the indicated quantities are represented in the 
form of an asymptotic expansion in reciprocal powers of the parameter 
“ka. A solution also is obtained in the case of glancing incidence 
of a plane wave. 
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Let us note that the asymptotic expressions obtained by the 
method of integral equations are distinguished by their considerable 
complexity, and frequently require tabuiation of the new special 
functions appearing in the expressions. 


In the recently issued volume of Handbuch der Physik [50], which 
is devoted to diffraction theory, the complex characteristic of plane 
wave scattering by a strip was studied directly, omitting the ealcu- 
lation of the currents. For this characteristic, a singular integral 
equation was formulated, the solution of which was sought in the form 
of an asymptotic series in reciprocal powers of vka. The first term 
of the series corresponds to Equations (6.14) and (6.16). The follow- 
ing term takes into account the interaction of the edges, and becomes 
infinite with the glancing incidence of a plane wave and also for 
observation points lying in the strip’s plane. Therefore, the simple 
expressions obtained in [50] do not allow one to construct the com- 
plete scattering characteristic. In [50] diffraction by a disk, a 
sphere, and an infinite circular cylinder was investigated, and also 
a review of the general methods of diffraction theory and a biblio- 
graphy encompassing a large number of works (mainly German and 
American) were given. 


The book of King and Yu[592 presented (as a rule without deriva- 
tion) a series of asymptotic expressions relating to a slit and a 
circular hole and also to other diffraction objects. Here, however, 
equations from which one would be able to construct the scattering 
characteristics of a strip and a disk with any incidence of a plane 
wave also are missing. 


Works on diffraction by three-dimensional bodies having edres 
are comparatively scarce. In the paper of Siegel et al. [41], the 
effective scattering surface for a finite cone with the incidence of 
a plane wave on it along the symmetry axis is calculated from clemen- 
tary arguments. The expressions obtained here do not fully character- 
‘fze the fringing field, and are suitable only for sharp cones to which 
we already referred in § 17. In the papers of Keller [44], the 
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diffraction ray concept is used for calculating the scattering of 
scalar and electromagnetic plane waves by a finite circular cone with 
a flat base and also by a cone having a spherical rounding off instead 
of a flat base. The resulting expressions :are not applicable in the 
vicinity of certain irradiation and observation directions. In § 17 
we showed that the field scattered by a cone and by certain bodies of 
rotation is not expressed only in terms of the functions f and g, 
which refer to diffraction rays diverging from a wedge edge. This 
result evidently attests to the impossibility. of complete calculation 
of the scattering characteristic with the diffraction ray concept. 


Diffraction problems arising in antenna theory are usually dis- 
tinguished by their great complexity, since the corresponding metal 
bodies (mirror, horn, etc.) have a compiicated shape. Since the dimen- 
Sions of these bodies and the dimensions of the radiating apertures 
are considerably larger than the wavelength, the application of 
physical diffraction theory to antenna problems is very promising. 
only the first steps have been taken in this direction. Thus, Kinber 
[60, 61] performed a calculation of the decoupling and lateral radia- 
tion of mirror antennas. The feature specific to mirror antennas is 
that diffraction rays arising at the mirror's edge undergo multiple 
reflection on its concave surface. This multiple reflection was 
studied by Kinber in more detail as applied to the concave surface of 
a cylinder and sphere [62, 63]. 


Diffraction problems relating to an antenna dipole —a thin 
cylindrical conductor — are investigated in Chapter VII, and 
references to the literature are also given there. 


In conclusion, let us say a few words about diffraction of short 
waves by smooth bodies. The basic principles relating to such pro- 
blems were set forth in the fundamental works of Fok and Leontovich. 
These principles were established by the following methods of 
mathematical diffraction theory: 
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FOOTNOTES 


Footnote (1) on page 138. These calculations were performed 


under the guidance of P. §. Mikazan. 


Footnote (2) on page 150. See the footnote on Page 86. 
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CHAPTER VI 


CERTAIN PHENOMENA CONNECTED WITH THE NONUNIFORM 
PART OF THE SURFACE CURRENT 


In the previous chapters, a theoretical investigation was 
conducted of the field radiated by the nonuniform part of the current. 
In this chapter we will discuss a method for measuring this field 
(§ 26) and we will investigate the phenomenon of the reflected 
signal's depolarization (§ 27). 


An experimental method for measuring the field from the nonuni- 
form part of the current was first proposed for bodies of rotation in 
the paper of Ye. N. Mayzel's and the author [12]. Later it was shown 
that this method has a universal character, and is suitable for 
measuring the field from the nonuniform part of the current excited by 
a plane wave on any metal body [13]. 


§ 26. Measurement of the Field Radiated by the 


Nonuniform Part of the Current 


Let an ideally conducting body cof arbitrary shape be found in 
free space. A surface element of this body is shown in BPirure 69, 


The coordinate system was selected in such a way that its origin would 
lie near the body, and the source Q would be located in the plane 

x #0. If the distance between the body and the source is a great 
deal larger than the body's dimensions, then the incident wave in the 
vicinity of the body may be investigated as a plane wav-. iet us 
reorssent it in the form 


E,= EL, eter rescon p60, (26.01) 


Here y is the angle between the normal N to the wave front and the 
2 axis. 


Now let us place in front of the source, parallel to the radiated 
wave front, a polarizer P which transformed linear polarized radiation 
into a circularly polarized wave. Let the wave passing through the 
polarizer with an electric vector En lag in phase by 90° behind the 
wave with an electric vector E, (Figure 70). In this case, the pcolar- 
izer achieves a clockwise rotation i » AS a result, the incident 
wave field at the coordinate origin will equal 


i 
i (26.07) 


The field scattered by the body may be represented in the wave 
zone .n the following way: 


° late et OO ltt 
p= Et Se Oe 
E, Ss 
iaE,, & e 
E,=4,=—9 yz Yon y. (26.03) 


where a is a certain length characterizing the body's size and £(4,9) 
anc (8,9) are unknown angular functions. In the general case, the 


fry 


wee 


Poctnote appears on page !74. 


field (26.03) is an elliptically 
polarized wave. In the direction 
toward the source (0=*—1, 9=—+) ; 
this wave passes through the 
polarizer and creates behind it 

the field 


Figure 69. The problem of elec- 
tromagentic wave diffraction by (26.04) 
an arbitrary metal body. 


aS - is a surface element of the 


body where 
N - the normal to the incident 
wave front, 
Q - the source, v =t/y Vf). 
P - the polarizer converting m=, 2 a (26.05) 
linearly polarized radia- . 
tion to a wave with 
circular polarization. If the source radiates a 
wave of another polarization 
(Hy | yoz), then the wave reflected 
a. by the body and passing through 
@, ee the polarizer is described at the 
fro ; point Q by similar relationships 
7 ink i= 
iaHes é T 
7 H,=£,= 2 Re oo ' 
Figure 70. inn 
jaH e t 
H,=—E,= “ia Vo. } 


(26.06) 


Now let us investigate in the physical optics approach the 
diffraction of a plane linearly polarized wave by the same body. 
According to definition (3.01), the uniform part of the current exci- 
ted on the body's surface by a plane wave with E-polarization of the 
incident wave (E, | yoz) equals 
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f=- a E.2(n, siny +n, cos pe”, | 


R= se Eys-me-siny ec, 


f= ge Euan, cos ye’, | (26.07) 
and with H-polarization (Hy | y0z) 
fi =0, | 
=r Mus nye”, 
¢ 
f=-M n, e | (26.08) 


Here Foy and Hox are the electric and magnetic field amplitudes of 

the incident wave with E-polarization and H-polarization, respectively; 
yokiv'siny--z’cot) 18 the incident wave phase at the point (x’, y’, 2°) 
on the body's surface; ny» ny» n,» are the components of the normal 

to the surface at the same point. 


Furthermore, calculating the vector potential in the far zone 
on the basis of this current and substituting its values into the 


equations ° 


od 


E, =H, = ikA,, (24.09) 
we find the fringing field. With E-polarization, it equals 
A iaR 
E,=—H,= - Exe ftn, sin { cos @-+ 

L(n,siny-+ cos yisin 9] e*ds, (26.29) 

; eink : , 5 . 
E,=H4,= 35 Eus'-y-" {fa (sin {cus sin p — cos7sin 4) — ‘ 
7 . . (26.19) 

—(n, sity ‘=n, cos yy eos ¢ cos A) oS, 


and with H-rolarization 


UO De a, Wat I-72 


'= 
> 
“yr 


\ R 
. B= H,,C0s@ ef n,e'*ds, 


. a . 
B= Ha rhe le sin® + n, sing cos 9) aS, (26.11) 


Here R, ® , $ are the spherical coordinates of the observation point, 
®=y—krcossM , and integration is carried out over the 1lluminated 
elements of the body's surface. In the case of radar when the obser- 


vation and irradiation directions coincide (9==—-1, 7=—-¢), 
Equations (26.10) and (26.11) yield 


i ak : 
£,=—H,= -= 20a Jinysin t+ 
+1, c0s y) e* dS, (26.12) 


F£,=H,=0 \ 
and 


ek 
Eye = ie Mee FSi sin + 


++ n, cos y)e* dS, 
£,=H,=0. ; 


(26.13) 


Furthermore, assuming the incident wave amplitudes ar. specified by 
Equation (26.02), let us write Expressions (26.12) and (26.13) in the 
following way: 


% 
4 YY 
Ee=—H, = SB, | 
EH, = 12k ot eae, | 
eae yz 8k" | (26.14) 
where 
B= — St ZF [inysiny + a cos yeas. (26.15) 


Now let us represent the angular functions of fringing field 
(26.03) in the form 
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f-S48. | 


Ya w+, (26.16) 


where the functions aa f° ana rl fi refer to the field radiated by 
the uniform and nonuniform part of the current, respectively. Substi- 
tuting these expressions into Equations (26.04) and (26.06) and taking 
into account relationship (26.15), let us find the fringing field 
passing through the polarizer P toward the source Q. In the case of 
E-polarization, 1t equals 


Rik 


E, =-A, = fakes (gt + 81) cen 


£E,= Ze iat (2° ee Le (26.17) 


and in the case of H-polarization 


(eHow 


; . waeoue. 
H,= E, 2: Mer gt p y) “es ‘ | 
H,=-—£,= (se E'.. By ele | 


(26.18) 


The physical meaning of the result obtained is as follows. The 
field scattered by the body at the point Q is the sum of two waves 
polarized in mutually perpendicular directions. The reflected wave 
which is polarized the same as the primary radiation of the source 1s 
determined by the function = (St +E) » and is created only by the 
nonuniform part of the current. The reflected wave with the perpendi- 
cular polarization is described by the function %.==2%?-]-E'—%', and 
is the field radiated by both parts of the current. Let us note that 
in the general case the functions si and zi do not coincide, and 
therefore they are not balanced out in the expressions for Z_. In 
other words, the field radiated by the uniform part of the current in 
this case may not be separated from the fringing field. 


Thus, the investigated method allows one to separate from the 


total field scattered by any metal body of finite dimensions ‘hat part 
of the field which is caused by a distortion of the surface (the 
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curvature, a sharp bend, a point, a bulge, a hole, etc.). One chould 
note that, in the case of electromagnetic wave scattering by a system 
of separate bodies, the separable part of the field is due not only 
to the surface's distortion, but also to the diffraction interaction 
of the bodies. 


It is necessary, however, to keep in mind that it is possible to 
realize the indicated fringing field distribution not in an arbitrary 
observation direction, but only in a direction for which the condition 
E = aye is fulfilled — for example, in the direction towards the 
source. 


Consideration of the nonuniform part of the current also enables 
one to explain the reflected wave depolarization which we will inves- 
tigate in the following section. 

(2) 


Figure 71 presents the results of measurements and calculations 


of the effective scattering surface 


at = atl, [P= eg EPP ENS, (26.19) 


which is dependent upon the nonuniform part of the current excited by 

a plane electromagnetic wave on a disk. The disk's diameter equals 

2a = 5A (A is the wavelength). The calculations were performed with 
consideration of the secondary diffraction on the basis of the approxi- 
mation equations for the functions £ and £ which were derived in § 24. 
Since it is difficult to prepare a thin disk with a sufficiently flat 
surface, the measurements were performed with an obtuse cone close to 
the shape of a disk and having a height approximately equal to one 


tenth of the d.ianeter. 


As is seen from Figure 71, the theoretical and experimental curves 
are fairly close together. A certain divergence between them, 
especially in the region of y values close to 90°, may evidently be 
explained both by the model's conical shape and also by the 


(2) pootnote appears on page 174. 
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Figure 71. Diagram of the radiation from 
the nonuniform part of the current flowing 
on a disk. 


approximation character of the computational equations. The value 
y = 90° corresponds to the direction along the disk's surface, and 
the value y = 0° — to the direction normal to the disk. 


§ 27. Reflected Wave Depolarization 


Let us again return to the problem of scatterinr of an electro- 
magnetic wave by an arbitrary metal body. The relative position of 
the source Q, of a surface element of the irradiated body, and of the 
coordinate system 1s shown in Figure 69. Let us recall that the 
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source Q is in the plane y0z, and 

radiates a linearly polarized wave. 

Furthermore, we shall assume that 
the polarizer P which is shown in 


* 
Figure 69 is now absent. 
& Let us designate by a the 
angle between the plane y0Oz and 
Figure 72. the incident wave electric vector 


Ey (Figure 72). The field of 
this wave will be represented in 
the form 


é, = H,, = E,,e""" slay es ere a 


(27.01) 


H, a an E,, = 4,0" a+ztos 9) 
where 


; . E.. 
E,, = E,sin2. Hua =— E, cosa, 5 = tgs. (27.02) 
The field scattered by the body is determined in the wave zone 
by the equations 


E,=— Hy =F Eh (1 §,9)+ 
ett 


+ HysS, (y, 3, 9)) R’ 


E,=H,= B Ets (t. ®, 9+ 
an 
+H Blt 8.9) Ge 
(27.03) 


Here a is a certain length characterizing the body's dimensions, R, 
® , ¢ are the spherical coordinates of the observation point, 
E,.(7.% 9) and £,,(y,8 9) are unknown angular functions. 


It is obvious that the fringing field polarization — that Is, 


the orientation of its electric vector in space — depends in a com- 
plex way on the observation and irradiation directions. In the 
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direction toward the source, it may not coincide with the polarization 
of the wave radiated by the source. Such a phenomenon is called 
reflected wave depolarization. 


It is easy to establish the reason for depolarization, if one 
investigates tne fringing field as the sum of the fields radiated by 
the uniform and nonuniform parts of the current. According to § 26, 
the uniform part of the current radiates the following field in the 
direction towards the source (92-1, =—+) 


£,=—H, = “G2 5-¥, 


‘\ 


iaE,, v** = 


icH,, #8? 
E,=H=—* a (27.04) 
The functions a and 29 satisfy the condition 29 = ae’. and are de- 


scribed by Equation (26.15). From Equation (27.04), let us immediately 
obtain the equality 


fe 
Fy = 8% (27.05) 


which means that in the physical optics approach the reflected wave 
does not experience depolarization. Consequently, the reflected wave 


depolarization is caused only by the nonuniform part of the current 
or, in other words, by the surface distortion. 


Let us derive an equation for the magnitude of angle 6. This is 
the angle by which the electric field vector of the reflected wave is 
turned in respect to the electric vector of the wave radiated by the 


source. For this purpose, let us represent the functions T1(2) and 
E402) in the form 


Fo oe 
B= Ey + Sy | 


=) =Sy a’ | (27.06) 


where the terms Dee and aes correspond to the field radiated by 


the uniform part of the current, and the terms een and ere 
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correspond to the field radiated by the nonuniform part of the current. 
Comparing Expressions (27.04) and (27.03), we find that 


y= 2°, z =0, ) 


= =0, re | (27.07) 


Therefore, the field scattered in the direction towards the 
source (®=«-1. =-+). will equal 


E,=—H,y= BiB E+5)+H,.3) SS, 


Ey= Hy = “9 [Ey8) — Hee (EP — SS. | (27.98) 


This field's electric vector forms an angle 8 with the yoz plane. 
The angle 8 is determined by the equation 


_ Ee _ S48 —Hctge 
I= Podetige F* (27.09) 


As a result, the desired angle 6 which characterizes the depolarization 
magnitude will equal 


b=a—p. (27.10) 


Thus, the field from the nonuniform part of the current, separatle 
"in a pure form" by means of a polarizer (§ 26), leads to depolarizatior 
of the scattered radiation. 


Specific results from the depolarization calculation of waves 
reflected from certain bodies may be found, for example, in the works 
of Chytil [75 - 77] and Beckmann [78]. In particular, in Reference [77 
it was shown that the depolarization effect on the effective scatter- 
ing surface of convex bodies in practice may be neglected only with 


the condition ka 2 4, 
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Footnote (1) on page 164. 


Footnote (2) on page 169. 
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FOOTNOTES 


A system cf metal rates parallel 
to the e, vector may serve as the 


simplest example of such a 
polarizer. 


See the footnote on page 86. 


174 


CHAPTER VII 


DIFFRACTION BY A THIN CYLINDRICAL CONDUCTOR 


In almost all the works devoted to the diffraction of plane 
electromagnetic waves by a thin cylindrical conductor, the current in- 
duced in the conductor was studied, and then, by intecrating this 
current, the fringing field in the far zone was calculated. However, 
in view of the complexity of this problem, they succeeded in obtain- 
ing relatively simple equations only in the particular case when the 
observation direction and the direction toward the source coincided, 
and was perpendicular to the conductor axis. In the general case 
when these directions did not coincide and were arbitrary, the expres- 
sions for the fringing field became very complicated and unsuitable 
for making calculations. Since they were obtained by integrating 
approximation expressions for the current, it turns out that they have 
still one other shortcoming — they do not satisfy the principle of 
duality. 


In this chapter, explicit expressions are obtained for the 
fringing field which are suitable for making calculations with any 
direction of irradiation and observation. We shall consider both the 
primary edge waves excited by the incident plane wave and also the 
secondary, tertiary, etc. edge waves. The total fringing field ie 
found by summing all the diffraction waves. 
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§ 28. Current Waves in an Ideally Conducting Vibrator 


The electrodynanic problem of determining the current in thin 
cylindrical conductors (vibrator) usually is reduced to an interro- 
differential equation. The latter is derived by means of boundary 
conditions on the conductor surface, and is substantially simplified 
in the case of thin conductors when the inequalities 


+ <1 and ha <1, (28.01) 


are fulfilled, where a is the radius and L is the length of the 


conductor and kat =— 5 


Its solution may be found, for example, by the method of succes- 
sive approximations [79, 80] or by the perturbation method [85]. 
Recently, Vaynshteyn [81, 82] proposed a new solution for this equation. 
Since we will subsequentl; base our work on the results of References 
(81, 82], let us discuss them in more detail. 


Let us assume that the vibrator's symmetry axis coincides with 
the z axis, and its ends have the coordinates z = zy and z = Zo 
(L = Zz z,)- In the case of excitation of the dipole by a concen- 
trated external field 


E* = Gi(z) (28.02) 


the current J(z) in the conductor may obviously be written in the form 
of the sum of the waves travelling along the conductor with a velocity 
ec from the excitation point z = 0 and the ends z = zy and = = Z5° In 
Reference [81] it was shown that the complex amplitudes of these waves 
are slowly varying functions of the z coordinate. These functions 

may be approximately expressed in terms of the function w(z), so that 
we obtain the following cxpression for the current J(z): 


J(zy=Jfo(zyer tt ay(z—2 jeer 
+A, 4 (2,—2) MO), (28.03) 
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Here the quantity 


A= ce. » Y= 1,781... (28.04) * 
determines the initial value of the current wave propagated from the 
excitation point +), The function w(z) is the solution of the inte- 
gral equation, and in addition to the variable z it also depends on 
the parameters k and a. We will not list here all the properties of 
the function y(z), but let us note only that it satisfies the 
conditions 


4(0)-=1, p(oo)=O, (28.05) 


and its absolute value monotonically decreases with an increase of z. 
This decrease, which is rather slcw and does not have an exponential 
character, is due to radiation. 


The constants Ay and A, determine the initial vatues of the 
current waves originating at the points z = Z4 and z = Zo» respectively, 
and travelling in the direction towards the opposite end of th 
conductor. These constants wre found from the conditions at the 
conductor ends 


Ha aed =P (28.06) 
and equal 
A= - 6 [4(- 2.) — 4lz.bg(Lye™ Ye 
A= — Fy Ile) #(— ap Lye he, coe 
where 


D==1--yt(Lyen. (28.08) 


et es ee eee Se 
“)pootnote appears on page 216. 
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Considering that the quantity 1/D is equal to the infinite geo- 
metric progression, 


=the +y nye +..., fee309) 


Expression (28.03) may be written in the expanded form 


Jz) = Aste (iz pe?'"! —9(— 24) oo [9 2, 
— 4 (Lye gle, — zy e+ gt Lye gle — 
ze"... Jy ede (ee, — 2) Om 
—4 (Lie y(z — ze M4 
+ yh (Lyee(z, — z) eM... (28.10) 


The physical meaning of Expression (28.03) is seen from this. The 
first term in Equation (28.10) is the primary current wave which coin- 
cides with the wave excited by a concentrated emf in an infinitely 
long conductor. The second term (in all brackets) corresponds to the 
current resulting from the reflection of the primary current wave 

from the conductor end z 3 Za» and as a result of subsequent reflec- 
tions from the conductor ends which arise from this wave 
—Ao(—z)e "9(2— 2,)c*-" | The third term (in all brackets) corre- 
Sponds to the current resulting from the reflection of the primary 
wave from the end z 2 Zo and as a result of the subsequent reflections 
from the conductor ends arising from this wave --/epir.)e™ p(z,-- ze, 


It also follows from Equation (28.03) that external field (28.02) 
excites in the semi-infinite conductor (2:52.20) the current 


Jzy=Aletiz pe?! —4(—2,ye yz — 2,pe yh, (28.11) 
and in the semi-infinite conductor (--%<2<2,) the current 


Jz) =Slp iz yen’! —v (ze y(z, — zp er, (28.12) 


Comparing these expressions with the proper terms in (28.19), we see 
that the reflection of all the current waves at the end of a finite 
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length vibrator occurs in the same way as at the end cf a semi-infinite 


conductor. 


In the case of a passive vibrator (z,<27.). excited ty the piane 
wave 


fea Eye. wean — & cost, (25.17) 
the current also is represented in the form of the sum of waves (see 
[82]) 


J(z)=S (ew? —4F(z—2,) eftetenza 
—. (2, —- 2) eirstit (fy— a) + A, oz ~2,) om eng 


+A, $(2,—2) el", (28.14) 


where the first term corresponds to the current excited ty a plane 
wave in an infinitely lone conductor. Its complex amplitude S equals 


S= iwE,, 


2k atn* 3 in 


ytosind . (25.15) 


The second and third terms are primary edse waves arising as 4a consé- 
quence of the cut-off of the current sett | They are expressed in 
terms of the functions y + (z) and w - (z) which depend, in additicn 
to tne variable z and the parameters k and a, on the angle ® . These 


functions satisfy the relationships 
72M)=1, 4, (oc) =O, 
¥, (2) gas = ge (2)| ga = (2) 


4. (2)! sao= 9 (ace =1- corel) 


The initial values of the primary edge waves are sucn that tieir sum 
wz 


with the wave Se gives a current equal to zerc at the conductor 


ends. 


The last two terms in Equation (23.14) correspond to seeendary, 


Yaa 4 


tertiary, etc., edge waves, and have the same form as they da 
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transmitting vitrator (compare Equation (26.03)]. The unknown coeffi- 
efents x and k, are found from Conditions (28.06) and equal 


D 


A, oe a elt toler, (L) -- 
--b.(L) p(Lye 


(a- 2) odes, 
8 


Je 
A, — + Pa ial) (e (L) me 


— 5. (LY (Lp ery et, (28.17) 


Using equality (28.09), Expression (28.14) may be written in the more 
graphic form 


J(2) = S{e™ ~$.(2 —2,) efzetthtems) 4 
fe. (Lye tt yz, — 2) eer 
(Lye? (2 — 2, 
Hey (Lyet’ y(z,—2) Mr 
~—. (2,—2) eben anh 
ty, (Ly efter fy z —2,) eft. 
—y% (Lye'** y (2, —2) Pein 
yt(Lye az — 2, yy. (28.18) 


Here besides the wave sel¥z 


and the primary edge waves, which we 
talked about in connection with Equation (28.14), the secondary, ter- 
tiary, etc. waves diverging from the ends 2 = 24 and z = Z5 are 
explicitly written out; they correspond to the first, second, etc. 


terms in the graphs. 


Passing to the limit in Equation (23.14) when Zo * @, we find 
the current in the semi-infinite conductor (215 @) 


D2) 2: S-feiet gy. (2 ne zen ney (28.19) 

and, similarly, we find the current in the semi-infinite conductor 

(-2, 25 ) 
\ 


S(2) = S-feie? — 4, (2, -— 2p eer een, (23.20) 
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It is not difficult to see that In the case of @ passive vibrator 
the reflection of current waves at its ends occurs In the sume way as 


at the end of a semi-infinite conductor. 


Thus, the complex amplitudes of current waves ina thin, finite 
length conductor are proportional to the function: (2) and wy, (2) 
which monotonically deercase with an increase of 2 as a consequence 
of radiation. Let us note several properties of current waves in a 
vibrator. Each advancing wave in sum with the reflected viave excited 
by it gives a zero current at the conductor's end. In the case 
pa =1,2.3...) and D «© 0, a current resonance besins in 


L=2,--2, =a- 


the vibrator. 


The precision of Expressions (28.03) and (78.14) obtained by the 
method of slowly varying functions is different in various sections 
of the conductor. It is comparatively low near the conductor's ends 
(and in the vicinity of the point z = C of a transmitting vibrator) 
where the current waves arise, and where their complex amplitude varies 
rather rapidly. As the distance from these vibrator elements increases, 
the precision of these equations increases without limit. 


It should be stated that with a more rigorous atproach [7°, 20] 
the amplitudes of all the reflected waves will be determined by dir- 
ferent functions; however, the difference between them rapidly decreases 
with an increase of.the reflection number. The functions y(s) and 
v, (2) only approximately descrite these current waves, but on the other 
hand they allow one to effectively sum them and to obtain closed 


equations. 


Using the variational method ror the functions w(s) and y iz), 


we obtained the approximate, but on the other hand, simple cyuntions 
(see [83]) 
-! 
la = 
v9 
9(2) = az" 


la ae — E(2gxrie- ties 
1 (28.21) 


(Equation continued on next vase) 
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Yq (2)=: —y-——""8 
: * >" E@agsye {2 (28.21) 
where 
In(—t) = ie, Ini Hig. (28.22) 
& UF eosé ‘ 
. sae, g = (ka)’, 1,=9- 5 (28.23) 
and 
“ “ 
E(y)=Ciytisiy=—( Sat. (28.24) 
g 


The integral cosine Ci y and the integral sine si y are determined by 
the relationships 


@ 
Gy=— [Sita siya [Pat (28.25) 
8 


and are thoroughly tabulated functions. 


The equations written above for the current in a finite conductor 
are distinguishable by their visualicability, and they enable one 
to liken the conductor to a section of a transmission line in which, 
however, the attenuation of the current waves takes place, not accord- 
ing to an exponential law, but accordins to a more complicated law 
which is determined by Equations (28.21). In addition, the diffraction 
character of the problem is reriected in the equations. The conductor's 
specific features as a diffraction object are included in the very 
slow attenuation of the current waves. As a consequence of this, it 
is impossible to limit oneself to considering only secondary and 
terviary diffractions, and it is necessary to sum all the reflected 
waves. As a result of such a summaticn, a “resonance denominator" 
D appears which takes into account the resonance properties of a 
thin, finite length conductor. 
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§ 29. Radtation of a 'Transmittin Vibrator 


The radiation characteristic of a transmittinys vibrator tay be 
calculated from a known equation by integrating the currents in it. 
However, such an approach is not advisable because, as was indicated 
above, the precision of Equations (28.03) is different in different 
parts of the conductor, and is low near its ends (z = Z, and z = Zo) 
and tne point z= 0. The principle of duality gives more precice 
results. This principle leads to the following expression for the 
radiation field in the far zone [82]: 


é ith 
E, =H, =—-~-£-5-— S10), 
2 sin 8 In ‘fran 
E,=-H,=0. (29.01) 
The function 
1(9) =1-—~ Yo (2,) effet sens 8p }. (= z,) e ites it +cos 6) 
+ By, (Lye 9 4 By (Ly et" (29.02) (29.02) 


is connected with the current (28.14) excited in a vibrator by plane 
wave (28.13) by the relationship 


J (0) -= Sf (2), (29.03) 
The coefficients By and Bs do not depend on the angle @®. 


Expression (29.02) enables one to trace the formation of the 
radiation. The first term (one) is the radiation field of an infinitely 
long conductor excited by a concentrated emf. Propagating in the 
direction 09 = 0, this field reaches the conductor's end z = =, and 
— being diffracted by it —— generates a primary edge wave (the zeeond 
term). Ina similar way, the primary edge wave diverging from the 
conductor's end z = Za is excited (the third term). The last two terms 
in Equation (29.02) determine the waves arising as a result of subkse- 
quent diffraction (secondary, tertiary, etc.). The amplitudes By and 


B, of these waves may be found from the conditions 
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(9) -= f(#)-+9, (29.04) 


which means that the radiation of a finite conductor in the direction 
of its geometric extension must be equal to zero. These conditions, 
together with a consileration of the relationships (28.05) and (28.16), 
lead to the system of equations. 


B, + By (Lye = 9-20, 


By} (Lye + B, = y(z,) 0", (29.05) 
from which we find without difficulty 
B= tel—20 ihr gtze™ pe, 
By = 4, [4(2,) —~ 9th) (2, per ee, hanes 


Keeping in mind (28.09), let us represent the functions f(%¥) ina 
more graphic form 


f(B)==1— ya (z,yes 44 (2,) eX 
x (3. Here | 4 (L) elt! o, (L) pe itteces ® ie. 
+ o*(L) etl (L) ent b ia 
— Jo (zen Ml tern Yt (2 yettt- x 
x {v. (L) ea fhteens . = w(L) et 4 (L) e ite cos v4 
$ t(Lye* eg, (Ly enh, (29.07) 


where the secondary, tertiary, etc. waves corresponding to the first, 


second, and following terms in the trackets are explicitly written 
out. 


Thus, the field radiated by a transmitting vibrator arises as a 
result of multiple diffraction of edge waves at the vibrator's ends. 
Let us note in connection with this that the edge wave is diffracted 
by the opposite end of the vibrator in the same way as at the end of 
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a corresponding seml-infinite conductor. It fc not difficult to 
establish this’ by investiyating the radiation of a semi-infinite 
conductor excited by a concentrated emf. 


§ 30. Primary and Secondary Diffraction Ly a 


Passive Vibrator 


Let a plane electromagnetic wave fall at an angle % on a thin 
cylindrical conductor of length L = z, - Z4 and radius a (Figure 73). 
For purposes of generality, we will consider that the incident wave's 
electric field E, forms an angle a with the plane of the figure. 


) 2 
Then, its tanzential component on the conductor surface will.equal 


Et = £,,-e"*, (30.01) 
where 


Eu==Esin®,, El=E,cosa, w, =='— k’cosd,. (30.02) 


The current induced in the vibrator by this field was investirzated 
by us in § 28. As was already indicated above, Expression (28.14) 
which was obtained for it has a relatively low precision near tne con- 
ductor ends. Therefore, it is inadvisable to seek the fringing field 
by integration of the current. Let us also note that the frinzinr 
field found by such a method does not satisfy the principle of duality. 


We shall seek the scattering characteristic of a passive vibrator 
by starting from the follow.ng scattering picture which naturally 
follows from the previous results. An incident plane wave, beinz dif- 
fracted at the conductor ends, excites primary edge waves which are 
radiated into the surrounding space. Being propagated alone the con- 
ductor, each of these waves experiences diffraction at the oprpesite 
end of the conductor and excites secondary edge waves. The latter, 
in turn, generate tertiary edge waves, etc. The total fringins field 
is comprised of the sum of all the edge waves being formed during 
sequential (multiple) diffraction atthe conductor's ends. 
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In § 28 and 29, we notea that current waves are reflected fror 
tie ends of a finite leneth eonduct r the same as from the end cf a 
Semieinfinite conducter, and that the diffraction of these wavec at 
euch end takes place in tne came way ac at tre end of a semi-infinite 
ecnluctor. Therefore, tne prinary edee waves may be found from the 
croblem of: scatterine of a plane wave hy the semi-infinite conductor 
(21> @) and the conductcr (-", 2,). The sum of such waves gives tne 
primary JjJiffraction field 


ear 
EM H = -E se ,). (30.03) 
where 
re, ap 
. » a, 
t one z te? . Ege as dep cos *,) =. 
ar ose Peace Foust, - Feast) 
note 
i 'g 2 ‘2 2 _ e iteteur® + cos 8) 
—~"F (cost + cus U,) D(a cos #, s cos B,) a (30.04) ° 
<a f The function (D(w, wg) may be 
~" calculated by means of the rigor- 
5 % & ; ous solution to the problem of a 
.. semi-infinite vibrator (see [82] 


§ 3 and [83] § 4), and in this 
case, it satisfies the relation- 


Fizure 73. The incidence of 2 shin 
plane wave on a thin cylindrical 
conductor; % is the incidert 
angle. 
ss a. 
Pla, EVD (- vw -—wed- tn “In a4 1 


{a yea" 


va JR, ge) ee, | (30.05) 


However, hencefcrth it will not be necessary to have the rigorous 
expression for the functicn 4. Let us note that Equations (29.093) 
and (30.04) are similar to Expressions (€.13) for a strip. These 
latter expressions do not take into account secondary diffraction. 


cn 
cy 
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The secondary ede wave propassated from the endl 2 = fn Lo exe 
elted during adiftractton at this end of the primary current wave 


= Sp? (2 —2,)eerrn tena, (30.06) 


where by (2) we mean the functions obtained from the functions », (z) 
by replaciny ® by 4. For the purpose of cnlculatine the desired 
secondary waves, it is necessary for us, first of all, to find that 
external field which, when applied to an infinite conductor 
(--m~wer-a), would excite the current (30.06) on its section (z,%2%%), 


For this purpose, let us study the current induced in an infinite 
conductor by the external field 


ES st E..e' 3 (2 --2,) #(2—2,) = 


lwitn 2<2, 

ee (30.07) 
Let us assume that Wo has a small negative imayinary part (Im Wy < 0). 
We may rerard the quantity Eye 8a as a concentrated emf which, in 
accordance with equation (28.03), creates in an infinite eonductor 


fe a % (t\z--f) ez tik tem) dt. 


ae (30.08) 


ihe 


Therefore, im accordance with the principle of superposition, the 
total current created in the region *,%z<e by tie externol field 
(30.07) will equal 


° 


Eee! feet aber -t 
Jaye Pe NC ae yee NR. 
419 ‘, 
so « 
we - fF se. glitter, fa cgy ett elk yt 
= tha P--2, (30.09) 


The resultiny interral may be expressed in terms of the funetieons 
wiz - 2) and vo (z - 2))5 and the corresponding relationships derived 
in § 2 of [82]. As a result, we find 
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J(2)= ee ae $(z —2,)e 
Sh sta > Ie — 
cet of (2-2) 9 
21h sie? 0, le 
y4a cos > (30.10) 


Thus, it turns out that external field (30.07) excites, in 
addition to the wave y’, also the wave y. In order to excite a "pure" 


y? wave, it is necessary obviously to apply an additional external 
field 


Ei = @,8(z —2,), 


(30.11) 
such that 
be gk, i 
¢ “ ; y(z—2,)e8 9 + 
Sik sint 2 In yko 
Es y (ez =0. 
Nae (30.12) 
Hence, 
Eyse Tee 
yn— . (30.13) 
2ik sin Tr 


In order that the sum of external fields (30.07) and (30.11) would 


create the current (30.06), it is still necessary to fulfill the 
equality 


Ey 


eget = SOR Ne oc oes cat 


27k sin? 0, Io arene ii 2H sin Mla (30.14) 


MOE 


which determines the quantity 
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4 


Snake 
tha cos 


ai 
le qka sinh, 


Ey =— Ew 


Consequently, for the excitation in an infinite conductor (with 
) of current waves (30.06), it is necessary to apply the external 


Z> 2) 
field 
é 
a——; 
qha cos + 
EL = Exe os (z—2,)— e'"e(2z -2) , 
fa qha sia, 2k sin? -3- 


° . z—2)={ Iwith z2<z2, 
; Owith z>2,. (30.16) 


In a completely similar way, one may show that the external 


field 
Ia : , 
ya sin-> nye : jes 
Ef = Eu ~— 5" [| 812, -2) — ele, —2) 
“In jiasin *, H 2th eos? > 
eiz,.—2) <] | with z>2, 
; (Owlth 2<02, (30.17) 


excites the following current in an infinite single-wire line (with 


Zz < Zo) 
oa nea (30.18) 


ere 
= Se i Os ab ear 


Now let us study the diffraction of current waves (30.06) by the 
For this purpose, let us use the 


semi-infinite conductor (-«, Zo). 


Lorentz lemma [4] 
(a7. ian =o. (30.19) 

Here ji == —iop,:R -~R’) 18 the current of the auxiliary dipole with the 

moment p, which 1s located at point 1 with the coordinates (R, 9); Hy 
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is its field on the conductor surface, where the external currents 
55 are specified; E, is the field created by these currents at point 1 
(Fisure 74). 


The external current 35 is 
determined by the well-known 
equation 


es 


die 
IP = — > [nt] (30.20) 
Pigure 74. 
in terms of the electric field E 
on the conductor's surface. In 
view of the boundary condition 


E,-E,=0 (30.21) 
we have 
yea ae Ea (30.22) 


Furthermore, defining the dipole moment Py in terms of its field 
in free space (at the point x » y = z = 0) 


‘ 1. te 
Fae hPa see (30.23) 


and changing from the magnetic intensity Hy to the total current 
J= Aye (30.24) 


induced by the dipole in the conductor, we obtain from the Lorentz 
lemma the following relationship: 


‘ ptsinn ef@® 
E,, =. = —. : | Ef J(z) dz. 
) 


“we, OR (30.25) 
If the dipole p, is moved to a distance R >> z, - Zo» then the 
field radiated by it may be investigated on the section Za - Zy of 
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the semi-infinite conductor (-a2, Za) as a plane wave. Then the 
current induced in this section of the conductor will be determined 
by the equation 


J(2z) = S [et — ef*8, (z, — 2) ef -], (30.26) 
where 
Sen es, , w=--keos. (30.27) 


2é 
2k? sin? Ala yea sia’ 


We will select the quantity Zo in such a way that, at a distance 

Zn - 2 from the conductor end, the reflected current wave would be 
practically equal to zero (%+(z,—2)>90). Substituting the function 
(30.26) into the right-hand member of Equality (30.25) and taking for 
the quantity Ey the external field (30.16), we obtain 


Eg = Hy =——-g— X 


2i 
2 sin 9 Io qyhasin® 


- 


x i" +f ES fe ae (2,—2) elt tet ih (eral | dza . 
a 


[es ae %, (L) eleeaninty + 
tae tio oe, 


& 
a f (wtere) % 2 2, ie 
+ Eu “lwePor _ E,,e" pita § Shel at ter-a)ée} (30.28) 


An important feature of this relationship is that the interration 
is performed here not along the entire conductor (-@, Zo)s but only 
along part of it (-«, 24), where the function ¥+(2;—z) describes the 
current with good precision. The integrals here are calculated the 
same as in Equation (30.09). As a result, the field radiated by the 
semi-infinite conductor (-9, Z) will equal 


Ew = H,, = 


2100 Ia - as sind 


(Equation continued on next. pare) 
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E,,e! (34% 


: ise, ‘ Po ae 20 4 
x hi [ee Le ics wp eos a” 
sin?-) la ---— — 
4! BE, ,eith sites seen (L) yka sin 5 (L) 
y Relired Horas od AS eceed -- =< . 
“ik(cos #-P cos i,) | *+ ope ae t= (30.29) 
yha cos + 


The terms in the braces having the phase factor elW correspond to 
the desired secondary wave diverging from the conductor's end z = Z5* 
Using Equations (30.15) and (30.15), thls wave may be represented in 


the form 
(2... tar 
ED (e,)== HO (z,) = -— 2d Semen, (30.30) 
2 sin '3 ka aad ; 
where 
fergt +a 
g (2, j= =—- —-1eo-—_,— x 
Rin h,(cosd + cos*,)in 3 Slag, 
x cost * > cos’? =z 1a——_,- ag We ()— 
7 yha — 
— sin® *e. sin * 7 la —y¥ y°(L)] . 
15a sia (30.31) 


In a similar way, let us find the secondary diffraction wave 


beine proparated from the end zg = 4° In order to do thls, it is 
necessary to investirate the diffraction of primary wave (30.18) at 

the end < = 2) of the semi-infinite conductor (z,*72< 0). In this case, 
the principle of duality leads to the followin=s relationship 


Sk ; isin @f@® e 
Bay Ha ie Re [f J(e)dz. 


(30.32) 
which, after substitutine the funetion (30.17) and the current 


ink? ae = 
J(2= in Sen umn ta -fe7—p_(2-- z,) eiv't pik tz a 
2k? sin? Ala yha siod ; (30.33) 


, 
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in it gives us “he field radiated by the semi-infinite conductor 
(2)> o), The wave radiated by the conductor's end is the desired 
secondary edge wave and may be represented in the form 


: (2) (2 ian : 
ED (2,) = 12 (2,) = -- EEN, Sg iteveon (30. 34) 


2sto 0 In mae R 


where 


ry} Eelt ete Hal ’ 


x 


8 (2,)=—- 2 
A sind, (cos ® + cos 3,)lo eT Sing, 


x ee Se. sin® + In —+~— y- (L)— 


qh sia 
— cos? 22. cost—"- in —— 9° (L) | 
me “- r tha cos-5- a J (30.35) 


Otherwise, this expression may be written directly by replacing, in 
Equations (30.30) and (30.31), z, by 24> 8 by z—8@ and & by =—%. 


§ 31. Multiple Diffraction of Edge Waves 


The secondary waves (30.30) and (30.34) which were found are the 
waves diverging from the ends of the semi-infinite conductors (-°, Z5) 
and (Zy5 *). If one excites an infinite single-wire line bv the 
external field 


E* 22: ig (2 24), (31.01) 


where 


&y'== ge (2,) | 7 TNS ee 2k ape na ett — ikz, cos &. 


then a spherical wave arises which with @®=g coincides with wave 
(30.30). With the excitation of an infinite line by the concentrated 


emf 
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Et= G8 (z —2,), (31.03) 


where 


t 
la $5.0) 
g?= gi (2.)|,_ = E ae = ett —its,eon me 


sla, la as aah, (31.04) 


a wave arises which coincides with wave (30.34) when @=0. It is 
not difficult to see that these external fields actually excite in an 
infinite single-wire line current waves which are equivalent to the 
secondary current waves in a passive vibrator [that is, equivalent to 
those.waves which are expressed by the first terms in the brackets of 
Equation (28.18)]. Therefore, the tertiary waves may be investigated 
as edge waves radiated by the semi-infinite conductors (Z,> @) and 
(-®, 25) with their excitation by the external fields (31.01) and 
(31.03), respectively. From Equations (30.25) and (30.32), we find 
without difficulty the total field radiated with the indicated 
excitation by the conductor (Zy> a) 


sf fat 


E = H ooo a” [er $_( L) elmerbiaty 

, asad tay & (31.05) 
and the total field radiated by the conductor (-=, z,) 
Ey= y= panies cee [elm 4, (L) efttet!3L), 

‘ala Mle eosin (31.06) 


As a result, we obtain for the tertiary waves diverging from the 
ends Zy and 25 che following expressions: 


G0) (2) elt en itt con 8 


EO (2,) =H (2) = 


a (31.07) 
asia din eo siae 
Ev) (2,) taR 
Eo 2 y= HM (2 = a art xe sahtecos © 
(20 » (2,) Sia bia atte P (31.08) 


where 
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Gi (2,) = — Bo (Lye, | 


8 (z,) = — 8-9, (Lye. | (31.09) 


In tne dlrections toward the cpposite end of the conductor, these 
waves are equivalent to the radiation of an infinite sinzleewire line 
excited by the external fields 


E* = &.3(z —2,), Br = B0N(2,)| os ; (31.10) 
EL = Gy -8(2,— 2), 6 = 82,)], (31.11) 
Consequently, the quaternary waves again may be investigated as 
edre waves radiated by the semi-infinite conductors (-», Za) and 


(zy5 ©) with their excitation by the external fields (31.10) and 
(31.11). Using the reciprcecity principle, we easily obtain 


ise 


E*"(2,) _ H(2,)= cM(2,) em ite cos . 
®. r= O's = 2 R 
2stn Btn Yasin 
(*) (2. lth —ikr, - 2 
EQ) (z,) =H (2,)= ean eS Sete, (Si.12) 
din In So sinw | 
where 
6 (z,) = — &-9_(Lye™, 
(24) = — 6-9, (Loe, (31.23) 
In a completely similar way, the nibh order edge waves 
Ey eH (2) = — Site _ aR ite cos 
J me ey 34 R ° 
osiad lo Yeo sin & 
abe + ier ‘ 
é (z;) a H'(2,) elie — ee = a iRéaces & 
2sin 9 lo~ RS cin t (31.14) 


are found. Here 
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6 (z,) = — a (Ly e!**, ) 


GI (24) — G9, (Lye | (31.15) 


and 
= 8'"(z,) on Bir) = Bt” (z,) te” (31.16) 


Thus, the field arising with multiple diffractions (starting with 
the second) may be represented in the following form: 


SAEs e+ Ey ed = 
aw? 
’ oft nl a —ibz, cos 8 
=, Lam eae + 
2sla Pia Tha sine nw? 
m he 6a) “_ - ke, cos 8 . 
. +2 (z,)e j | (31.17) 
where 
se F 
¥ ore) = 
awd 
= BN (2) +189 (Lyett— Arye, ater) 
yee) = 
a=? ' a 
= GI 2) (Ep 9 (Lye — BJ Be, (31.19) 


and the functions &(z,,) and 8 are determined by Equations (30.31), 


(30.35), (31.02) and (31.04). We will not write out here the rather 
unwieldy final expression for this field, but we wil] proceed with 
a calculation of the total fleld scattered by a vibrator. 


§ 32. Total Fringing Field 


Before beginning the derivation of the expression for the 
scattering characteristic, let us make the following observation. The 
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functions 9 which enter into Equation (30.04) satisfy relationships 
(30.05), and may be found by factoring. However, our investigation 

of the Successive waves arising with diffraction at the conductor's 
ends was approximate. Therefore, it makes no sense to use the precise 
Expression (30.04) for the primary field. We shall use the approxi- 
mation expressions for the function ¢. 


®(— kcus >, — & cos %)=In ace may, 
tha sin ~3- sin > 
P(kcos), kcosd,) =: In -——--. ee 
ha cos “9° cos eS | (32.01) 


which were obtained by the variational method and have a precision 
which is sufficient for our purpose (see [83]). More precisely speak- 
ing, we will use approximation Equations (32.01) in conjunction with 
the rigorous Expression (30.05), and we will set 


é 


la eons Ya 
\ , yea cos “x cos - 
Di Feas’, easy 
1 10 ckasin® '" Yea sing, 
In -— -- -ee ' 
qha sin > sia > 
D(kcosd, Reus) Ni. DAE ia 
* jka sia’ '" yéasin', (32.02) 
Then the primary field will equal 
ED) |. ge EE a 
2(cos9 ens t)in mera tn AR 
Aa sing’ yka stn &, 
x ctg al) cty + In ele tb —+- eae (cos Bees %) 
tha cos ~~ cus ig 
mtg bg fin fg grit tent 
q4asin-> sin + (32.03) 


Now summing Expressions (31.17) and (32.03), we find the total fteld 
Scattered by a passive vibrator in the form 
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aha 8 F(0, %), (32.04) 


where 


2i 
F(8, 8,)= —_____4____...—-,,— X 
(cos ® + cos 9,) sin sin &, In -saae waek 


ite, (cos 64<00 8.) 


® 
X ]cos* > cos? $s. In 
yha cos “yr 608 —H ‘ 


-~sin? 5 8 - sin? 3. ‘In = jee Sree 9 1 
jha sin = sn = 
wet Poa 4 a, i 
‘Fe | sin? => sin? } In nga 
yaa sla =~ 2 


® 
a costs cos? In ——— 
tea cos 


eh s,cowD, + 2, cus br 4+ 


+e" |sin® > sin oe ig 1  — 
qhe sia > 
af In i Y \" (z,cos@+s, cco Og 


id 
qta cor 7 


cos @ + cos &,) la — 
( + Y 7 iad 10 pith — its, cosy 
2D 


a e elt cos ey o enits cos o_ 


(cos 8 + cus 8,) In yra 
ae 40 


10 ike, 6, ~ike, os 8 
—Ye cos Iv, e i ¢ }. . (32.05) 


Nhl 50. dal—ike, cos 9, 
er [pve ~~ 


in which all the functions ¥, and ve have the argument L. The result- 
ing expression, despite its complexity, has a clear physical meaning. 
Actually, the first term in the braces corresponds to the primary ede 
wave radiated by the conductor's end Zz = 243 the second term corre- 
sponds to the primary wave radiated by the conductor's end cz = Zo 

The terms included in the first pair of brackets refer to the secon- 
dary wave departing from the end z = Zi» and the terms in the second 
set of brackets refer to the secondary wave departing from the end 

Z = Z5. The remaining terms describe the sum of all subsequent waves 
arising with multiple diffraction and have a resonance character. 

The resonance begins with L = 2, - 2, % n-A/2(n = 1, 2, 3...) when D2 0. 
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Another important feature of the scattering characteristic is 
that it satisfies: the reciprocity principle — that is, it does not 
change its vaiue with a mutual interchange of ® and . One may 
also show that the vibrator does not radiate in the directions along 
its axis, and it does not scatter electromagnetic waves with glancing 
lrradiation, that is, 


F(0, 8))= F(z, %)) = 7 (8, O)= F(%, =) == 0. (32.06) 
Furthermore, using representation (28.25) for the functions y 


and ¥,, we obtain the following expression for scattering character- 
istic (32.05) in the direction of the mirror-reflected ray (0=1—0} : 


k 
F(z—9,, 8.) = — —“—+ 
in yee slo a, 
(VE ( 2k! sin? +) +E ( ORL cos* +) 
+klL —— -———,, le + 


‘(10 rasa) 
‘ ( 
= fin + 


+ eel A aan sin in 8, 
(sina, a, In Tea sin am) 1 


® fs 
+ —v cus* =—~In ann tt ; yarn! cos ody. 
tka sia 


J 40 cae & i 10 qth (feos 8) 
+/ 3 ¥2 sin’ — In ——5— | 4_-e 
pie CoB Zo > 


yea 2thh 7.0, ikk (= a +0 7,0 rd 
Dp ® [v,ve oer %— 2 14,— 


i 
WO ea ies 
my ee gt yet Ee ah. (32.07) 


With glancing irradiation of a vibrator, when 9=0 or %=r, it 
follows from this that F(z, 0)=F(0, 2):-=0 


Now assuming that = in Equation (32.07), we obtain the 


relationship cgi Nsg “a PERL 
ee rere t or NG 
BET ee rn = 2(in 3) 
“™ ya k 7a 
4 t . ‘ i oat 
+ pclae |e tae yen 
8 yaa 


Equation continued on next pare 
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(32.08) 


which characterizes the reflected field magnitude with normal 
irradiation. 


Let us also write the expressions for the function F(®, 8) which 
corresponds to the radar case when the observation and irradiation 
directions coincide (8=) 


F(8, 0) = 
a eee i ae j rote 


2i 
(sin ain Jka sin z) cos 8 
f ike, CO 4. 


4 A 
—sin'--In----- =-e 
. ysasin®-,- 


-f2 i pln —+-5-p_(L)— 


yea sin ~~ 


tha cos? ~~ 


— cos! + In ee ty Ee (pee) 08 Opie 
yea cos > 
— 2st lo atl) ty, (L) y_(L) e ih (2, +20 cos ide 


TSE nie th (Lye 


yt (Lye atest 4 ey . (32.09) 


We may show that when Iz Equation (32.09) leads to Expression 
(32.03). 


The scattering characteristic (32.05) which was found above was 
octained by summing all the waves formed with multiple diffraction. 
Such a method is very graphic, but somewhat lengthy. Cne may arrive 
at the same result more quickly if it is assumed vhat the edge wave 
diffraction process at the passive vibrator's ends takes place, start- 
ing with tertiary action, the same as at the end of a transmitting 
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Vibrator.  Pherefore, the pasttive vibeater’s: seartertiae abaya merry 


be ddreet Ty soupht ti the Corn 


K (9, 9%) -- a E08, 8). 


LLL? Sate Bue Viet TB 


Wheres 
119. 9) x 
woe cos 8-f case, 
~ Jeus? e-custe In 4 ge Hee gos Bjvoe By 
“ yhe con Pa oe 
wae ain? é sin? ”* in ‘ . v Try Oo). 
2° 2 ®, 
jtaska-., sta . 
‘ ; | a 
-fe "| sin! yo sin? In ' a Y_tL)-- 
yaa sin 2 
we 7 wt te I i a? 1 2 Rie cow Ot owen 8 
eon cosy Ta iy ¥,! at - 
yha cos “2: 
et i | ra) é ° 
‘fe’ [sine y MAP? In » th) 
yaa slo 2° 


® 8, 7 . rr) «te 
“tus? scons? S In - : ; ¥, (L) |e Mitatee/ @ ved = |- 
yha cos 2 


| Cy : Le te, cond | Cy, (1) eoobes ces . 


The bast twee Germ: tao Bquatlom qatsd dd sree tte sum on: 


(se.te) 


eben : 


Waves start fires wlth the tert Lary waver whdeby cpte Peopreited Prseer ttre 


eortduetorls: emda ow = fy atid soe yy Penpeet tvedy. The 


rbP Co rime betermntned from dhe cond it fon 


Pay pee ay 0, 


Whteh Tends te the aystem or eo quart lor 


’ ‘ . 1 r Aad the ost 
Cy e fy ? UM as rire 7 


ba teak, 


a5 ott t : e coe 
©, bP eyglbe vin i bbe 


Cron wirbeh 


eetechintis 


! 
(ny 1) 
toes ha} 


. kd e al — brgosd, 
Cyst age in fp ye’ atti 


“D yha 
ee re eof ticoste }. 
t cee ‘ ©, sal~ike cord 
C.=-a5 0" ‘lay lve °— 


— yf ere onte (32.14) 
§ 33. A Vibrator Which {ts Short in Comparison With 


the Wavelength (a Passive Dipole) 


The theory of plane wave scattering by a thin cylindrical 
vibrator which is discussed In this chapter is based on a number of 
physical considerations. One good aspect of this theory is the fact 
that its precision increases with the length of the vibrator, since 
the current waves whose diffraction we are investigating are expressed 
more clearly, the longer the vibrator. SHowever, one may also show 
that for short vibrators, the length of which is small in comparison 
with the wavelength, the equations derived by us have good precision. 


It is clear that a vibrator which is short in comparison with 
the wavelength acts as a dipole, creating a fringing field 


(ag 
E, =H, =~ bp, ge sin 8, (33.01) 


where the dipele moment p.. may be enleulated by solving the electro- 
static problem. This dipole Moment depends on the dimensions and 
shape of the vibrator. In accordance with [92], the dipole moment of 
a cylinder in a uniform electrostatic fleld E, equals 


L 8 
pe BiN gs (33,02) 
where PCJ) $s a dimenstonless funetton 72 = £/3n whieh is shewn Jn 


Firure 75 by the continuous curve, With 2 >> L, one may ealeulate the 


function D(z) by means of the asymptot Le expansion 


Dil) = t(atart \ Q, = 2( inal 2 x): 


my 


PTD =H =2 32059271 


If in this expansion one limits oneself to the first term, then 


Dij)=——_—.. 
3(tu—-5) (33.04) 


The results of numerical calculations based on Equations (33.04) are 
shown in Figure 75 by the dashed curve: we see that the latter 
equation gives a good precision already with 2 2 9. 


an ¢ 
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Figure 75. Graph of the function D(z) which deter- 
mines the cylinder's dipole moment. 


Thus, the dipole moment of a vibrator which is short in compari- 
son with the wavelength equals 


RSet —a tr {'+0[ (is) }}. (33.05) 
a 3 


and its scattering characteristic must have the form 
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A°L® sindsind, -? 
F (8, ne aT rt 1-of(ine)- I}. reer 
In this section we find the first two terms of the expansion of 
the vibrator's scattering characteristic F in reciprocal powers of 
the large parameter L/a (with A + @), and we compare them with 
Expression (33.06). We shall limit ourselves to the case §=% = +, 
when the function F is described by the simpler Equation (32.08). 


With small values of the argument z, the functions yw(z) and 
+(2)= 94), <= #| : [see Equations (28.21)] may be represented in the 
form 


$(2y21— £960 + faq). | 


b(2) = 1 ~£(2— £0) oe 
Ha) 1-H FO) 40(sig). | (33.07) 


The functions g and g included here are determined by the equations: 


g(2)—g (0)=I0 AF pel As, #(0)= in (33.08) 


a 


@ 

Pe ye | Lae To Oe tials A Oe 

€(2)— (= In eM | Eads, FIN FE. (33.09) 
& 

Let us note that Expressions (33.07) completely agree with the 
corresponding terms of the aysmptotic expansion for the functions ¥ 
and wy, which may be obtained from the initial integral equations which 
determine these functions (see, for example, [81], § 4). 


Limiting ourselves in the expansion for the functions (z) and 
$(z) to terms of the order of (xz)3, we have 
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(2 = 1— -—-— {ike ae 1)-+ 


ln . 
ya 
yee 3 / A bd iy 
-+- eer cins ey ae ap Ret Int) (33.10) 


and 


a= I rae [ine (103 ase = 1)+ 
+*F (in2— +) — it (in MA ~))- (33.11) 


i \-2 
In addition, terms of the order (19 3) are omltted in Expressions 
Now if we substitute these expressions into 
i\~9 
7) 


(33.10) and (33.12). | 
Equation (32.08), then one should omit terms of the order ( 


Therefore, the function a may be represented in the form 


it. 
gikl 
r(8-9)= le e004 
waibaie en ~ 29 (L)1o va et 
i 
2ta = — 
Tee asp tt}, 
Fite Pe (33.12) 
Furthermore, taking into account Equations (33.10) and (33.11), 
we find 


ue [ts as ~ E(t |= 


=F (10 —ine SF), (33.13) 
nde : 
sit i. (33.14) 
25 yale at 2 rence 1 BS tn at 
~ikk (in — t) aL? (pine 
ms +ikte (fine se). meee 
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and finally 


i 
2la-—— 
. yka _ i 
Pawent FL = 
=(1 i-e— — BL i RL) in 


— Fite vas 2) EE pn ME OE (1410 2) 


pi ge Ln HL ines i *) . (33.16 


(33.16) 
Using these relationships, it is not difficult to show that 
F(>. $)=*o ryt O[(in Sy 
($4) fut zptolesr yy asan 


The equation which has been found may be rewritten in the form 


F( +)= 24{ In a —¥y {i+0|{n- z) }}. 


(33.18) 
It completely agrees with Equation (33.06) which follows from [92]. 
This result confirms the correctness of scattering characteristic 
(32.05) calculated by us, and shows that it is applicable for vibrators 
of any length. 
§ 34. The Results of Numerical Calculations 
The function F(0, 0) enables one to calculate the integral scatter- 


ing thickness S and the effective scattering area o of a passive 


vibrator. The integral scattering thickness is determined by the 
relationship 


hd P 
ieee) (34.01) 
where 


(34.02) 
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is the energy flux density in the incident wave averaged over an 
oscillation period, and 


‘ 
P= Re [ntenejas— +. Esin®, Rel J(aye'te dz (34.03) 
&% 
is the value of the energy scattered by the vibrator into the surround- 


ing space averaged over a period. Since one may represent the fring- 
ing field in the far zone in the direction ®=2—® by the equations 


aa — ik. ft ° ZC 
Ey, =— sind 7 0% te (34.04) 
ay 
and 
yd . 
Ey=H, =~ ESF (a=, 8), (34.05) 


then, having determined from this the integral 


fuera = aie F(x—8,, 4,) (34.06) 
we obtain 
=e cost a-lin F (2 — 8. 0,). (34.97) 


Calculations of the quantity S/L? (with a=0, %=:-)- ), performed by us 
for vibrators with the parameter ae errr taking the values yxy = -0.05 
and x = -0.1, are found to be in agreement with the results of 
Leontovich and Levin [85]. With x = -0.1, our curve (the dotted line 
in Figure 76) 1s only slightly displaced in the direction of longer 
wavelengths and gives slightly higher resonance peaks. 


The effective scattering area o according to the definition 
equals 


(8, 9,)—= See (34.08) 
where p is the known quantity (34.02) and 


Figure 76. The integral scattering thick- 

ness of a vibrator as a functior. of its 

length (with normal incidence of a plane 
wave). Curves 1 were calculated by Leontovich 
and Levin [85]. Curves 2 were calculated 

on the basis of Equation (34.07). 


F (8, Ra)? 


represents the average value of the energy flux density scattered by 


the vibrator in the direction f . Consequently 


FTD-HC-23-259-71 


Eytan get (34.09) 


.2(0, 0) =" cus? a-| F(, 8) /?. (34.10) 
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If the receiving antenna operates with the same polarization as 
the transmitting antenna, then Phe corresponding value of the effec- 
tive. area aac equal 


ee ; . | 9,00, 9) = 7 casts. | F(8, 84) (34.11) 
ae ee mG 


In } ‘the case pf adie: when the transmitting and receiving antennas are 
combined arid: the polarization is arbitrary, the vibrator's scattering 
properties ‘are characterized by the average value 

ne ere ae ee i 


t 
hate ee ood 


fi bia be P tany. a, (2, 0)da= FF oo (34.22) 
In Figures 77 and 78, the results of calculations performed on 
the. basis. of Equations (34.12) and (32.08) for the case of normal 
irradiation (*=+) are represented by the ‘dotted lines. Figure 77 
illustrates the ee of the function o on the quantity kL with 
a given. value of Qa, =2ia = 156, —- that is, when the ratio of the 
vibrator's length to its diameter equals L/2a = 452.. In Figure 78 
the graph of the function o is constructed as a funciion of the fre- 
quency f = /d+107° (in megahertz) for the prescribed parameters 
L = 5 cm and a, = 15. Here the curves plotted by Lindroth (79] are 
drawn with a continuous line, and the curve in Figure 77 calculated 
by Van Vleck et al. [86] is traced by the dash-dot-curve. 
® 


The curves of Lindroth and Van Vleck were calculated by inte- 


Ce ia 


Pe a oe : e 


. 


grating the current which is found as a result of the approximate 
solution of the integral equation. However, this procedure was per- 
formed in [79] and [86] in a different way. Lindroth obtained an 
expression for the fringing field in the form of an expansion in 
reciprocal powers of the parameter 2_. The expression includes terms 
of the order of a In [86] a different kind of approximation was 
used which led, as can be seen from Figure 77, to rather rough result 
especially in the resonance region. Our curve (the dotted area) 
agrees almost everywhere within the limits of graphical precisicn wit 
the curve of Lindroth. A noticeable divergence is observed only in 
the magnitude of the first resonance peak. 
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Figure 78. The effective scatter- 
ing area of a vibrator as a 
function of the frequency 
f = c/r-10~° (in megahertz) with 
normal incidence of a plane wave. 
The designations are the same 
as those in Figure 77. 


2 
4 In Figure 79 and 80 radar 
diagrams are constructed for vi- 
brators of a length L = 0.5. and 


Figure 77. The effective scatter~ Ll = 2\ with the specified value 


a na fie aS rpg ete os L/a = 900. Curves 1 were calcu- 
unction o s length w 
normal incidence ge plane wave, lated by Tal using the variational 


Curve 1 was calculated by method [87]. Curves 3 were ob- 
Lindroth [79]; curve 2 was 


calculated by Van Vleck [86] by tained by the method of induced 


means of the method of integral emf [86]; curves 4 were obtained 
equations; curve 3 was calcu- 
lated on the basis of Equation by the above-indicated method of 


(34.12), Van Vleck. The results of calcu- 
lations based on our Equations 
(34.12), (32.08) ana (32.09) -age 


a shown by curves 2. 
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Figure 79. A comparison of the diagrams 
for the effective scattering area of a 
half-wave vibrator calculated by various 
methods. 


Curve 1 was calculated by Tai [87] by the 
variational method; 


Curve 2 was calculated on the basis of 
Equation (34.12); 


Curve 3 was calculated by the method of 
(Beis: emf (in the work of Van Vleck 
> 


Curve 4 was calculated by Van Vleck [86] 
by the method of integral equations. 


In the cited references, the fringing field was calculated by 
the direct integration of the current. In order to determine the 
current, various approximation methods were used. In the variational 
method [87] a functional was constructed for this purpose which was 
stationary in respect to small current variations. Then the current 
was sought in the form of some function containing undetermined con- 
stants. These constants were found from the condition of the func- 
tional's stationarity. This method enables one to rather easily 
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Figure 80. Diagrams for the effective 
Scattering area of a vibrator calculated 
by various methods. The designations 
are the same as those in Figure 79: 


L 2 2A 


obtain the first approximation; however, its results, especially for 
long conductors, may depend in a substantial way on the form of the 
trial function. In the induced emf method [86], the current is 

sought in the form of a combination of trigonometric functions with 
unknown coefficients. These coefficients are determined by using the 
law of conservation of energy. This is the simplest method, but it 

has a number of serious defects. Thus, as a consequence of incorrectly 
accounting for the current component having the incident field phase, 


FTD-HC-2 3-259-71 rad 


it leads to inaccurate results in the case of odd resonances (cspe- 
cially for long conductors), and it does not give the displacement of 
‘the resonance peaks from the values \ = 2L/n (n = 1, 3, 5....) in the 
direction of longer wavelengths. 


The results obtained by us are also approximate. However, our 
Equation (32.05) satisfies the reciprocity principle, and is applicable 
for any length vibrator. For very short vibrators L << \, it changes 
into the asymptotic expression for the scattering characteristic of 
a dipole (see § 33). For vibrators with a length of several wavelengths 
(Lx nk, n #1, 2, 3, 4), Equation (32.05) gives satisfactory results. 
Calculations performed on its basis for radar reflection with normal 
irradiation agree with the results of Lindroth. Good agreement is 
also observed with the results of Leontovich and Levin for the inte- 
gral scattering characteristic. With an increase of the vibrator's 
length, the precision of this equation increases, and in this way it 
is favorably distinguished from the equations proposed for the 
scattering characteristics by other authors. 


Moreover, the divergence between the various approximation 
methods indicates the necessity of performing rather detailed calcu- 
lations based on precise methods, for the purpose of evaluating the 
actual error of the approximation methods. Such calculations may be 
performed, for example, by means of the method discussed in 
References [88, 89] or [91]. 
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FOOTNOTE 


Footnote ‘1) on page 1/7. Let us note that one may r2fine 
Equation (28.04) by multiplying its 
righthand member by the factor % 


(usually % x 1) caleulated in 
Reference [84], 
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CONCTUES LON 


In this book, the gotutton of ao number of di rftract Lon prob tems 
wis obtained based on the approximate: constderit ton of the bela pere 
furbation fn othe vietolty of oa sharp bend of the surtaee oon athicurp 
edyve. Equations were derived for the seatter diay chiupaetertot fess. or 
In certain cases (Chapter WV), for the radar rerleetion thleknegses 
With a speetffred Lrradtation direction. The expressions whteh have 
been found have acelear physteal meantiag. They satdery the ree dproe fly 


prineipte, and Chey are eonventent for mak dae ea leulat toms. 


The results obtatned enable us to form a more complete: concept 
oP the applfeabliity Tintts of the physteal optles approaeb. ot ts 
Ustlly eustomary ftooassume that this approaeh wtves veltable resutts 
ently PP the bedyts dimens fons are Larye da eompaeisen whth the waves 
lemeth.  Sueh oan optiden ds based on the Collowlime: areument. The 
Physteatl opties approach eonsfders enty the radtat fom trom the 
Unt form part oof the current, and does not Pnelude in the eateulat fore 
the nonunfPorm part of the eurrent whieh ds eoneentrated tao the 
Vietnity of the bends and the sharp edmes.  ‘Pherefore, when the body fs: 
dimens fons are consdderabhy larrer than theo wavetensth, Che nonmuat form 
part of the current: cecupfhes aorvedativedy amath pert or the body fas 
Surface. Therefore, ome would think that bts Pot bucmee woudd be 


smell. 


Mut fn oaetuabit¢y Tt turns out that the rPediatildty of phys tect 
opt les. results depends substan! tally, not ondy on the bodyvfs dimer 
Sfons, but alse oon the bodyfa shape And the Peird dat fom and: oteseryert pon 
dfreet homes.  Forcexample, with the atlanme tie fie dderee ef a oqweaywe cor 
the Plat Pace oP a body, the edee some ocenpled Dy fhe pretrun d Ceeper: 
part oof Che current fs constderably broadened aad the oetfeet oor this 
current: becomes substantial. Therefore, plysteal Opi des rivet quadli- 
tatively fnecorrect results: for the Cheld seat Cered Py Pdat ptsrtesns 
with plane tae Treeadtiation fndependent of the eatho between their 


dTitens fons and Che waveleneth. The effeed of the nmonnund form pat oot 


MPD = tee 82 07] 


the current becomes noticeable also in those directions where, accord- 
ing to physical optics, the frinving field must be equal to zero or 
have a small value. 


The problem of diffraction of a plane wave with its incidence on 
a cone along its axis (§ 17) serves as a clear example of how impor- 
tant the above-indicated factors are. Although in this case the non- 
uniform part of the current, concentrated near the cone’s vertex, has 
practically no influence on the scattering, nevertheless, the physical 
optical approach gives values for the radar thickness which are tens 
of decibels smaller than the experimental values, even with large 
dimensions of the cone. The deciding factor here is the nonuniform 
part of the current flowing in the vicinity of the sharp circular base 
rim of the conical surface; the nonuniform part of the current has an 
especially large value for sharply pointed cones. 


Another interesting example of a similar nature Is the scattering 
of a plane wave by a finite paraboloid of rotation (§ 18) where the 
physical optics approach leads to qualitatively incorrect results. 

The effective scattering area calculated in this approach turns out 
to be a periodic function of the paraboloid length, and with certain 
lengths it becomes zero which most certainly does not correspond to 
reality. 


The investigation of the diffraction of ede waves shows 
(Chapter V) that for flat plates one may limit oneself to considera- 
tion of secondary diffraction, if their linear dimensions are one- 
and-a-half to two times larger than the wavelenrth. 


Let us note that we attempted te obtain equations for the 
scattering characteristics which would possess physical visualiza- 
bility and which would be convenient for makin calculatlons. In 
keeping with this, we were oblised to Introduce various kinds of 
interpolation equations and simplified equations which satisfy the 
formulated requirements, but In return are not tn the general case 
the dominant terms of the rtyorous asymptotic expansion itn powers of 
the small parameter \/a. 
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Our purpose was not to calculate the current on the body's sur- 
face, the field In the near zone, or the Interral scatterine thickness. 
These questions are investigated in a number of other works based on 
the physical theory of diffraction which were already enumerated in 
§ 25. In them, in particular, the first terms of asymptotic expan- 
sions in powers of A/a were obtained for the integral thickness which 
characterizes the total power scattered by a body. !lowever, in these 
works, as a rule, equations are missing for the scattering character- 
istics wiiich are valid with any direction of irradiation and observa- 
tion. Therefore, the results of this book and the indicated works 
mutually supplement one another. 


At present, only a limited number of diffraction problems have 
ylelded to theoretical studies, as a result of which experimental 
studies of diffraction by various bodies have taken on sreat importance. 
In Chapter VI an experimental method was discussed which enabled one 
to Isolate in a "pure form", and to measure, the field from the non- 
uniform part of the current excited by a plane wave on a metal body 
of any shape. _ In the same chapter, it was shown that the well-known 
phenomenon of depolarization of the wave reflected from a body which 
is found in free space is produced by the nonuniform part of the 
current, or, in other words, by the surface distortion. 


The investigation carried out in Chapter VIT for the problem of 
diffraction by a thin, finite length cylindrical conductor represents 
a natural development and completion of the method of considerins the 
multiple diffraction of edge waves which was applied in Chapter V. 

In Chapter VII equations were derived for the scattering diarram whieh 
are suitable for vibrators of an arbitrary lenrth with anv irradiitien 
and observation directions. 


The results obtained in this book show the frulftfulness of 
physical diffraction theory, and enable one to arrive at the solution 
of other more complicated problems. Such problens may be divided 
into two classes, froblems which may now already be solved on the 
basis of the known results of diffracticn theory are related to the 
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first class. AS an example of such a problem, one may point to the 
sroblem of diffraction of a plane wave by a frustum of a cone or by 

an infinitely long cylinder with a polygonal transverse cross section. 
Those problems whose solution requires obtaining (using the methods 

of mathematical diffraction theory) a whole series of new results 

nust be referred to the second class. In particular, in order to give 
a complete solution to the diffraction problem of a finite cone, it 

is necessary to have more precise knowledge on the diffraction laws 

of a semi-infinite cone. 


Summing up, one may say that physical diffraction theory aids 
one in analyzing and sorting out the diffraction phenomena for complex 
bodies, poses problems for treatment by mathematical diffraction 
theory, and enables one to effectively apply the rigorous results of 
mathematical diffraction theory for the solution of new problems. 


In conclusion, I express my deep thanks to L. A. Vaynshteyn for 
his valuable advice and regular discussion of the questions to which 
this book is devoted, and also for his attentive reading of the 
manuscript and for a number of useful remarks. I also take this 
opportunity to express sincere thanks to M. L. Levin for his interest 
in this work and his helpful remarks. 
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